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PROGRAM  BTLS  DESCRIPTION 


I.  INTRODUCTION 

This  is  a  description  of  the  computer  program  entitled  PROGRAM  BTLS  provided 
by  Science  Applications,  Inc.  (SAI)  of  McLean,  V a.  It  is  based  mainly  upon 
the  principles  discussed  in  an  article  by  C.W.  Spofford,  R.R.  Greene,  and 

J. B.  Hersey  of  SAI,  "The  Estimation  of  Geo-Acoustic  Ocean  Sediment  Parameters 
from  Measured  Bottom-Loss  Data,"  to  be  published  in  the  Journal  of  the 
Acoustical  Society  of  America,  under  contract  to  Naval  Ocean  Research  and 
Development  Activity,  contract  number  NORDA  SEAS  N00  1 4-78-C-021 1 . 

Individual  statements,  computations,  or  formulas  used  in  the  computer  program 
and  subroutines  are  discussed  and  analyzed  in  order  to  relate  them  more 
clearly  to  the  physics  and  assumptions  made  by  Spofford  et  al.  (1982).  This 
is  done  to  enable  the  user  to  more  readily  have  access  to  an  understanding  of 
the  program  and  its  proper  utilization.  It  may  also  enable  such  modifications 
(if  any)  as  may  be  desirable  to  be  made. 

The  geo-acoustical  model  upon  which  Program  BTLS  is  based  is  shown  in  Figure  1 
which  is  reproduced  below.  The  sediment  is  characterized  by  a  compressional - 
wave  sound-speed  profile  given  by  (equation  13  of  Spofford  et  al .  (1982)). 


where  gQ  is  the  gradient  at  the  top  of  sediment,  Z  is  the  depth  below  the 
top  of  the  sediment,  C  is  the  sound  speed  (with  C(o)=CQ),  and  B  is  a  para¬ 
meter  which  controls  the  shape  of  the  profile.  From  this  equation  an 
expression  for  the  sound  speed  gradient  as  a  function  of  sound  speed  C  (Z) 
is  found  to  be  (equation  14  of  Spofford  etal.(1982)) 


Figure  1.  Simplified  geo-acoustic  model  used  to  simulate 
bottom-loss  measurements. 
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From  Spofford  et  a1.(1982) 


These  two  equations  are  used  for  ray  calculations  involving  the  sediment 
and  are  of  major  importance  in  the  model. 

It  is  assumed  in  the  model  that  a  thin  layer  characterized  by  a  high  sound 
speed  overlies  the  sediment,  and  at  the  bottom  of  the  sediment  there  is  a 
basement  with  acoustic  properties  different  from  the  sediment. 

Various  possible  ray  paths  for  which  computations  are  done  by  Program  BTLS 
are  shown  in  Figure  2. 

Pages  showing  the  GEO-PARAMETER  INPUTS  and  INVERSION  PARAMETER  INPUTS  follow 
Figure  2. 

A  copy  of  PROGRAM  BTLS  is  provided  for  convenience  in  Appendix  A,  a  computer 
derived  flow  chart  is  in  Appendix  B,  and  Appendix  C  contains  a  computer  run 
stream  example.  The  computer  program  and  subroutines  are  discussed  in  the 
order  in  which  they  appear  in  the  printout. 


310NV  ONIZVUO 


Ray  paThs  (a)  and  angle-versus-range  curves  (_b) 

for  sediment  reflected  ( - -)  and  refracted/ 

basement -ref 1 ect ed  (— — )  paths. 

From  Spcfford  et.  al.(19S2) 


Figure  2. 
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TABLE  1  .  GEOPARAMETER  INPUTS 


INPUT  CODE 


PARAMETERS 


UNIT 


T 

NNFILE 
10  PT 

ZB 

ZS 

ZR 

CS 

Cl 


RATIO 

G 

BETA 

GAM 

REF 

AL  PH0 

RH2 

RH3 

D 


Two-way  travel  time  thickness  of  sediment  in 

tenths  of  second  0.1s 

Arbitrary  file  or  station  number 

Input  Option:  0  for  geo-parameters 

1  for  inversion  parameters 


Water  Depth  m 

Source  Depth  m 

Receiver  Depth  m 

Sound  speed  at  source  or  receiver, 

whichever  is  shallower  m/s 

Sound  speed  at  bottom  of  water  m/s 


GEO-ACOUSTIC  PARAMETERS  follow: 

Ratio  of  speed  in  sediment  to  speed  at  bottom 


o  f  wa  te  r 

Gradient  at  top  of  sediment  1/s 

Cononical  curve  type 

2 

Attenuation  profile  gradient  dB/m  /kHz 

Basement  reflection  coefficient 

Surface  sediment  attenuation  dB/m/kHz 

3 

Thin  layer  density  g/cm 

3 

Sediment  density  g/cm 

Thin  layer  thickness  in  meters  m 
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TABLE  2.  INVERSION  PARAMETER  INPUTS 


INPUT  CODE 

PARAMETER 

UNIT 

T 

Two-way  travel  time  thickness  of  sediment  in 
tenths  of  second 

0.1  s 

NNFILE 

Arbitrary  file  or  station  number 

IOPT 

Input  Option:  0  for  geo-parameters 

1  for  inversion  parameters 

ZB 

Water  depth 

m 

ZS 

Source  depth 

m 

ZR 

Receiver  depth 

m 

CS 

Sound  speed  at  source  or  receiver, 
whichever  is  shallower 

m/s 

Cl 

Sound  speed  at  bottom  of  water 

m/s 

DBLOSS 

FREQ 

A  value  of  loss,  frequency  and  angle  - 
example:  DBL0S=7 
frequency  =  200  Hz 
angle  =  20° 

dB 

Hertz 

degrees 

ALOSS0 

Low  frequency  loss  level  at  high  angles 

dB 

ALOSSM 

Minimum  loss  level  at  high  angles 

dB 

FM 

Frequency  of  minimum  loss  (of  ALOSSM) 

Hertz  (Hz) 

TIR 

Angle  of  total  internal  reflection 

degrees 

THC 

Angle  of  caustic 

degrees 

RATIO 

Ratio  of  speed  in  sediment  to  speed  at  bottom 
of  water 

- 

G 

Gradient  at  top  of  sediment 

1/s 

BETA 

Canonical  curve  type 

Sound  speed  profile  parameter 

- 

GAM 

Attenuation  profile  gradient 

dB/m2kHz 

REF 

Basement  Reflection  coefficient  in  pressure  units 

. 
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II.  THE  PROGRAM 


We  now  consider  PROGRAM  BTLS.  If  IOPT  =  1,  inversion  parameters  are  computed 
as  follows: 


OPTIONS  AND  ERROR  CHECKING 


found  by  use  of  Snell's  law, 
CosGj  =  Cos©2 


(3) 


Sediment  C3 

Figure  3.  Ray  Path  for  Total  Internal 
Reflection 

where  Ci  and  C2  are  speeds  of  sound  at  the  bottom  of  water  and  in  the  layer 
respectively  while  ©j  =  TIR  is  the  grazing  angle  in  the  water  at  the  water- 
layer  interface  and  ©2  =  0,  i_.e^  the  refracted  ray  is  parallel  to  the  water- 
layer  boundary.  Thus  cosTIR  =  Ci  and 

q" 

RATIO  =  C2  =  l/cos  TIR  _>  Cv)  if  nR 

is  expressed  in  radians  or  Cx 

C2  =  l/cos(TIR  x  ir  ^  ,  TTD  . 

^  '  180/  w  3  for  TIR  input 

in  degrees.  Then 

C2  »Cj  x  RATIO  C6) 


and  since  it  is  assumed  the  sound  goes  into  the  sediment  below  the  layer,  the 
assumption  is  made  that  the  speed  of  sound  in  the  sediment  below  the  layer  is 
the  same  as  in  the  layer  so  C3  =  C2. 


IF  GRADIENT  G  IS  NOT  SPECIFIED,  CALCULATE  IT  FROM  THE  CAUSTIC  ANGLE,  THC 

If  G  f  0  go  to  300 

If  THC  _<  0  let  THC  =  30  degrees.  This  is  a  starting  estimate  of  the 
angle  in  degrees. 

Call  GCOMP  (the  last  subroutine  in  the  program). 
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SUBROUTINE  GCOMP 


This  subroutine  is  used  to  compute  the  sound  speed  gradient,  gQ,  at  the  top 
of  the  sediment.  Quantities  are  defined  as  follows: 


Figure  4.  Ray  path  showing  angle  designations. 


THRC  =  30 
THC0  =  THC  =  0 

c 

B  =  Beta 

THRC  =  30  X  tt/180  converts  angle  to  radians. 

THINC  =  1  an  angle  increment  of  1  degree. 

D1  =  2*ZB-ZR-ZS  is  twice  the  effective  depth,  D,  given  by  equation  (10)  of 
Spofford  et  a  1 .  ( 1982 ) . 

TH3  =  ©=  cos  c_^  cos0  from  Snell’s  law  and  figure  4.  £7  J 

cl 


FACT  1  =  FI  =  2BC3 
1+B 


FACT  2  =  F2  =  1  +  J_  (1 
2B 


+ 


_ _ ) 

sin03cos03 


(S'  ) 
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In  order  to  obtain  the  equation  for  GZERO  =  gc  in  the  program,  we  need  the 
equation  for  the  range  as  given  by  equations  (16)  and  (21)  of  Spofford  et  al . 


(1982)  which  can  be  obtained  as  follows: 
SOURCE 


RECEIVER 


D1=2D 


Figure  5.  Diagram  showing  range,  depth,  and  apparent  bottom  grazing 
angle  relationship. 

.  "  _p__  _  =  -gL 

t-o.  v\  &io  ~  K/JL  ^  ^ 

where  D1  =  twice  the  effective  depth. 


( /o  ) 


For  sound  refracted  in  the  bottom,  an  increment,  Rs>  must  be  added  as  shown. 


Figure  6.  Range  increment  for  refraction  in  sediment. 
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S'  O 


In  the  sediment  tan  &  - 


on 

a?-) 


©I  Z  = 


7^*  <5>^> 


oo)  eT.  l/t)  of  S^Sf£orfLS±«l-  's 


k, --  a. 


•z. 


d  ^ 


w  ^  ere 


/  s 


,-fu 


e=>  r* 


the  turning  point  depth  for  the  refracted  ray,  or  the  basement  depth  for  rays 
which  reflect  from  the  basement.  Snell’s  law  applies  in  the  sediment  where 
the  sound  speed  gradient  is  assumed  to  be  vertical  so  (letting  subscript  o 
refer  to  the  top  of  the  sediment) 


f  PS  &(*)  ~ 


C  oS 


Cc 


C<*) 


(») 


and  taking  a 


derivative  with  respect  to  z  we  get 


?  i  * 


©<*) 


^  &(?-)  _  C  OJ~  &o  cX  (T  ^0  5  Qo  Q  ( z)  fl'-ZJ 

<M  ZL  C  o  c>t  2.  C  o 


where  g(z)  is  the  sound  speed  gradient  in  the  sediment. 

iav,  =  (-  c  ^Oo/Co) 

{a«9(Z-)=-  0(2.)  ’ 


Ccte(z)z  C(2.)^-fA'Z~ 

- IsLSu -  3- 

_  ~  g  ft;  H2 

K-  s.  j_\  ct&c^oiz-  -  a. 


cCzj  J  6<z^~  d?- 

9  ^  Zv)  d  2. 


In  the  model  it  is  assumed  g(z)  =  go ( 1 +B ) Co  and  we  have 

C ( Z )+BCo 

C<X)  -  C°  -  ^  o 

C  ©  -S  ©O 

&a)/ccn  + 

<3o(  n-*X*J 

'BCoCosBC*)  •  CCzJ  COS  6*1*1. 


c  cz.)  _  ^rP  _ 

g  f  z^)  “  c  o  i  ^e> 


(  \-h3jcosOo  9o  ( 


13  Co  £<=>? 


ecz.) 


.Co  cos-e^) 

4  — i - ; - “TT  SI/ 


g0(  i~T  bJc^sB^  <3c  ( t+ b)c 

C_a  fjLLeseiSl  + 


g0C'+?) 


C  OS  ©o 


C  <S>0 


-  >  //D 

f 'O 
r/7j 

dt) 

(/?) 

(zc) 

au 

(z  z  ) 

(L3) 
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Changing  variable  for  integration  from  z  to  @  gives 


Rc  = 

5  g„( 1+B ) co  s 


re 


B  cos9(Z)  +  cos_g (Zll  dG. 

'O  I  -  C0S^o 

e0 


7 


v  ) 


Interchanging  limits  changes  the  sign  to  give  upon  integration. 

X  C~  r <■  s;*XB 


Rs  -  X  - r - —  Ps  S  ?  H  ©  +  (-§- 

s  Jo  ( i+3Jc*s@°  L  '  ^ 


l)-L  J 

J  c  os  Sbj 


€>* 


&  -  & 


/  ") 


us 


Co 


fop  .  r\  ,  ^  -j  S'  *  9  c  os  ^>o 

-  2Bs"'fc)  +  +  -OJ  Sc 

“We' 


s  <30  (  i+B)cos&c 
Letting  9 ’  =  0  in  this  equation  for  a  path  refracted  in  the  sediment  without 
reaching  the  basement  gives 


(  ZL 


R ^ 


XT* 


9^  +do(  l-hB)cos&s>i  (9c 


Oc 


Si»&0(  ;+  2B )j  ■  <27-l 


Now  substitute  as  follows:  2D  =  D1 ,0^  So  ~  92  and  find  ^"0-^3 

D 1  .  C 

R  - 


f  ^3  4-  (  i  +  2B)7 

Ore.  gc  (  )  -hB)cos&j  cos&s  '  > 


K  “  ^  o{  _ 

R-  £k  •^Mfl  **  *  asS^D. 


I?  - 


D 


3o 

F 1  F  2-  Yr1^  . 

«c  6«> 


-F 


(3c; 


t&.\\9rc 

This  is  the  same  as  the  equation  for  R  in  the  program.  Take  the  derivative 
of  R  with  respect  to  0j  which  can  then  be  set  equal  to  zero  since  we  expect 
dR_  =  0  for  the  caustic  ray. 


d  Q 


<3°  \\^BX 


C 

o) 


Find  dGrc  as  follows: 
d90 


Snell's  law  is  cos^c  =  C  i  cos  Q? 

c3 

and  taking  a  derivative  gives 
-sin^rc  d&c  =  -Ci  sin^-a 

355-  cj 


d&rc  =  Ci  sin©3  and 

dBf  CjIT^c 


d  Fg  =  if  f 
dfc|  2B  /  sin^cosC^ 


+  ^jse^jtan^i 

sTn^T 


C5  3) 

(343 


By  use  of  identies  we  can  write 


dFo  =  1 

dfc^j  2Bsind^cos£% 


/  1  +  &  /sin<§  -Z0T&3 

'[  ^(rinfcv^r " 


Substitution  of  (2)  and  (3)  into  (1)  gives 

-p.  1  _  _  JR  zs  ?  C i  N  y  i*  ^  ^>3  ,  Ft  ^  _ 

_Vlc.se  Bee  (  c3)s,^  3 o  C  2.B s^G^cos  e>3 


-h  Ft 


-^Sec 


- 


o  . 


C*"1) 

C jy  ) 
Ui ) 


identities  gives  the  result 

w>> 


which  is  the  same  as  the  equation  in  the  program  for  the  sound  speed  gradient 
at  the  top  of  the  sediment  GZERO. 


The  program  follows  an  iterative  procedure  as  follows: 
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MAX  I T  =  100  is  maximum  number  of  iterations  beginning  with  computations  of 
g0  and  R  for  THRC  =  30.  Then  THRC  is  incremented  by  1  degree,  converted  to 
radians,  TH3  is  computed  from  Snell's  law  and  new  values  of  ^Q,  R,  and  THC2 
are  computed. 


If  |THC0-THCl|<_.01,  the  computation  is  finished  and  THC=THC2.  Otherwise 
iterations  continue  with  the  increment  changing  sign  and  being  divided  by  2 
until  convergence  occurs.  When  convergence  is  obtained,  the  output  is  printed: 
THC  =  THC2  and 
G  =  GZERO. 


Nonconvergence  is  indicated  by  the  output  statement 
ITERATIONS  AND  STILL  NO  SOLUTION. 

RETURN 

END  (of  Subroutine  GCOMP) 


Write  statements  start  printout 


CHANGE  DEGREES  TO  RADIANS 

THETA  1  This  is  done  for  these  angles. 

THC 

TIR  j 

CALL  ROUTINE  WHICH  CALCULATES  THE  REMAINING  GEOPHYSICAL  PARAMETERS  FROM  THE 
INVERSION  INPUTS 
CALL  GEOPHYS 

SUBROUTINE  GEOPHYS 


COMPUTE  DENSITIES  RH2 ,  RH3 

2. .  =  Pi 


(b't) 


7~)  -  -dll -  -  ——  =  quarter  wavelength  in  layer  at  top  of  sediment.  (^1  2  ) 

4  m  4 

A/ofe:  speeJi  o  f  souJ  i  *  *  U/Cr,  CZ  >  here 

•  ^  s  ~hea.c(  a  ^  s  /  «o  w  <*,  i~er ,  d  j ,  f  h  crMer  f  o 

f  orrecf-  /  e» y  e*~  ~t~  gts  23) 
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'  R  =  EXP(-ALOG(1O)*ALOSS0/2O  can  be  written 


(4^ 


R  *  exp(-lnlO)  ALOSS0 

20  .  We  can  see  this  as  follows: 


The  reflection  coefficient  R  =  (pressure  of  reflected  wave)/(pressure  of 
incident  wave).  Thus  ALOSS0  =  -20  log^R  "is  the  low  frequency  loss  levels 
at  high  angles. 


-ALOSS0  =  log _  R  =  ]n 

- 20  10  ft 


nlO 


In  R  =  - (lnlO)  ALOSS0 
20  . 


Then  e  raised  to  power  InR  gives  the  equation 


for  R. 


Zo  =  ( 1+R)/ ( l-R)Zj  .  This  is  the  impedance  Z3  of  a  medium  into  which  a  wave  (^£>) 


is  traveling  from  another  medium  which  has  impedance  Zj  and  is  given  by  eq. 

(2.14)  of  Brekhovskikh(1980)  with  subscripts  changed  to  correspond  to  our 

situation.  Normal  incidence  as  assumed. 

-(lnlO)  ALOSSM 
20 


R  =  e 

ALOSSM  is  the  minimum  loss  level  at  high  angles.  This  should  occur  when  the 
layer  thickness  D  =A2/4.  For  these  conditions  according  to  eq.  (3.23)  of 
Brekhovshikh(1980) , 


R  =  zf  -Z,Za 


(M*) 


If  ALOSSM  is  a  minimum,  then  R  is  a  maximum.  When  R  is  maximum,  then  Z^is 
greater  than  Zj  and  Z3  and  we  have  a  high  acoustical  impedance  layer  which 
acts  as  a  good  sound  reflector.  If  the  value  of  Z^is  between  those  of  If 
and  Zj  we  have  good  transmission  and  a  small  reflection  coefficient.  When 
Zz=V2/Z3  ;£=Oand  the  incident  ray  is  completely  transmitted.  In  this  ^ 
case  there  is  no  minimum  value  of  loss  at  a  particular  frequency  for  high 
angles.  Solve  the  equation  for  Z^. 
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Rzf  +  *Z/Z3  -  z|  'ZiZJ 
2l|  =  -Z,Z3(K->-i) 


-m 


(so) 

(SI) 

(*z) 


j-  iAj'tlly  the  secro-wW  expressi 


O'*. 


i  n  0.3  rer*^ei 

for  Zp  '  ^  +Vw's  s«l>r • 
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SPRAT  =  C3/C1  is  the  sound  speed  ratio. 

RH2  =  PP  =  Zo  is  the  density  in  the  layer. 

c2 

RH3  =P_  =  Z_3  is  the  density  in  the  sediment. 

u  Cj 


(s-i) 

O'4 ) 

(.SO 


COMPUTE  DISSIPATION  CONSTANT  ALPH9 

THETA  =  cof'.S-  c  OS  (THETA)  =  £^is  the  angle  in  the  layer  corresponding  to  (J>  L  ) 
the  bottom  grazing  angle  THETA  which  is  input  data.  The  sound  speed  Cj?  is 
assumed  to  be  the  speed  at  the  top  of  the  sediment  as  well  as  in  the  layer. 


(SI) 

CiSO 

(S<f) 


BS  =  sin 
BC  =  cos 

V32  =  Z3-Za\used  to  calculate  the  layer  reflection. 

z3+  Z^ 

V2/  =  Iz.-  Z  /[and  transmission  coefficients.  These  represent  equations  (31)  f 

za+  2iJ 

and  (32)  respectively  of$Dofford  et  al.(1982).. 

AUX  =  41TDF  sin  6*9  =  o  gives  twice  the  phase  length  for  propagation  of  a  wave  (<=>') 

~TT 

through  the  layer  as  is  shown  bel ow. ’■  C/  should 

in  \+r*r'i'«r  ke  U  S  f  ** P  '  ^  r  ' 

\Izt  +■  ^32  (  C  05  S  +  i  sih  <f) 


Vc 


Ctz) 


I  +■  ^2/  I62.  ( 

This  is  the  reflection  coefficient  for  a  plane  layer.  It  corresponds  to 
equation  (3.20)  of  Brekhovskhi k ( 1980)  with  appropriate  changes  in  subscripts 


giving  equation  (30)  of  Spofford  et  al.(1982)  and  use  of  the  relation 
Co -5  &  -t-isi'ncT  =  Q, ' 

/  fcZ,Z3  (~C2L,-rZz.)rz^-kZ3) 


T  = 


(U) 


[I  +  V2.1  V/52.  (c  o  S  £  +■  I  S  '  *1 
This  is  a  two-way  transmission  coefficient.  It  corresponds  to  the  use  of  the 

transmission  coefficient  equation  (3.21)  of  Brekhousk.hik(1980)  two  times; 

once  for  transmission  from  water  (1),  to  layer  (2),  to  sediment  (3),  and  once 

for  transmission  from  sediment  (3),  to  layer  (2),  to  water  (1). 
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For  a  wave  propagating  downward,  we  have  from  equation  (3.21)  of 
Brekhovskhik(198Q) 

1 


W  =  4Zi  Zo 

(Z  +Z2_)(Z2.+ZJ)  * 


-  i  kp  2.  d  i  kp  y  d 

e  +V32V2I e 


(6M  ) 


Taking  account  of  coordinate  axis  directions  and  layer  numbering  for  our  case 
we  should  have 


1 


W  ->  T.  =  4Z^Z?. 

1  (z“^7Tz^+z JT~  i«c  -i,c 

e  +Vj2_V^3  e 

and  for  a  wave  traveling  upward 


^-41^ 


1 


(Z,+Z2)(^+Z3)  -irf. 


icC 


((,(■  ) 


+vasv)2.e 

The  two-way  layer  transmission  coefficient  is  then  the  product 
T  =  T.  TZ  =  I6Z1Z2  /  Zfe 

M  (t 

X  1 


iot.  -  icA.\  /  -ioC  i<£\ 

[e  +VieV23e  e  +V32 V2/e  ) 

T  *  16ZfZ3  f  X2- 


(e 


(7-i  1 2-2 )(Z2+Z$)j 


T  = 
>where 


I  +  (e.72*  f  e“7^  Vu'j1 

i^iz3  (1 


_ -U^YZz^)) _  , 

(  +•  2^32  V^/  C  i>  s2«C  ^(V,32V/2))i 


(6>  9  ) 


sjle^s;»Q_  ,  (7o"> 


oC  -  =  ^2  d  si*  ©2  = 

the  phase  length  for  travel  through  the  layer,  and  we  see 

AUX  =  77"  s .  >1  g,  _  ^  <?.*  ^  i  r\  't  b  e  j  ) 

^ubro^d:  .  fj  C/  s-koni&{  be  /o  <^2,  here. 

The  denominator  of  T  given  in  the  subroutine  is  the  same  as  in  the  last 

expression  for  T.  Sign  changes  of  the  Z's  due  to  changes  of  direction  which 
might  be  troublesome  cancel  out. 
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AUX  = 


C?2) 


10 


-P  bLPSS/,0 


-\r 


T 

als  =  —  10  /og/0  |- 


/O 


-mLos.s/w  .  ,2 


-  l/‘ 


-  ~DL  . 


T 


We  can  arrive  at  this  relation  by  starting  with  the  assumption  that  losses 
behave  linearly  in  sediments  as  expressed  by  equation  (13)  of  Spofford  et  al  . 
(1982) 

L  =  K  •  S  •  f^  where  K  is  the  absorption  coefficient,  S  is  path  length, 

and  f  is  frequency, 
k 


This  can  be  rewritten  in  integral  form.  From  Snell’s  law  and  differentiation 

we  fi nd  \ 

tosecz.)  -  cos  e*  T  ] 

C  c 

-  S  ( *  OCz)  q(2i)  r  3  (z.)=  <1^1  >  O*  ' 

VL-  -2 f2 C o f k j ko t kzZJ. d 3 (z.)s  ffv(  I  +  B)<.q  (7c 

Jo  cosBe^Cz.)  C  (z.)  -h  E  Co 

.  dl'  fff^^X(K^kJzXcCz)+BCo)^^’  07 - 

dl--  f^r^sik‘c^k>)c  cz)i3^  +k*Bc°*y9° 

The  sound  speed  profile  for  the  model  is  assumed  to  be,  as  given  by  equation 


(13)  of  Spofford  et  a  1 . ( 1 982 ) :  / 

CCz.)- /2g0  (l+3)CvZ  +  (  t+sfcl]  -SCc, 

from  which  _  o  o 

(C+ZC0f  2$e(  l+B)C»Z  *  (l+B)  Co 

(c+lC$-(\+3?Cl  =,  25o  (  )+B)C.Z  , 

2.  -  j  c  *  BCoi2-  (  I  isfc}-  _  CC  *  ZC^~-  (  KB)2- cl- 
_  2.  9o  C  l  +  B)Co  ~  '2-°~ 

where 

0-  -  3o  (  1+  B)  Co 


(7*t; 

(  *0 
(SO 

^3; 
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^  ^  Zm-h  pi**  .'fr'Q*  Pvy~  ^  L-  3'^®^ 

re! 

2>LzZ3=p£ - _ 

<3 oC  l+B )czx&i 


0  I 

\  ft~C*£os  &  ,  y°  Cr,  tA? g|/Co  CosO  n^j  _/ 1  j.-b)  i 
2<rc«®.  |[t«A  +  0<j  (  I  +  O/Loj 


+  2^C„  f  BC^Lf/T^  (  1+sfCo'j]  <*e 


2.0- 


(8*0 


Let  D-  -_*Jk 


T)L  -  D 


5o  (  \t  Z)ooS&e> 

r&'t 


~J 


kc>CoC-os  S>  ,  ko'Cc>cos@fc<>  coj’Q  2.BC&cosQ  v  ^2^r 
CO,  &o  2^-cos  oJ^Wo  rr,  &  — 


c?  -  J^rlj  v-  2 
zi£c~etgdr- 


1)  )  CD  ^  {^2LS*S-+  c£-  f  +2£&£l^~l/i+2s)£osgj 

])  L  =  C  O-J  \  j  2_<Tcos£)0\_  ccs2^0  O  <o<,G>o  1  ^  J 

J8v^ 

fW- -^j^®  (u 


DL  -  £<3 


ko£;«@  ,  K*cl_  f  1  fs;»&-  tS’* 2jg 

g>0  2(rc0s@4(ci>ss&(  ^  /  ^-os^(2  4  y 


Cc»s 


-(  l  +  2  s)  3  i  *  S’  +  3^<>  & 


<£>' 


•f 


%  (f £  jcJ  _ 

2.0“  IcoS^Oc 


f-  -  ^ +  J^s;^>28>)M- 
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Interchanging  the  limits  of  integration  and  factoring  gives 

DL  =  ^  ?  i yj^).  &  1  — 

‘  '  ‘  (  3  /2<3~Ct 


<%z>  (  l  +  d)coS&o 


2o~Cds3<9< 


JO  ^  -  - -  «—  ^ 

/  @  +  si  H&ces&d 2  BC?  ko  +  BCe  ^Cd 

+  (  2  A  z<r  co^e0  , 

/  . .  e)(Z£k°  -  ( i ±2S] kJ 42  ^Co  <cJ ) 

(  A  2<rc.0ie>0  / 

<•(  e)(W  gc|wr»2»l'1  • 


) 


This  is  equation  (23)  of  Spofford  et  a1.(1982).  Note  that  in  the  second  line 

of  equation  (23)  of  Spofford  et  al .  (1982)  the  power  of  cos  So  should  be  2 

instead  of  4.  Now  consider  figure  8. 

jB  *  VzZ . 

.o-rer- 


r,=  t*x. 

Fi'^ca.  re  Sr.  ■T  v*\  At gT«  /" ’ 

r-e  ■prevc-A 


Se<W..'iAAe»^  +  _ Fe.Fr ac-fe^Z  K  i  n  ie^.'^enT 

7ZZsk]ps  "pef“  r  e  P\e  c.i-zc>{  a.  v*-<zl 


Loss  due  to  attenuation  in  the  sediment  only  is 


DL  =  -10  log 


X2. 


/O  X 


/ 


Cs?) 


^  j_2.  where  Tj  and  are  transmission  coefficients  defined 


We  see  X,  =»  l^2Z~o  ^ 
previously.  The  intensity  of  the  reflected  ray 

7 

I^_  =  V*-  I D  where  V  is  the  layer  reflection  coefficient  defined 
previously.  The  intensity  of  the  sound  leaving  the  bottom  after  reflection 


from  and  transmission  through  the  layer  is 

x*  =  z,  +  , 

Isr  ■  TfU  i-VzIc 


To 

An  expression  for  is  needed, 


(9CJ 
(i 0 

OA 
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Now  /  0  /eg  -  10  /03T,  “  >  O) 

lo0  Tz  -  lo3T,  -  2> L/JO,  W 

■2^2.  —  ei^h'lo^  (I03  jTy  —  T>Ly/j  o) ,  (&) 

r^S  4f-'  VZ<-  (T£/I'>)a„i;i03(l^Tl-T>L//oX  CIO 

V2  ^  tilo3  ( 1 0^1,  -T>l//d).  M 


To 

The  measured  loss 


rea  loss  1 

p&LoSS  =  —  )0  Jt>3 


Is 

To 


cn\ 


—  -DB IPSS//0  =  lo$[V2  +  ("  T^/J0ja  V) f//aj (7iy  77^  ~X>L//djJ(jr 


Note  that 

10 


-VBLOSS//P  _  \f+(r%/l^CL*ti'lcd  (/oqT^Io  -  n/jo)  (/<*>) 

(j'l/lD)a.^Ilo^(}o^lflo  -It/10)  r/au 


l°3 


-V?J-oes//o_  y2} 


)  =  /^772-r=  - 


6^2-; 


~T‘ 

=_  3Z.//0  ■  T,rz  (*>*) 

idvujSSTjz^z'  =  >9is  C 

D 


'Vow  I)L  -  -  /  O  /«3,„  / 

as  given  by  the  expression  in  the  program.  /O  o1~-q  This  cxe^ow  — 

\s\t\i~or~  S  be.  c  U.et^\g^W  ^  7^  u/  /  /■  K  i~V\e 

|7re  ceolr^  €.^h*i^7ow 

Equation  (23)  of  Spofford  et  al.(1982)  for  the  dissipation  loss 
in  the  sediment  can  be  rewritten  in  terms  of  frequency  FK  in  kHz,  SIG 
AMU  =^ANU  =P,  ALPH0  =  ^o-ko,  and  GAM  as  expressed  in  the  program  as 


,  ,£S. _ _ 

>K  <3>(  l+B)  COS&2. 


zd~ r - 


[5  >  <n 


/  ^ 


f^M^jle*hS^Qzc°s% 

\~Z *=T°s  ©2  7  2<r  £ 


r/  i  s  f‘  H  ^  4-  ft  43  )  E  ^2. _ 

*‘l -**&*.  0eVi7 ’ 


( ><*>7 


2) 


wh  ere@-0,  4*3  the  grazing  angle  of  the  ray  in  the  layer  with  the  top  of 
the  sediment,  ko  the  gradient  of  the  absorption  coefficient  and  °£o—  ^ o . 


which  is  the  equation  for  the  dissipation  constant  at  the  top  of  the  sediment. 
For  this  computation,  it  is  assumed  that  the  speed,  C £_»  and  angle,  ©2  are 
the  same  in  the  layer  and  at  the  top  of  the  sediment.  Also,  note  that  Fk  is 
the  center  frequency  of  a  one-octave  frequency  band  used  for  processing 
measured  data. 

RETURN 

END  (of  SUBROUTINE  GEOPHYS) 

GO  TO  499. 

499  CALL  CANGLE 

SUBROUTINE  CANGLE 

D1=2ZB-ZR-ZS  is  twice  the  effective  depth,  "O’,  given  by  equation  (10)  of 
Spofford  et  al .  ( 1 982 ) . 

CALCULATE  SEDIMENT  THICKNESS 
B  =  beta 


CZP  =  2000 

Sound  speed  at  the  top  of  the  sediment  is  assumed  to  be  •  The  equation 
for  one-way  travel  time  in  the  sediment  is  needed.  The  necessary  equation 
can  be  derived  by  assuming  the  sound  is  propagating  downward  in  the  Z  direction. 

J  '2__ 

By  definition  C(2.)j  'Zit~  ^ e  Q/cl ) 


ait  = 

t  - 


I 


CCZ.) 
dL  ~2- 

C(z-) 


(/a) 

OoCt) 
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The  sound  speed  gradient  in  the  sediment  is  given  by  the  model  as 

a  eU,  _  go  7  !±B)Cz 

Let  <r  i  gof  i+  8^2.  .  .  . 

J  2  -  ±  ^  ^  8Cz)  a  Q 

d  Vod+BXi  <r 

‘ccz)  c±-k£x - JLc , 

c  ncu) 


["* J> 
o/j 
(">-) 
("*) 

(/>*) 

(,,r) 

Cnc) 
("V 

g.o+B>U"^  V  '  J-  (//p) 

This  is  the  one-way  travel  time  for  sound  in  the  sediment  and  is  equation  (3) 
in  "DOCUMENTATION  OF  BOTTOM  LOSS  UPGRADE  PARAMETERS."  This  equation  cannot 
be  solved  in  closed  form  for  depth  of  sediment  BZ=z  so  we  will  find  C(z)  at 
the  bottom  of  the  sediment  (basement  depth)  by  use  of  Newton's  method  of 
approximating  a  root  and  use  C(z)  to  find  z.  From  the  last  equation 

■  <3*  ( =■  -£^-7  +  B  In  CCz)  -Bln  C2  ,  0/9) 

Cz<3o(  )+%)t  -  Cte)  ~  CZ  0ZC>) 

Ct 2.  o  O^B)l t  -h  I  —  B  )nC 2.3  — '  C  ( z.)  j-  T&  Cg'liq  C.  (z~)  -  f  ( i- 1) 


t  -  ( 

^2,i)  (  1 4- s  +  B q,  <Uc][ 

j.  *  X  r’ffej  4-  24  li  O'*)  -  ^  -  B  C2  In  CzJ^ 

t^±[C(*-)~Cz  ^C2ln-^r]  , 

t-p[c-g-l 


(  (iO 


Thus,  in  the  subroutine  we  find 

)j3  e7  if  3.) zf  ^  )n C^J  as  on  the  left  side  of  the  previous  equation. 
The  formula  for  Newton's  method  is 


A.  =  A  -  f(A)  where /J  is  an  approximation  and  t)  >  is  an 

1  VTK) 

improved  (  |25) 

Let  A  =  CZP  and  A|  =  CZ,  then 

CZ  =  CZP  -  f(CZP)  where 

f'(czp; 

(i^) 

f(CZP)  =  A-CZP-BC2  In  CZP  and 

(jif ) 

f*(CZP)  =  -1  -  BC2.  gives  the  result 

CZP  ' 

0^) 
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CZ  =  CZP  +  (A-CZP-BCalnCZP) 

1  +  BC^ 

CZP 


(/27) 


as  in  the  subroutine.  Starting  with  CZP  =  2000  as  a  first  approximation, 
iterations  continue  until  convergence  is  completed.  This  is  chosen  to  occur 
(as  in  the  subroutine) 

IF|CZ-CZP|<.01  ,  let 
CZP  =  CZ  =  C(z). 


From  the  formula  g  =  dC  =  gn(l-t-B)C2 

d.  z  C  +  BC^ 

we  find  Z  by  integration  to  be 
Z  .jdz 

rCz  C4g c, — 

Z  J<;,  9  o  ( 1  's)  <-2 

'CZ. 

I  I 

2  =  ^o()rB)Cz 

I 


J 


CC4%C^)c\C, 
r^Z  iCZ. 


0**) 

(n’t) 

(  I 

('7') 


z  •  I  £-  +-2CzcJ^ 


H-B Xz 


(  **0 


(t  *y) 


,L  *-  J  > 

B^ez  4 B^cf  -  c£-2B 


70  Bz  =  z  s%(  n-s)Ct.^ 

as  in  the  subroutine,  gives  the  sediment  thickness. 


Vz.  i-  ZCzf-  Cl(  1  +  zT]^ 


( (  3  s) 
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COMPUTE  REFINED  CRITICAL  ANGLE  THRC  = 

CALL  TOTCRIT 
SUBROUTINE  TOTCRIT 
Check  that  C^C^ 

This  insures  that  if  total  internal  reflection  occurs,  it  will  be  from  the  top 
of  the  (high  speed)  layer  and  the  layer  is  distinct  from  the  sediment  beneath. 
IF  (C35-C2)  GO  TO  10  (where  angle  is  computed) 

Print  format  "C3  MUST  BE  LESS  THAN  OR  EQUAL  TO  C2" 

COMPUTE  TIR  =  ANGLE  OF  TOTAL  INTERNAL  REFLECTION 
TIR  =  6 

If  C2.£.CI  GO  TO  100.  In  this  case  TIR  is  not  calculated.  For  total  internal 
reflection,  it  is  required  that  CZ>CI . 

TIR  =  CoS"1  (l^c) 

from  Snell's  law  when  refracted  ray  angle  =  0. 

100  CONTINUE. 

COMPUTE  THC  =  CRITICAL  ANGLE  (CAUSTIC). 

This  is  done  by  an  iterative  procedure.  Start  with  definitions  as  follows: 


PI2  =  T r/z. 


(  I  37) 


D1  =  2  2  3'Z^-Z  S  - 
< \~yJL  re  ce  l^er 


DELT  =  TT/ 


D2  =  2  C3  S 


go  C  \fB) 

FAC  = 


(11O 

(i'll-) 


as  in  figure  9. 


S«. 
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.  For  the  first  approximation  use  9rc 


An  approximate  range  is  found  by  using  9rc  and  9^  in  the  range  equation  found 
by  combining  equations  (16)  and  (21)  of  Spofford  et  al.(1982). 

(H3) 

The  range  is  varied  by  incrementing  the  angle  9rc  by  del  t  -  radian: 

AUX  =  9rc  +  IT 
too 

It  is  assumed  the  incremented  angle  AUX  is  between  anc*  the  angle  of  total 
internal  reflection. 

IF  (AUX  >%)  GO  TO  50.0. 

IF  (AUX  TIR)  GO  TO  500. 

The  angle  9j  is  recomputed  in  terms  of  the  incremented  value  of  Grc  and  a 
.new  range  RAUX  is  computed  by  the  same  formula  that  gives  R  above.  If  the 
new  range,  RAUX,  is  less  than  the  old  one,  R,  the  procedure  is  repeated  until 
the  new  one,  RAUX,  becomes  greater  than  the  old  one,  R.  When  this  occurs, 
the  increment  delt  is  change  to  -  slslt  and  the  process  repeats  until  conver- 

3 -4 

gence  occurs,  i.e.  IF  |  delt  <  1  X  10  .  The  idea  of  this  iterative  process 
can  be  illustrated  by  figure  10  shown  here  adapted  from  figure  8.b.  of 
Spofford  et  al . ( 1 982 ) . 

B 

vwi'h  Voiluf »  R 

F,  10 ■  X  I  I  usll "r-A  tit>~  e  "P  »  v->  {  <-~v  VAILS’  O^C~  j  •e  ^  ~t  t  L  6? 

C  Aw.6if-;c  ^ 


R  =  D1  +  2BC3tan^  J  1  +  (1  +  ^  ")  1 

&rc  <gD(  1 +-3)  \  si^e^cpsGj  ZB 

*-  — 

Here  the  lower  limit  in  tVie  ion  is  zero, 

R  =  D1  +  Da  tan  03  fl  +  f  1  +  Bj  \  1  J 

~tcL*&rc.  /_  \  J 

MINIMIZE  EXPRESSION  FOR  RANGE  BY  CONJUGATE  DIRECTIONS. 
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The  angle  0£  is  desired.  We  find  the  angle  by  varying  the  range  R  as  a 
function  of  angle  0  until  R  is  a  minimum.  When  this  occurs,  0=©c  the  caustic 
angle.  If  convergence  does  not  occur,  the  procedure  is  stopped,  and  the 
statement  "CRITICAL  ANGLE  FAILED  TO  CONVERGE"  is  printed. 

BOTTOM  ANGLE  OF  REFRACTED  RAY  AT  CRITICAL  ANGLE 

THC  =  tan  Dl_  =  0  is  found  from  the  range  and  twice  the  effective  depth  (W&) 
R 

as  shown  in  figure  11.  This  is  the  observed  angle  of  the  caustic  ray  as  illus- 


The  speed  at  the  greatest  depth  in  the  sediment  reached  by  the  basement  grazing 
ray  found  from  the  formula  for  the  sound  speed  profile  as  a  function  of  the 
depth  assumed  by  the  model  is 


1+B)C3BZ  +  c23  (1+B) 


2 


-bc3 


(H7) 


THIB  =  cos'  Cl  =  ©iR,  the  grazing  angle  in  the  water  for  the  basement  grazing  //y'g) 
ray  is  found  from  Snell’s  law. 

RETURN 

END  (of  SUBROUTINE  TOTCRIT) 

We  now  continue  with  SUBROUTINE  CANGLE 
COMPUTE  CRITICAL  RAY  DEPTH  CRD 

From  Snell's  law,  the  speed  at  the  greatest  depth  of  the  caustic  ray  is 

CAUX  »  Ci  (  J 
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We  previously  found  the  equation  for  BZ  the  sediment  thickness  as  a  function 
of  CZ  the  sound  speed  at  the  corresponding  depth.  Now  let's  use  the  same 

r 

equation  with  BZ  =  CRD,  CZ  =  CAUX  and  -  Cj  the  sound  speed  at  the  top 

of  the  sediment  so  CRD  =  2.3 e(  TTbJc^  [f  ^  (/$£) 

as  given  in  SUBROUTINE  CANGLE  and  illustrated  in  figure  11. 

CALCULATE  VELOCITY,  VELOCITY  GRADIENT,  AND  ATTENUATION  CONSTANT  AT  DEPTH 
OF  500  METERS. 


Setting  depth  in  sediment  Z  =  500  meters  and  using  the  formula  for  the  assumed 
sound  speed  profile  we  find  the  speed  in  the  sediment  at  depth  500  meters  to  be 

BSPD  =  [2  (  /+b)(<t iX&o)  +  cf  ( 1  -+  b)2J^  -  (,sl) 

where  Cj  is  sound  speed  at  the  top  of  the  sediment.  The  sound  speed  gradient 
at  500  meters  depth  in  the  sediment  is,  from  the  assumed  sound  speed  gradient 


equation, 

<3  o  C  l’b'3  )  C3 

BGRAD  =  BSPD  f 

The  sound  absorption  coefficient  is  found  to  be 
ALB0T  =  K  =  KQ  +  k^  Z 
=  °Co  +  V(500) 

COMPUTE  TWO  WAY  TRAVEL  TIME  IN  SECONDS 
TWTT  =  2T 


( isx) 


SPRAT  =  C1/C1 

COMPUTE  BASEMENT  CRITICAL  ANGLES  THIB,  TH3B 
Sound  speed  in  the  sediment  at  the  basement  level  is  found  from  the  sound  speed 
profile  equation 

CB-  (dz)  f  (  H~B)ZCj'  -££3  ( /&) 
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£©uirce 


R/2. 


Snell 's  law  gives 

COS  Ox'S,  - 


£j, 

c9 


For  a  refracted  ray  with  its  vertex  at  the  basement  as  shown  in  figure  12, 

(/St) 

&ie.=  cos'1  (/sc) 

which  is  the  largest  bottom  grazing  angle  for  a  ray  which  is  refracted  above 
the  basement. 


-I  Cl 


®3S  “  CoS  1 

In  terms  of  these  angles  the  range  equation  becomes 

7)1  +•  I  +  /  '  &3B 


(/si) 


R  "  ta~e,B  ^  9,  (  i  +  b) 

as  in  SUBROUTINE  CANGLE 


2 .3 


03Scos  ez8/ 

COMPUTE  SMALLEST  APPARENT  ANGLE  FOR  BASEMENT  REFLECTION. 


r 


/ss) 


(/si) 


From  the  value  of  R  found  last 

&s= 

as  shown  in  figure  12.  If  the  angle  is  larger  than  this,  there  will  be  a 
reflection  from  the  basement. 

CHANGE  RADIANS  TO  DEGREES 

>< 

ft  =  TIK  X  ,  xo/v  --  «“3'e  ^  t'KI  I'*'™1  refle" 

n  Q  X  \&o/rr  r  *pr*re*t  c~~s1-rc. 

h'  “  x  f^/rr  Ut£.> 

'  ^  x  |  p/n-,  t -j  '  «3'£  1 11  tl,£  r 

/?3=C7lB  '  ref|*e.t,'®"  frc £t,e 
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PRINT  OUT  CRITICAL  ANGLES 
Format  "CRITICAL  ANGLES"  -  -  - 
RETURN 

END  (of  subroutine  CANGLE).  We  go  back  to  main  program  and 

PRINT  GEO-PARAMETERS  in  Format  "GEO-PARAMETERS"  - 

CALL  ROUTINES  FOR  FREQUENCY— INDEPENDENT  INTERMEDIATE  VARIABLE 
The  subroutines  REFL,  REFR,  and  BASE  are  used  to  compute  quantities  needed 
in  the  bottom  loss  subroutine, LOSS.  CALL  REFL 
SUBROUTIKE  REFL 

NSAP  =  smallest  angle  to  reach  surface  due  to  sound  profile. 

COMPUTE  NSAP 


NSAP  =  1  is  defined. 
IF  (CL>CS)  go  to  100. 


For  a  ray  path  as  shown  in  figure  1 3, 


S  eurtt 


(orr«ceiVerJ 


C 


Fi‘e«rr  3-  /  ^strati'0,  *  °  l[  ^  1  'A<J 

/  e  SAP. 

a_ray  may  be  horizontal  at  the  source  and  have  a  bottom  grazing  angle  SAP 

where  Snell’s  law  gives 
^5  _  C  / 

/  c«s  5 A'  ' 

SAP  =  Cos  1  If  C  /  Cs  this  calculation  will  be  skipped— the  angle  'K*} 

is  unrealistic. 


/Co  ^ 


NSAP  =  SAP  +  1  is  redefined. 

SAP  is  printed  in  format  "SMALLEST  ANGLE  TO  REACH  SURFACE."  - 

COMPUTE  REFLECTION  COEFFICIENTS  FOR  ANGLES  GREATER  THAN  TIR,  SAP 
NTIR  =  TIR  +  1  angle  in  degrees. 

N  =  MAX9  (NTIR, NSAP)  says  pick  the  largest  of  the  angles  expressed  as  integers 


in  degrees. 
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These  are  constants  to  be  used  in  computing  impedances 


12=02^2 


needed  to  find  reflection  coefficients. 


^3=P  3 C 3  ) 

DELT  =  n/]80  is  used  to  convert  angle  in  degrees  to  radians  and  as  an  increment 
in  angle  for  a  one  degree  change. 

TH  =  (NSAP-1  Wi80=  9  is  angle  equal  to  SAP. 

FAC  2  =  C2/C1 

FAC  3  =  C 3 / C 1  >  These  are  speed  ratios. 

FAC  S  =  Cs/C]  J 

AR  =  2ZR/CS  =  twice  travel  time  from  receiver  to  surface. 

AS  =  2ZS/CS  =  twice  travel  time  from  source  to  surface. 

D1  =  2D  =  twice  average  distance  from  source  arvd  receiver  to  bottom. 

DO  200  I  =  NSAP,  90  and  IF  (l<N)  GO  TO  200  insure  the  following  calculations 
are  performed  for  angles  varying  by  steps  of  1  degree: 

For  angles  greater  than  the  larger  of  TIR,SAP  Surface  Scattering  Travel  Time 
Differences,  SIN2,  Reflection  Coefficients,  Travel  Time,  and  Geometric  Inten¬ 
sity  are  computed.  If  SAP<TIR,  for  each  step  from  SAP  to  TIR  Surface 
Scattering  Travel  Time  Differences  only  are  computed  because  total  internal 
reflection  occurs  in  this  range  of  angles,  so  the  reflection  coefficient  is 
unity  and  special  calculations  of  the  other  quantities  are  not  applicable. 

For  each  step  from  TIR  to  90  all  quantities  are  calculated. 

SURFACE  SCATTERING  TRAVEL  TIME  DIFFERENCES 

The  sine  of  the  angle  made  by  the  ray  at  the  source  is  found  from  Snell's  law 


where  C}is  the  bottom  grazing  angle. 
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Travel  time  difference  can  be  found  as  illustrated  by  the  following  figure. 


r  + 


?c cittereM  p^tk  cK'ffene^cre  , 

Two  different  paths  are  designated  K,1  and  K,2.  The  assumption  is  made  that 

the  bottom  grazing  angle  is  approximately  the  same  for  the  two  paths.  Then 
the  surface  reflected  path,  K,2,  will  be  longer  than  the  path  K,1  by  approxi¬ 
mately  2ZS  sin0^  and  the  travel  time  difference  is 


(IOU) 


Interchanging  source  and  receiver  in  Figure  14  yields 


RSI  N2  ( I )  =  sin6^>the  grazing  angle  in  the  high  speed  layer  by  a  procedure 
similar  to  the  one  which  yields  S  above. 


REFLECTION  COEFFICIENTS 


CIO  1) 


These  are  the  acoustic  impedances 
of  the  water,  layer,  and  sediment.  ( 


icients  for  the  boundaries  designated. 


32 


RGMI)  .  O1^ 

is  the  factor  for  spherical  spreading  which  represents  decrease  of  amplitude 
with  range. 

RETURN 

END  (of  SUBROUTINE  REFL)  (Go  back  to  main  program.) 

CALL  REFR 
SUBROUTINE  REFR 

NTIR  =  TIR+1  gives  angle  in  degrees. 

NTHC  =  THC+1  Caustic  angle  +  1  degree. 

N  =  NTHC  -  NTIR  integer  giving  difference  between  the  caustic  angle  and 
angle  of  total  internal  reflection  in  degrees. 

DEL  =TT/180  is  an  increment  of  1  degree  expressed  in  radians. 

TH  is  the  caustic  grazing  angl e  expressed  in  radians  =  9rc. 

Z,1 

12  f  factors  used  to  compute  reflection  and  transmission  coefficients. 

) 

■  FAC2 

FAC3  /  sound  speed  ratios 
FACS  J 

SEST  =  Ore) 

/Starting  estimates  of  small  and  big  angles. 

BEST  =  Ore) 

D1  through  D14  are  factors  to  be  used  in  computation  of  geometric  amplitude 
and  other  quantities  to  follow.  They  are  expressed  in  terms  of  the  same 
quantities  used  before,  e.g.  in  SUBROUTINE  GCOMP  and  GEOPHYS. 

DO  10  1=1,90  sets  the  array  for  transmission  coefficients  =  0. 

DO  100  1=1, N  and 

J  =  NTHC-1  insures  computation  is  done  for  each  integer  angle  between  caustic 
angle  and  angle  of  total  internal  reflection. 
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R  =  Dl/tanTH  gives  the  range  not  including  that  portion  required  for  refraction 
in  the  bottom. 

CALL  REFRANG 


SUBROUTINE  REFRANG 

IF  TH>TIR  go  to  100.  Otherwise  print 

"SPECULAR  ANGLE  OF  REFRACTED  RAY  LESS  THAN  TIR" 
and  stop  computation. 

D1  =  2D 
PI2  =  TT/2 

D2  =  2C3B/g0(  I+-B) 

D3  -  1 

2B 

FAC3  =  Cj/C , 


-  ©3 


R  =  Dl/tanTH 

T3  =  C-&S  ^  ^  (sE5t)  g-jves  (by  Snell's  law)  angle  in  sediment  for  ^  1  *  J 

bottom  grazing  angle  SEST.  _  (  ni 

23. -  ^  U«eJ  l  +-' 


VAL 


2.3 


I  +77 


B: 


S  ih  E>2  c  os  6% 


~t<xv\  s^^>t  (  h-b') 

VAL  =  R1  -  R  as  shown  in  Figure  15,  the  difference  between  a  computed  range 


-R 


10.0001 
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The  angles  SA  and  BA  correspond  to  rays  a  and  b  for  range  R2  shown  in  Figure 

2. 

COMPUTE  SA 

DO  and  IF  statements  start  an  iterative  procedure  which  varies  SEST  and  ©3 
by  +  or  -  increments  until  | DEI_/ <1X10  while  insuring  TIR£SA*$"C.  When 
convergence  is  completed,  the  error  VAL  between  R*  and  R  will  be  very  small 
and  the  smallest  bottom  grazing  angle  SA  for  which  a  ray  refracted  in  the 
sediment  can  travel  from  source  to  receiver  is  found.  If  convergence  does 
not  occur  after  500  iterations,  the  statement  is  printed  "SMALL  BOTTOM 
REFRACTION  ANGLE  DOES  NOT  CONVERGE,"  and  the  procedure  stops. 

SA3=T3  gives  the  angle  of  the  ray  at  the  top  of  the  sediment  corresponding 
to  the  bottom  grazing  angle  SA. 

BA  INITIALIZATION 

T3  =  63  =  cos^C3  cos  BEST  is  the  angle  at  the  top  of  the  sediment  for  bottom 

Ci 

grazing  angle  BEST  (with  initial  value  Grc). 

VAL  has  the  same  significance  as  in  COMPUTE  SA  with  SEST  replaced  by  BEST. 

■  COMPUTE  BA 

DO  and  IF  statements  generate  an  iterative  procedure  which  varies  BEST  and  ©3 
by  +  or  -  increments  as  needed  until  /dEl/<  I X 10 
while  insuring  SAH  ' 

The  procedure  is  similar  to  the  one  for  COMPUTE  SA.  Upon  convergence  the 
biggest  bottom  grazing  angle  BA  for  which  a  ray  refracted  in  the  sediment 
can  travel  from  source  to  receiver  is  found.  If  convergence  does  not  occur 
after  500  iterations,  the  statement  is  printed  "LARGE  BOTTOM  REFRACTION  ANGLE 
DOES  NOT  CONVERGE,"  and  computation  stops. 

BA3=T3=  G^gi ves  the  angle  of  the  ray  at  the  top  of  the  sediment  corresponding 
to  the  bottom  grazing  angle  BA. 

RETURN 

END  (of  subroutine  REFRANG).  Now  we  go  back  to  SUBROUTINE  REFR  . 
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IF  ( B A>TH IB)  go  to  998.  THIB=cos,CjL  so  for  BA^THIB  the  ray  hits  the  basement.^ 

C8 

Otherwise,  compute  the  sines  of  grazing  angles  BA  and  corresponding  angles  for 
the  same  ray  in  the  layer  and  sediment. 

0740 


BC1  =  cosBA 
BS1  =  sinBA 


■F 


c  os 


BC2  =  cosBA  (Snell's  law) 

C| _ _ 

BS2  =  \j  1-(BC2)^  =  sin  ^2  (1n  layer) 

BC3  =  SCI  -  C  os  -  Cos  <9j 

<-/  C, 

BS3  =  sin  £2j(in  sediment) 

BIG  ANGLE  REFLECTION  COEFFICIENT 
A  i  =  Z,/s^bJ) 

t<2_  =  ^2./s '  ^^fAcoustic  impedances  used  to  compute  reflection  coefficients.  ({£ t} 
A3  •2$/s\'*6i) 

BV32( J)  =  ( A3-A2)/ ( A3+A2) 

BV21 ( J)  =  ( A2-A1 ) ( A2+A1 ) 


0?$) 

(176) 

C-t  77  ) 

(  I  7  Sr) 

C 

( (to) 


f(^) 

Os^) 

Big  angle  reflection  coefficients  for  bouncries 
designated  by  the  subscripts.  f/S^) 


Big  angle  transmission  coefficient  factor  (two-way) 


BTRF  =  J**-  ~f~ - — 


)J 


(iss) 


similar  to  the  one  in  Subroutine  GEOPHVS. 

BSIN2(J)  =  sin  (in  layer)  ‘ 

GEOMETRIC  AMPLITUDE 

The  geometric  spreading  loss  factor,  equation  (38)  of  Spofford  et  al.(l9S2)  is 


To 


( 1 57; 


We  are  dealing  with  the  amplitude  or  pressure  ratio  so  let 

r/v,  ■ 


‘B6KJ)  - 


Rt“"5f,l£h 


(tea) 
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We  need  the  derivative  dR  .  Taking  the  range  formula  from  Subroutine  GCOMP  or 

dBA 

equations  (16)  and  (21)  ofjjgfjWgW,  with  appropriate  substitutions  of  angles 

*3  (ISV 


gives  R  =  D1  +  Fj_  F-  tan  Q. 

tan#/  90 


where  F^  =  2BC3 


CK  \A  (?( 


-[ 


I  ± 
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(I  f  8) 

d  T2.  X>  I  F/  r  cJ  ^"2.  •/  o  1  7—*  1 

*-Fzsec 


(  I  +  si*&fc0s&^)7  .(7*0 


(??/) 


fj;  ^2.  _  <p(  ^2  c?l  ^3  , 

do,  ~  A  @3  AS/ 

From  Snell 's  law  C  o<,  &,  -  J^.L  cosS? 

c5  S 

-  s  i  w»  9,  -  -  O.-  s  ;y\  @2  d  0% 

C3  cT&t 

d  B>  3  __  ^  s  ;  vn  0  > 

c,  s  f  n  S3  ) 


C^/  Of 


a(F- 


7V9,  =  2Bs;«e3»s^[  ' 


(V?*J 

differentiating  gives^fjl 

(t<iu) 

(  f9S' ) 

;  fc_ 3)  s)*S*,  (,?c 


s: «  &3 


c\K  _  _  DJ _ _  .  hj  \  L.2±<n^, - __ 

Te,’  s\T9t  (2SC) 


Fj  ^CzsinOt 


I  x.  /Q  —J 

3  I  5 r\  &3  c~e>% 


(  n 


+  ®3 jj  <■ 


+•  F: 


^6>, 


r,  (  Cl.  S\n@!  _N 

(T/S/h^c^^ 


Dl  4- 

slT&i  *3o  . 

"l  (  4- 

"  2,3  "5^ ^ &3  2 B 


^  /??^) 


^3 


2B  2B.s  :  ^  &3,c<?£@£ 


c(  Ic 
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£j2By[0,s,'n&lcc>s&s 
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<1 
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J>l 


Pi 


B3o 

(  n-b)  +- 


Ci.  s'  *  B  i  t 

c, 


(**0 


Bz 


$  '  C©5  ©, 


vf" 


CoS'  @1 


? 


^9,“  +  B3o  V  CI  sTT^Z^-^j/L  j  ? 

oll^  'Ll  ,  [  i-f-B*)  /FtVCs  si*  @1  —  Vj  4^  ,^>?~^^3-)  ^2c  J 


siw\ 


2.  ^  B  l  *3  *>A  ^7  5  »  *  os2!^, 


£L1  ---M—  4. 

^9,  ^  0, 


^9/ 

5cUsf;tuti^  /Ue  express?;?*,  for-  F,  3'^-S 

r  >  y-B  V26CsV^S*^^i.Y  I  4.  Sj±**±)  (  *c 
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£43  f  A  »7  B43 

Z)7 


-)  ( ^7 


ok  £  •  K 


■f  K  e  SU^OUT >bJ£  FF/^F. 


38 


This  is  the  derivative  of  the  range  with  respect  to  the  bottom  grazing  angle 
evaluated  at  the  Big  Angle  in  the  water.  With  this  derivative  DRDT  we  com¬ 
pute  BGI(J)  as  in  the  subroutine  for  the  Big  Angle. 

ARC  LENGTH 


Our  model  assumes  that  loss  mechanisms  in  sediments  behave  linearly,  e.g.  loss 
L  =  K  •  S  •  f  by  equation  (13)  of  Spofford  et  al.(1982),  or  in  terms 
of  symbols  used  in  our  program  DL  =  aoSF^,  so  in  the  program  it  is  assumed 
the  frequency  averaged  or  effective  value  of  the  sound  path,  ARC  LENGTH,  in 
the  sediment  is  given  by 


S  =  DL  where  DL  is  dissipation  loss  in  the  sediment  given  by 
aoF., 


equation  (23)  of  Spofford  et  a!.(1982),  or  as  shown  in  Subroutine  GEOPHYS,  ac 

is  the  average  value  of  attenuation  constant  for  this  path  length,  and  FK 

denotes  the  center  frequency  of  a  one-octave  measured  bandwidth.  Thus,  let 

BLS  =  DL  and  in  terms  of  the  equation  for  DL  with  lower  limit  01=O,  upper  (2&9.  ) 
cxoFk 

limit  00=63  and  the  factors  from  this  subroutine,  D8  through  D14,  we  get 

R.  R  2  C_2x _ _  f  Xd  (si*  (BX_ C3 

"  (  lf0)ctfS©,  6<r  [cas&ij  4^  cos  &3 

<s>3  ,  ( I ±  29 CpO-  ±  g ;  n 

2c-t<,s30s  2(Tcos&3 

as  in  the  subroutine  statement  for  the  ARC  LENGTH 
TRAVEL  TIME 

We  need  the  total  travel  time  of  a  path  with  bottom  grazing  angle  BA  (Big 

Angle).  We  first  find  the  travel  time  in  the  sediment.  By  definition  C(Z)  = 

dS  where  S  is  an  element  of  path  length.  The  vertical  component  of  dS  is 
dt 


dz  =  sin©  dS  so  dS  =  dz/sin0,  and,  noting  that  0  =  G(z),  we  have 

dt  =  dS  =  dz  from  which 

c (z )  C(z)sinO(z) 

dz 


'’V 

Z 


Ts  =  2 


c  Cz^  s  /  ^  (Be  z.) 


(xn) 

( 


This  is  equation  (19)  of  Spofford  et  al  . ( 19S2)  where  z  represents  the  turning 
point  depth  or  basement  depth  for  rays  which  reflect  from  the  basement. 
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Snell's  law  gives  C(z)  =  c0  cosOU^  where  c  is  the  sound  speed  at  the  top  /2/3J 

COS0  v 


of  the  sediment  and  Q0  is  the  angle  at  the  top  of  the  sediment.  Differen¬ 
tiate  this  with  respect  to  z. 

J9 


—  S  }  G  ~ —  C  c  &  ol  B  C2.  ~)  .  Cej  Bo  Q  fj?  ) 

^  ~~c7~  dz.  ~~cl  y  ' 

-J-—  i  -  Co  dO 

sfv^  ^U)cos9o 

Ts  =  7  p  c»  d* 


'o 

- r  __  gc» 

^  C-osOo 


5  Cz)  C (z.)  cos  Be 

e‘  de 


i$  g Cz.)cu) 

r  -  2Cc  Cs'  d£ 

s  C~0°X,  kL!fZ2g>c<2) 

C  Cz)  +  -BG> 

T*-  -  —  -2__  ( C  (z)  +  ZCe>  )X  & 

e“* J  3o(  H-3)C(*> 


i<  ~  - 


5  0  (  )  f-B )  cos  Bo 


0  *i W9 


(2.H) 
Qzts  ) 
(ztc) 

(znj 

tub) 

(2/9) 

(2  icj) 


<9 


lJT 


Tc 


(  1  +  %) 

Z 


5  "  %Ot3)  [ 


y  Cos  Bo  c  OS  & 

r  x=>  e°  -e> 

I  — 2— -  +.3  X  (sec  e>  -hta^  e) 


de  c^O 

(2  z  2V 


Be 


Ut:*q  c\  tr[a  identity  tUe  /imt?  3  es 

J  J  /  r  ^  -  &o 

__  'P*  I  &  'D  I  -  (  ^lB~L  ®  ^  j 


g,<r /+ B)j  ces&‘> 


*-!>  In  A 


Ce>s  & 

I  -s; ^  S 


\&~Bc 


(  22  J 
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This  is  equation  (22)  of  Spofford  et  al.(1982).  The  time  for  travel  in 
the  water  is  the  slant  range  divided  by  the  sound  speed  in  the  water,  i.e. 
Dl/C/sinBA.  Adding  this  tp  Ts  where  we  substi tute  ^/=0  for  a  path  refracted 
in  the  sediment,  'feAS  ,  a.Y\&(  s!  ^  yieMs 

^  _ fj 

CiBSI  <3*  (H 


7~Z 


B/?3 


<35 


B03 


ffs)]  ^ 


fj  0<AJ 


Ci >s  & 


&  \'/z. 


S  O 


T  -  - 

'  C,B5I 

J2± 


.  (  Ltsj*  &_\ 

-  \ 

-Z3 — r  i 

Secs  B#3> 


&43 


-ir^[T= 


%S>  3 
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T  = 


er,  SSI 


V-. 


BTCr)r  -grp  +  ■’<-  Hi)  *■  J-(  IfcH] 


which  is  the  expression  for  the  travel  time  in'  Subroutine  REFR. 


DELTA  T  SOURCE  RECEIVER  (for  Big  Angle) 

8SS  =^l  -  g?-  Self  =  ]j  I  - 

BC1  =  cos  BA  yields  CS  cosBA  =  cosOs  so 

U 


(227j 

Snel 1 1 s  1  aw  wi th 

(ZZ  S') 


BSS  =  sin  0s  where  0s  is  ray  angle  at  source  (or  receiver)  (27<i') 

BDTR(J)  =  D6(ZR)sin0s  =  2^ZR\  sin0s  is  the  travel  time  difference  for  the  (^Z3o) 

Big  Angle  path  from  source  to  bottom  to  surface  to  receiver. 

Similarly  BDTS(J)  =  2  /ZS^  sin0s  is  the  travel  time  difference  for  the  Big  (_  2,3l~) 

cs 


Angle  path  from  source  to  surface  to  bottom  to  receiver. 


We  want  to  deal  with  rays  not  hitting  the  basement  -  refracted  in  sediment  only, 

so  the  statement  IF  SA>THIB(=cos"^  Cj_)  go  to  999  stops  computation.  Sound  ^232j 

C8 

speed  in  sediment  at  basement  level  is  CB  and  THIB  is  the  bottom  grazing 
angle  of  the  refracted  ray  which  reaches  the  bottom  of  sediment  at  basement 
level . 

AUX  =  FACS*SC1  (  733 ) 

=  “cf  cosSA=cos0s  where  0s  is  angle  at  source  (or  receiver)  for  a  ray  (  *)  3  <-l ) 
with  bottom  grazing  angle  SA.  \L  ' ' 
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IF  RAY  DOES  NOT  REACH  SURFACE  NOT  COMPUTED 

IF  AUX>1  go  to  999  stops  computation  if  cosft^l  which  is  not  physically 
realistic. 


DELTA  T  SOURCE  RECEIVER  (for  Small  Angle) 

SSS  =  y  l-(AUxj^'  =  sin0s  for  bottom  grazing  angle  SA. 


(23.5) 


SDTR(J)  =  2  (1R)  sin^s 


SDTS(J)  =  2fZSjsin0s 


IS) 

is: 


Travel  time  difference  for  indirect  path  forfeit 
surface  reflections  at  receiver  and  source  .(2 }  7/ 
Use  of  Snell's  law  and  trig  identity  is  made  to  define  the  following: 

O 


SSI 


■f 


(SCI)  =  sin  SA 


SC2  =  C i  cos^A  -  CcS  where  ^is  the  angle  in  the  layer  for  ray  with  /i  ^ 

Ci 

bottom  grazing  angle  SA, 


SS2  -  Ul-(SC2f  =  sin  &  , 


(2V»J 


SC3  =  C^(SClj  =  £3  cqsSA  =  coj£^  where^  is  angle  at  top  of  sediment  for  bottom  few> 
C-i  Ci 

grazing  angle  SA, 


-  (SCzf-  =  I//- 


SS3  =  |J  / 

SMALL  ANGLE  REFLECTION  COEFFICIENT 
A  ,  r 

A 2_-  '  *  &L 

A3; 

SV3Z(J)  =  (A3-A2)/(A3+AZ)">) 

SV2](J)  =  ( A2.-A  / )  /  ( A2+A  J ) 


Impedances  corresponding  to  Small 
Angle  ray. 


(z^ ) 


) 


STRF(J)  =  16A| At  /  A  a. 

UA3+Ai)(A2+A,  ) 


kz 


SSIN2  =  sin <9^  for  small  angle  ray. 
GEOMETRIC  AMPLITUDE  (for  Small  Angle) 


Reflection  coefficients  for  boundaries 
designated  by  subscripts.  (  2  </</  ) 

Two-way  travel  layer  transmission 
coefficient  factor.  ^‘/J  ) 

Ml  ) 


DRDT  -  same  as  DRDT  for  Big  Angle  case  with  BA.BA3  changed  to  SA,SA3  etc. 
This  is  the  derivative  of  the  range  with  respect  to  the  bottom  grazing  angle 
evaluated  at  SA. 


42 


sgi  (  j)  =  \r 


cos  S/ 9 


Rsiw,  S'/HdR 

TO 


p0 


the  pressure  ratio  ^2*/ ? 
for  the  Small  Angle  path. 


This  is  the  geometric  spreading  loss  factor. 

ARC  LENGTH 

SLS(J)  is  the  same  as  BLS(J)  with  angles  changed  to  those  appropriate  to  the 
Small  Angle  path. 

TRAVEL  TIME 

ST ( J )  is  the  same  as  BT(J)  with  angles  changed  to  those  appropriate  to  the 
Smal 1  Angl e  case . 

999  TH=TH-DEL  increments  the  angle.  The  computations  are  made  in  steps  of 
1  degree  for  all  applicable  angles. 

RETURN 

END  (of  Subroutine  REFR)  (Go  back  to  the  main  program.)  CALL  BASE 
SUBROUTINE  BASE 

The  statements  NTIR ,  DEL,  TH,  N=90-MAX9(N,NTIR) ,  and  DO  10  1=1,90  cause  calcu¬ 
lations  to  be  made  in  steps  of  1  degree  for  angles  from  the  angle  of  total 
internal  reflection  or  the  apparent  bottom  angle  for  a  ray  vertexing  at  the 
bottom  of  the  sediment,  whichever  is  larger,  to  90°.  This  insures  reflec¬ 
tions  at  the  basement  are  accounted  for.  The  Z's  and  FAC's  are  the  same  as 


used  before . _ 

CB  =V2  3»(H-bXc3)8z  f  Cf(n-B)2  "  BCj 

gives  the  sound  speed  in  the  sediment  at  the  basement  level 
D1  through  D14  are  factors  used  later. 


(Z¥9) 


BSA  ^'rr/z  ~l  X/O  &  (used  in  Subroutine  BASANG)  (  2 

TH3  =cos  ~&3  is  angle  at  top  of  sediment  corresponding  to  bottom  ^  < 

grazing  angle  for  a  ray  which  reaches  the  bottom  of  the  sediment. 

R  is  the  range  for  a  ray  refracted  in  the  sediment  with  vertex  at  the  basement. 
This  expression  is  found  by  substi  tuti ng0e=C^,  &'=0  i nto  the  equation  for  Rs 
g  i  v  i  ng 
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Kc  = 


<-3 

60  (  )  +  5)  1  os  8$ 

|^Ctfs6>3 

C-i 

r  9j 

<$o  (■  14  ^3  , 

|  Cos  £3 

I 

CD 

< _ 1 

So(  &},  1 

Ct  I 

(93  -k 

g^>C  l+B)c.e>S  <5^1 

t^r-asf 

:?3  + 

f? 


+  B  s  *2  &3J 

i  *2 &3  -f-B)(z)] 

IT  D/*~  L.  - 

-  »i  — f-§Jt(c.2  <r  +3)s;*ZG’3 

U~  ©m  Cas2<©s  [  2  v 


) 

(Zi’2.) 


AUX  =  tan  "*  Dl_  This  is  the  apparent  angle  for  a  ray  tangent  at  the  bottom  (z£2>) 


Program  modified  by  addition  of  these 
statements  to  insure  completion  of  a  run. 

clears  registers 

orders  computations  for  angles  in  the 
range  mentioned. 

TH=TH-DEL  increments  angles 


DO  10  1=1,90 
BSS/N2( I )=0.0 
BSV 32 ( I )  =0.0 
BSV21 ( I )  =0.0 
BSGI(I)  =0.0  J 

10  BSTRF( I )=0 
DO  100  I  =  1 ,  N 
J=90-I 
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CALL  BASANG 


SUBROUTINE  BASANG 

IF  TH>TIR  GO  TO  100  allows  computations  to  proceed  for  angles  greater  than  TIR. 
Otherwise  print  "SPECULAR  ANGLE  OF  BASEMENT  RAY  LESS  THAN  TIR"  and  computation 
stops.  CB  is  same  as  in  Subroutine  BASE — sound  speed  at  basement  level  of 
sediment.  The  next  several  factors  are  like  ones  used  before. 

BEST  =  BSA  is  estimate  of  basement  reflected  ray  angle  in  water. 
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AUX  =  cosBSA 

THB  =  cos  -jangle  ra-X  at  basement. 

TH3=cos  ^  cc>s  ^  angle  of  ray  at  top  of  sediment. 

Now  find  the  range  in  terms  of©0and£^as  limits  from  previous  equations  for 
R  and  Rs. 


) 

USS) 
(2 SC) 


Ql 


4.  Sh 1  ^3  cosgj  +£  & 
Cos  Sh,  c  °  £ 


%l  >**)“*$!  t. 

_  Sb  _  i9g-  _  2.  B  st  *  ®bJ 

c  OS  £h>  c  Gs  ^-3  * 

©D  —  ^^3  OUa<?I 

I>2 —  r  ^  Cox-(^)?~.s-.’»i<^gggs^s 


(2X7J 


/?  r  — — 

3  Ce?s2Q3l  Ce>s 

+  g  ,  .^2-  Cos  ZSA  3  -  st*  &jsJJ 


Rs* 


■h 


(— c 


:  1  5  ;  *  2  £>3  -  S  I  ^ 

osEssfl  2  4 


f  2-i?) 


£  OS 


55/?  /s  /  h  Bb  -  s  5  *  & 


B 


/  U-'s  /'£ 


the  second  term  in  VAL.  From  the  formula  for  R  we  see 
P) 


VAL  = 


t CL  va  s  / 


^2  r  02-&G  ,  s'1^2. Qz-si^^ @b  / 
-f7Z^f[~2~  4  r 

^  ,  L  ^  \1  T>1 


ZG  CJ 


+  jj£>  COS  BS  A  (s’ In  Q3  -  s- 1  *  £2 


s 


i  T~H 

where  TH  is  the  apparent  bottom  angle.  This  VAL  is  similar  to  the  one  used 
in  SUBROUTINE  REFRANG  except  basement  reflection  occurs.  See  Figure  15. 

Now  the  IF  and  DO  statements  start  iterations  which  adjust  the  angles  to  make 
VAL  approach  zero  until  convergence  is  obtained  when  JdElJ  4  lXlO'^whereupon 
BEST=BSA.  If  convergence  does  not  occur,  print  "BASEMENT  ANGLE  DOES  NOT  CONVERGE." 
and  stop  computation. 
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The  basement  reflected  ray  angle  at  the  top  of  the  sediment  is 

BSA3=TH3=  ©3  and  BSA  is  the  basement  reflected  ray  angle  in  the  water.  We  (2C/J 
use  BSA  and  BSA3  in  SUBROUTINE  BASE. 

RETURN 

END  (of  subroutine  BASANG)  Go  back  to  SUBROUTINE  BASE. 


IF  (BSA3.EQ.0.0)  GO  TO  100-added  to  insure  completion  of  run. 


BC1  =  cos  BSA  (  z  L  Z  ) 

BS1  =  sin  BSA  C Z- ?  ?  ) 

BC2  =  C2  cos  BSA  =  cosOo  ,  i s  in  the  layer.  (  2  C  <-/  \ 

C|  J 

BS2  =  sin  ©£  in  layer  £  2.LS  J 

BC3  =  C3_  cos  BSA  =  cos  @3  »  ©J  is  at  top  of  sediment.  (llc.') 

Cl  ' 


BS3  =  sin  ©3 


C  2.07) 


BCB  =  £B  cos  BSA  =  cos  &Q,  ^jis  basement  grazing  angle  in  the  sediment.  ^ 
C  f 


BSB  =  sin©0 

DELTA  T  SOURCE  RECEIVER 


(ZC9) 


AUX  =  Cs_  cos  BSA  =  cos&s  where  ©5  is  angle  at  the  source  for  the  basement  ) 

C| 


reflected  ray. 

IF  RAY  DOES  NOT  REACH  SURFACE  NOT  COMPUTED. 


IF  AUX>1  go  to  999  stops  computation  if  cos©j>  1  which  is  not  realistic 
in  nature. 


SS  =\j  I  =  fl  ~  Co ~?Gs  * 


(Z  7() 


S  !  va  ©  * 


BSDTR(J)  =  2lcs 
BSDTS(J)  =  (  *  &<=> 

LAYER  REFLECTION  COEFFICIENTS 

A  |  r  Zi/f  1  ^  $ 

A  2  -  Z-2/ 3 '  *  &2. 

A  3  -  Z,/sU  ©3 


Travel  time  differences  for  receiver 

and  source  respectively.  (zm) 

Acoustic  impedances,  medium  designated 
by  subscripts.  (173) 
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BSV32(J)  =  (A3-A2)/(Ai+Az)  1 
BSV21(J)  =  (ArAi)/(A2+A/)  J  . 
BSTRF(J)  =  16  A,A3  /  A 2. 


Reflection  coefficients  for 


^  2^  ) 


boundaries  designated  by  subscripts. 


This  is  a  two-way  layer  transmission  coefficient  factor  for  basement  reflected 


ray. 
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C  2  7c) 


BSSIN2(J)  =  sinO^,  6g  is  angle  in  layer. 

GEOMETRIC  AMPLITUDE 

BSAB  =  cos"*(cos9g)  =  Og  basement  reflection  angle  in  sediment  at  basement. 
We  now  find  the  derivative  of  the  range  with  respect  to  the  bottom- grazing 
angle,  9^ . 


r  =  +•  _  f  S->'o&jcosB3  -  c  c>s£faj  (Z-7  7) 
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This  is  the  expression  in  Subroutine  Base  for  DRDT  where  the  quantity  inside 
•braces  is  DRDT3,  6j  =  BSA  and  other  angles  and  their  functions  are  as  mentioned 
above  in  discussing  this  subroutine. 

R  =  Dl/tanTH  gives  the  range  in  terms  of  TH. 

C  25-sO 

ARC  LENGTH 


BSGI(J)  = 


Cos 


is  the  geometric  spreading  loss  factor. 


is  found  in  a  manner  similar  to  that  in  Subroutine  REFR  but  with  use  of  Qg 
for  the  lower  limit  in  equation  (23)  of  Spofford  et  al . (1982)  corresponding 
to  the  basement  reflection  angle.  This  gives 


BSLS  = 


D  L 


2  C3 


0C0  ccs  & 5 


zZSt 


IkXXi.  tX  TX.  £2.  f  Sr  >1  c  os  —S’!  *  &g  os  &e) 

4c~  CosXOj  \  ' 

f  ref  +  ; * e3 -s ; -  £s) 

(2 crcos^Oj  2.0“  cos&g  J '  ' 

,  ( iXL ^3  ±2jycz_  ^  b y ^ gcC^jf 
+‘  v  4-cr  cosz  2.  <y~  yJ 

TRAVEL  TIME  2  8  W 

is  found  as  in  Subroutine  REFR  with  Og  substi tuted  as  the  lower  limit  in  the 
equation  for  TS,  equation  (22)  of  Spofford  et  a1.(1982) 


(28C) 


T  = 


2>/  _ _ 

C  <  Si*&!  t3o(l+S)l  C0&&3, 


-f 


t 


*&3Csi*e&  n 


5) 


l Si  H<^3j 


and  by  the  identity  mentioned  in  Subroutine  REFR 


Li 


BST(J)  = 


C  Os  ©3 


8  /  / 1  ±s :  VI  £3/ 1  ~S l  VI  (9g^ 

2.  ^  f-S  i  *&Bjj 


0 


The  computations  in  this  subroutine  are  made  in  1  degree  steps  for  all 
appl  icable  angl  es . 

RETURN 

END  (of  SUBROUTINE  BASE)  (Go  back  to  main  program.) 


i 
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SUBROUTINE  LOSS  (JJ) 

Discussion.  This  supplements  section  III.C.  Simulated  Bottom  Loss,  as  set 
forth  in  Spofford  et  al . (1982) .  , 

Other  subroutines  provide  for  computations  of  travel  time  for  base  path, 
angles,  reflection  and  transmission  coefficient  factors,  surface  scattering 
travel  time  differences,  etc.  which  are  needed  in  this  subroutine. 

For  a  base  path  travel  time  of  T^  j  the  three  other  paths  have  travel  times 

+  (2  zs/cs) 

Tk,  3  =  *TT,  I  F  (2  2.X./CS)  s-,*9s 

Tk  ^  -Tk  ,  +  (2J2LS/cs )  3  <  (jL'Z.'R/cs) S>'  *\&s  (Z'rc  ) 

The  travel  time  differences  are  converted  to  phase  shifts,  and  the  paths 
acquire  additional  phase  shifts  for  surface  reflections  of  -IT',  ITj  3-7T 

radians  respectively.  The  subroutine  accounts  for  all  these  effects  by 
treating  the  surface  scattering  travel  time  interference  "artifact"  as  a 
separate  factor. 


It  is  assumed  the  frequency  averaged  intensity  ratio  is  given  by 

C  F2. 

P*P  C IF 


7-o, 


Fa-F, 


I ) 


r, 


where 


P  represents  the  sound  pressure  and  F  is  the  frequency. 


Ordinarily  one  assumes  a  superposition  of  waves  results  in  pressure 

j<t>K 


so  one  would  write 


K  '>  .  L 

but  in  this  model  it  is  assumed  'P=  ^  ^  ***.».] 

o* 

where  some  justification  is  by  Spofford  et  al.(1932). 


--  (f  a  eitpLj, • 

(Vr'i  J 
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If  we  include  the  frequency  dependence  in  yk  j  so  this  represents  phase  per 
unit  frequency  we  get 


J 


and 


(4-)  =  <«**) 

V  e/  a  4  J 

Now  let  ^A*  A,^  be  the  frequency  averaged  products  of  amplitudes;  they  can  be 
factored  out  of  the  integral,  and  note  F^-Fj  =  F^  for  1  octave  bandwidth,  so 

^  w  a  \  I  •/  A.  J _ U  \  — 


W  ■  *■  $.<«**  C  «' 


■4£)p; 


1  <  ‘f’K!  - 


(vr?) 


I  IVKj 

For  paths  in  the  sediment  the  phases  are  assumed  to  be  complex  and  written 
(tyk*Xir(t+iA.»s) wi th  t  the  time  in  the  sediment  (plus  surface  scattering  travel  (zj$ 
time  difference),  «{©the  absorption  coefficient,  and  S  the  effective  path 
length  in  the  sediment.  The  exponential  terms  account  for  absorption  in  the 
sediment.  But  for  layer-reflected  paths  it  is  assumed  there  is  no  absoption 
and  the  imaginary  part  is  zero.  Thus  the  phases  are  real  and  contributions 
from  these  paths  are 


(  2^) 


fy 


<L 


j 


Consi der  the  first  sum.  For  layer-reflected  rays  of  unit  incident  amplitude, 
the  reflected  amplitude  is  simply  taken  to  be  the  reflection  coefficient,  V, 
and  intensity  is  V*V.  The  frequency  averaged  value  of  V*V  is  desired.  It  is 
obtained  by  computing  the  values  of  V*V  at  eight  different  frequencies  over 
a  one  octave  range  and  finding  the  numerical  average.  The  number  of  values 
used  in  the  average  evidently  is  arbitrarily  chosen  to  be  eight. 
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Consider  the  second  sum.  Four  bottom  interacting  reflected  pulses  reach  the 
receiver  from  the  source  in  some  order  arbitrarily  chosen,  for  example,  on  a 
time  scale  as  illustrated. 


Pt 


VR  — >|^- 

"uV - 2V 


PS 


t: 


T>S- 


— * V 


2.T rt 


The  sum  includes  all  combinations  of  travel  time  differences  such  as: 


K,L 

Phase  Shift/P 

Surface  reflection  at 

1-2 

DR 

receiver 

1-3 

DS 

source 

1-4' 

DS  +  DR 

source  and  receiver 

2-3 

DS  -  DR 

source  and  receiver 

2-4 

DS 

source 

3-4 

DR 

recei ver 

The  frequency  averaged  sum  is 

sj wi  FjJvvi  +15) -si*  Fi_(VR+ds) 

&-R.  +  DS 


suifz,  {P£iPgl^siri  El  [21d>El 

VS  ~VR 


2.F/  L  &K+-DS  VS-VR 

o  _  q  2 i  h  E^.  &_R  zJL!j±5l HR. 

c  r >s  DR  J 


for  the  layer  reflected  paths.  The  terms  corresponding  to  only  one  surface 
reflection  are  given  negative  signs  to  account  for  a  phase  shift  of  T f  radians 
while  those  corresponding  to  two  reflections  have  positive  signs  to  account 
for  a  phase  shift  of  2ir  radians.  The  terms  corresponding  to  no  phase  shift 
are  accounted  for  by  adding  the  numeral  4  to  this  sum.  The  sum  just  mentioned 
is  multiplied  by  the  frequency  averaged  value  of  V*V.  This  accounts  for  the 
laye^  reflected  ray  intensities. 


Next  the  layer  transmitted  rays  are  accounted  for  by  computing  the  layer  two- 
way  transmission  coefficient  at  the  center  frequency  of  the  octave  wide 
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band  and  multiplying  by  the  basement  reflection  coefficient  (where  appropriate) 


and  the  geometrical  spreading  loss  factor.  This  gives  the  amplitude  of  the 
transmitted  rays.  No  further  frequency  averaging  is  done  for  transmitted  path 
ampl itudes . 

The  phases  of  the  four  layer-transmitted  signals  are  accounted  for  by  changing 
the  signs  of  two  of  the  amplitudes  which  is  equivalent  to  changing  the  phases 
by  Tf  radians  to  account  for  one  surface  reflection  each  and  no  change  of  the 
other  two  to  account  for  phase  change  of  zero  for  the  bottom-reflected-only 
path  and  phase  change  of  2 Tf  for  the  twice-surface-reflected  path.  In 
addition,  phase  shifts  and  attenuation  in  the  sediment  are  accounted  for  by 
forming  complex  frequency  dependent  phases  as  follows: 


C3*i) 


where  t  is  base  ray  travel  time,  tp  and  t^  are  surface  scattering  travel  time 
di fferences ,  <£0  i s  the  sediment  attentuation  coefficient,  and  S  is  the  sedi¬ 
ment  path  length.  These  phase  factors  account  for  attenuation  in  the  sediment. 
Products  of  amplitudes  are  formed  as  are  the  exponential  phase  factors  and 
their  integrals  are  summed  as  needed  to  produce  the  transmitted  ray  portion  of 


The  reflected  and  transmitted  frequency  averaged  intensities  are  combined  in¬ 


coherently  to  find  the  total 
equation  (43)  of  Spofford  et  al.(1982). 
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The  computed  spreading  loss  is  given  most  simply  by 
L°~  1^)0  (sja^fra^e) 


*  Z  --  -JOh,0  -4  —  •  (*<* ) 


(  k  y 

c<ps&) 


which  represents  a  sum  of  4  unit  initial  intensity  rays. 


Therefore,  bottom  loss  is  found  by  subtraction 


C  3  0  1) 

(-50?) 


Suitable  variations  in  computational  procedures  are  made  in  order  to  suit 
the  various  domains  to  be  accounted  for. 


Now  consider  SUBROUTINE  LOSS ( JJ) . 

ALPH0=  1 nlO  ALPH9  is  used  to  convert  the  usual  logarithmic  expression 
4  if  U  oooc?) 

in  dB/m/kHz  from  the  COMMON  block  to  suitable  units  for  use  in  an  exponential 
form. 

EX=log,^e  is  defined. 

Aim=(o+i)  used  to  provide  imaginary  quantities. 

F.  =F/ |/“F is  the  lower  limit  and  (  30<i  ) 

F^=F  if  2  ' is  the  upper  limit  of  a  one  octave  bandwidth  centered  at  F.  C  "3/0^ 

FPI2D  =  47TFD  Thickness  of  the  layer  is  D  and  sound  speed  in  the  layer  is  C ^  (3//  ' 

C2. 

NAVI =8  is  defined. 

FPI  INC=FPI2D(  l JT  -1 Al)  =  4tTD  (F^-F,  )  is  an  increment  l7-) 

NAVI  C2  8 

FPI  IN  =  41TD  F  -  (8-1)  4fl*D  (Fa-F,  )  L  3  /  3  J 

C2  "2  8C2' 

=  4tT0  (F-  7  (F- -F.  ))  ( 

16 


57 


This  is  an  increment  to  be  used  in  computing  frequency  dependent  phases  for 

reflection  and  transmission  coefficients.  Consider  a  one  octave  bandwidth 

centered  at  frequency  F(JJ)=F  with  lower  and  upper  limits  Fj  and  F^.  The 

term  FPIIN  is  used  to  compute  a  quantity  at  the  frequency  Fc  =  F  -  ^(F^-F^) 

16 


as  illustrated  in  figure  17. 


p,  '<j  /  7.  *?or  C  t<  »  K  ®  ~F  i  'O 

Ok_  pne  -  c  r  el  e  b  a.  . 

This  is  the  starting  frequency  for  computing  frequency  averaged  reflection 
and  transmission  coefficients. 

NTHC  =  ©cthe  caustic  angle  in  degrees 

n  =  90  -  ec 

.COMPUTE  LOSS  FOR  ANGLES  ABOVE  CRITICAL 

(the  caustic  angle).  These  calculations  apply  to  the  domain  of  range  R^as 
shown  in  figure  2. 

DO  10  I  =  1,90 

10  ALOSS( JJ,I )=25  This  puts  in  an  assumed  value. 

DO  100,1=1 ,N 

J  =  91-1  Computations  are  made  for  angles  in  steps  of  1  degree  from  <9C  to 
90°. 

FP I AUX  =  FP 1 1 N  Use  in  starting  the  frequency  average  calculation. 

DO  20  NAV=1,NAV1  Calculation  is  made  for  8  frequencies  within  the  one  octave 


bandwidth. 

F  -£(Fi-p.)]s:«ea  (3lfc) 

is  a  frequency  and  path  dependent  phase  angle  for  the  layer. 


4-tTD 

ALPH2=  cC^  -  -£T 
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C=cos«>C2.*''|  Used  in  computing  the  square  of  the  layer  reflection 

S=sin*^2.  J  coefficients. 

The  square  of  the  layer  reflection  coefficient  is  needed  bepause  we  want  the 
intensity  of  the  reflected  wave  which  is  proportional  to  the  layer  reflection 
coefficient,  or  equal  to  it  for  an  incident  wave  of  unit  intensity. 


The  layer  reflection  coefficient,  equation  (30)  of  Spofford  et  al . (1982)  is 
2.  j  <£.  d 

R  =  V|2.  +  V23  e _  and  if  we  let  oCg,  -  2e*-4  f  3ll) 


R*R  =  ( v.zW'*2  Vv.,  4  Vg^e1^  \ 

-  ^  V3Z  V)2.  f  g  2) 

I  *■  ^>2^2  3  (e e  ^z) 

_  \/j2=.  ^3  2^)  Z_C_£^S£^Z  jj  V2.3 

+•  2  VV2A/23  cos^2. 


2l  V 

 (  +-  Va3  C  4  ^ I/2J  ^  /  noC2/ 


0  +  ^ ^2.) “V{V, 2 V2 3, 


.5  > 


(  3/8  ) 


Interchanging  subscripts  changes  signs  of  the  factors  but  has  no  effect  on 

the  value  of  R*R  so  doing  this  and  introducing  prefix  R  we  find 

VV=R*R=  +  %l/32.  4-  (  W32  S/n  eC2,?~  /  3/^ 

as  in  Subroutine  LOSS.  The  boundary  reflection  coefficients  RV21,  RV32 
are  computed  in  Subroutine  REFL. 

AUX=AUX+VV  The  value  of  VV  is  computed  and  the  results  summed. 
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FPIAUX  =  FPIAUX  +  FPIINC  causes  increment  in  frequency  to  the  next  step  shown 
in  figure  17  and  computation  of  VV  again. 

20  CONTINUE.  Procedure  is  completed  8  times  yielding  a  su^,  of  8  values  of  VV. 

VV  =  AUX/NAVl  =  VV(Fc)+ — +VV(Fc8)  gives  an  average  value  of  VV,  the  squared 

8 


reflection  coefficient,  over  the  1  octave  frequency  band  as  illustrated  in 
figure  17. 

DR  =  2ftRDTR(J)  ~)  converts  surface  scattering  travel  time  differences  so 

/  (32-cO 

DS  =  2ttRDTS(J) J  when  multiplied  by  frequency  gives  phase  shift  in  radians. 

DT  The  first  four  consecutive  DT  statements  provide  the  sum  of  frequency 
averaged  phases  for  the  layer  reflected  waves  as  described  in  the  Discussion . 

The  fifth  DT  statement  provides  the  sum  multiplied  by  the  frequency  averaged 
reflection  coefficient  to  obtain  the  intensity  of  the  layer-reflected  rays. 

For  basement  reflected  paths  the  separate  ray  contributions  are  summed  inco¬ 
herently  on  the  assumption  that  the  time  delays  between  layer  and  basement 
reflected  paths  results  in  incoherent  arrivals  at  the  receiver. 

ALPH2=FPI2D*BSSIN2( J)/2-47TF( JJ)  sinQ^  =oCgives  phase  shift  for  1  way  travel 

c2.  2  2  (  32  l) 

in  the  layer. 


C  =  cos^ 

S  =  sin  «<-£ 

The  two-way  travel  layer-transmission  coefficient  can  be  obtained  as  shown 


for  SUBROUTINE  GEOPHYS: 

T  =  /  8  ft  i  [  $  2. _ \ 

_ ( _ 

1  +-2K52.V2.J  TosZ^2~i- 

where  we  have  in  the  numerator  the  expression  for  BSTRF. 


( 3  it) 


Consider  the 


denominator  of  this  expression. 


It  has  the  form  1  +  2M  cos2c<2_  +  M^,  where  (^2 3/ 


M  =  V32V21,  so  by  using  identities  we  can  write: 

S  \  +.  Cc'i'eL^  +  i  M  (oosV’i  -  s  <  i  c  ^  j 

-  (  /+  ZAI-f-AI-2")  1-2- i  £  ) 

=  (  4-  (M  -l)ZS.'  ^<£a.  (HG) 

-  [{  l+rt)  Cc  4-  1  (M -l)  (yZ-l) 


Substitution  of  this  expression  for  the  denominator  gives 

V  =  BSTRF _ _ 

£  C  i  -f-S-SV^S  BsvFf)  771*^  +  ]  ( WsvJzBsV2]  -  TJZ 


CHS) 


The  quantities  with  prefix  BS  come  from  SUBROUTINE  BASE. 

The  amplitude  of  a  ray  reflected  by  the  basement  after  passing  through  the 
layer  into  the  water  is  given  by  VV  = 
where  BSGl ( J) 


is  the  geomatric  spreading  factor  and  (REF) 


i  *1  ©j 

is  a  numerical  basement  pressure  reflection  coefficient. 


HI'-  wu  I,  u  i  Up  i  (-'u  j  j  i  t  w  I  I  f  w  y  wtlt 

C3«> 


(3 


AUX  =  0  is  defined. 

IF  VV=0  GO  TO  120  accounts  for  the  case  of  no  basement  reflection. 


The  amplitudes  of  the  basement  reflected  waves  are  designated 

BR( 1 ) =VV 

BR(P)=-VV 

BR( 3  )  =  -VV 

BR(4)=VV 


The  negative  signs  are  to  augment  phases  by  it  radians  for  1  surface  reflection 
and  the  positive  signs  account  for  zero  and  2ir  radian  phase  shifts  for  no 
surface  and  two  surface  reflections  respectively. 


'  b>3/) 

V= (BST( J)  +  ioc0BSLS( J))=(t+ioc0  S  )  is  used  to  form  the  complex  frequency 

dependent  phases  mentioned  in  the  Discussion  as  follows: 

PH  ( 1 )  =  tj?(  *  ZjT  (  t  +  ic^&5) 

PH  ( 2 )  =( pz=Zir(t+ts+'*°s) 

PH ( 3 )  =  -hiAeS) 

PH  (4)  =  4^-2  Tr{t+tR  +  ts  +  i<*-oS) 


Here  BST(J)  =  t  is  the  travel  time  in  water  plus  time  in  the  sediment,  tp 
and  tc  are  surface  scattering  travel  time  differences,  and  ALPH0*  BSLS(J) 
is  the  attenuation  loss  per  unit  frequency  in  the  sediment  =c£0S  where 
S  is  the  "effective"  path  length.  These  are  complex  phase  factors. 
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V  =  0  is  defined. 


Zo 


DO110  K=l,3  to  110  CONTINUE  produces  AUX=BR( K)*BR( L)  i.e.  cross  products  of 
the  amplitudes  of  the  transmitted  rays  for  all  permitted  combinations  and 
AUX1 =( PH( K)-C0NJG( PH(L) )*AIM 


=i  (01.-0*)  the  corresponding  complex  phase  differences. 
n  L 

IF  (REAL ( AUX1  )*F?<(-40)  GO  TO  110  causes  terms  like  to  be  ignored 

unless  greater  than  -40  (Note  that  10"^)  in  which  case 


V  =  V+BR(K)*BR(L) 


<*>£) 


C* iL. ) 


Here  we  have  terms  applicable  to  the  transmitted  rays  like  those  in  the 
formula  for  the  frequency  averaged  intensity  ratio  shown  in  the  Discussion 
of  this  subroutine.  Terms  of  this  kind  are  generated  and  summed. 

AUX  =  2*REAL(V)=2ReV  gives  the  effective  value  of  the  intensity  of  the  cross 
product  terms  for  the  four  contributing  paths. 


AUX2=  -2[2£t)<££ 

IF  AUX2*F2  =  -47T<XbSF2  Z.  -40  go  to  120  causes  terms  like  -417<<>SF2  to  be 


ignored  unless  greater  than  -40  in  which  case 

-4-rrf.oS  F\ 2.  -ATTt£.»S 

AUX  =  2ReV  +  4(  e _ -  g _ J 


— 4- tt<£>c  -S 


Cm) 


The  last  term  accounts  for  dissipation  loss  along  the  path  in  the  sediment 

iZ  .  2. 


for  four  rays  with  initial  intensities  of  the  same  amount,  ^BR(ljj)  =(VV)  , 
combined  incoherently. 

r  „  /  -4 

|2ReVf  4l.e  -e 


~47T£Ct,5‘^  j-ATTctcSF/ 


D  T—  PT  4- 


fi 


■4.tt<<o  s 


)2.‘ 


C  OS 


(&8) 


R 
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On  the  right-hand  side  of  the  last  equation  DT  represents  the  frequency 
averaged  effects  of  the  layer-reflected  rays  and  in  the  second  term  2ReV  gives 
the  effects  of  the  cross  product  terms  for  the  transmitted  rays  while  the  last 

term  gives  the  effects  of  the  frequency  averaged  attenuation  in  the  sediment 

2. 

on  the  4  rays  with  intensities  (VV)  .  In  the  denominators  the  factors  FI  and 

2. 

4tto(6S  are  needed  for  the  frequency  averaging  and  the  4  (cos0/R)  is  the  geo¬ 
metric  intensity  sum  for  4  rays  with  R/cos©as  the  slant  range.  The  last  step 
gives  the  bottom  loss  for  angles  above  the  caustic: 

ALOSS ( JJ,J)=-10  logJO(DT). 

CONTINUE 

NMAX=91-N  is  definition. 

COMPUTE  LOSS  FOR  ANGLES  BETWEEN  TIR  AND  THS . 

This  applies  to  the  domain  designated  by  Rj  in  figure  2,  i.e.  for  angles 

between  THS=SAP  the  smallest  bottom  grazing  angle  or  TIR  the  angle  of  total 

internal  reflection,  whichever  is  greater,  and  the  basement  limiting  angle. 

Previously  we  had  N=90-NTHC 

NMAX=91-90+NTHC=l+NTHC 

NTIR=TIR+1  gives  angle  in  degrees. 

N=NMAX-NTIR 

N=1+NTHC-NTIR-1=NTHC-NTIR  redefines  N 
DO  200  I  =  1 ,N=1 , NTHC-NTI R 

J=NMAX-1=NTHC  is  maximum  value  of  J  while  the  minimum  value  is  1+NTIR  or  NSAP. 
IF  (J<NSAP)  go  to  800  causes  computation  for  angles  between  TIR  and  THS  as 
stated.  Otherwise,  for  rays  in  domain  R2  shown  in  figure  2,  the  following 
is  done: 
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COMPUTE  REFLECTION  COEFFICIENT 
COMPUTE  FREQ  AVERAGE  ABSOLUTE  REFLECTION 

From  here  to  the  statement  VV=AUX/NAV1  the  frequency  average/!  squared  reflec¬ 
tion  coefficient  VV  is  computed  as  for  the  case  of  angles  above  the  caustic 
in  this  subroutine. 

ALPH2  =  FPI2D*RSIN2(  J)=  4ffDF(  JJ)sinfi>  = 


frequency  to  find  a  new  V. 


Compute 


C=cos 
S=sin  <^2 

V=RV21+RV32cosoC2.+i  RV32sinQC'2 _ 

( /  -A  EV2 /  ^32  rp^a^- 1  7  **£2)  •»  /  Cos  &,)  (33ct  ) 

I  ±  [ 

I  |  f-  RUZl  RVjZ.Ccs^i+I  RVZl  ^  “^2.  / 

Here  the  first  term  contains  the  new  calculated  for  the  band  center  fre¬ 
quency,  and  it  gives  proportional  values  of  real  and  imaginary  parts  of  the 
complex  layer  reflection  coefficient  while  the  second  term,  (f VV J  gives  the 
magnitude,  and  the. last  term  is  the  geometric  amplitude  factor.  Thus  V  is 
the  complex  amplitude  of  the  layer  reflected  ray.  This  is  used  to  obtain 
a  coherent  sum  of  contributions  from  the  4  paths  for  this  angular  domain. 


These  give  amplitudes  of  the  four  layer  reflected 
paths  including  effects  of  phase  shifts  upon  surface 


reflection. 
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This  applies  to  path  R2.C.  of  figure  2. 

AR(5)=0 

AR(6)=0  i  Sets  array  =  0 

AR(7)=0 

AR(8)=0  J 

The  layer  transmission  coefficient  factor  BTRF  for  the  Big  Angle  ray  is 
needed  as  computed  in  SUBROUTINE  RE FR . 


IF  BTRF=0  GO  TO  205  Skips  Big  Angle  computation 

IF  BT(S)  RT(J)  GO  TO  205/  for  physically  unrealistic  cases. 

Here  BT  is  total  travel  time  and  RT  is  travel  time  in  water. 

ALPH2=  ^3=  4ifDF( JJ)  si nd^.  where  Op  is  the  angle  in  the  layer  for  the  big 

cz  '  e 

angle  case  and  <^2.  is  The  phase  factor  for  transmission  one  way  through  the 
1 ayer . 


V  = 


f  A 2-  'f’ 

j  c  os  $$ 

/J/-fBv/32BV2./ )coscCx  4-  \  (BVmVZI  -l)si 


This  is  the  two-way  transmission  coefficient  for  paths  refracted  in  the  sedi¬ 
ment  multiplied  by  the  geometric  spreading  factor. 

AR (5 ) =V  ^  These  are  amplitudes  for  sediment  refracted 

AR(6)=-V  >  paths  including  effects  of  surface  reflections 

AR(7)=-V  |  on  phases. 

AR(8)=V  j 

V=BT(J)+ioCeBLS(J)  >) 

V=t+i^CoSc  the  total  travel  time  for  path  with  bottom  grazing  angle  BA  =  Big^/^t/^ 
Angle  plus  dissipation  loss  per  unit  frequency  in  sediment. 
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PH  ( 5 )  =  2tfU+tfoSc) 

PH(6)  =  2»T(t+tstic<DSc) 

PH ( 7 )  =  2T(t+tRti^DSc) 

PH  (8)  =  2rT(t+ts+tR+idCoSc) 

These  are  the  complex  phase  factors  including  effects  of  surface  scattering 
travel  times  and  dissipation  loss  factors  for  this  case. 

CONTINUE 

COMPUTE  SMALL  TRANSMISSION  COEFFICIENT 

This  applies  to  path  R2.b.  of  figure  2.  Here  the  computation  is  exactly  the 
same  as  for  the  Big  Angle  case  but  using  appropriate  factors  for  the  Small 
Angle  case.  One  exception  occurs.  The  transmission  coefficient  is  made 
negative  imaginary:  _ _ 

V=-V*SGI(J)*AIM  (m~>) 

This  should  account  for  a  (controversial)  phase  change  due  to  the  caustic  in 
the  sediment  shown  in  figure  2  for  ray  R2.b.  and  mentioned  on  page  56  of 
Spofford  et  al . ( 1982). 

The  amplitudes  and  phases  are  computed  as  in  the  Big  Angle  case.' 

COMPUTE  BASEMENT  RAY  TRANSMISSION  COEFFICIENT. 

This  applies  to  path  R2.d.  of  figure  2. 

The  same  procedure  is  used  as  for  the  Big  and  Small  Angle  coefficients.  We 
get  the  amplitudes  of  the  transmitted  rays  by  multiplying  the  transmission 
coefficients  by  the  geometric  spreading  factor  and  the  basement  reflection 
coefficient  REF.  This  accounts  for  ray  R2.d.  in  figure  2. 

Amplitudes  and  phase  factors  are  computed  as  before  in  this  subroutine. 

The  rest  of  the  computations  are  done  in  a  similar  manner  as  before  except 
for  differences  in  notation. 


G3V3  ) 
CWJ) 
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DO  210,K*1 ,15 
NA=K+1 


DO  210,L=NA,16 

9 

AUX  99=AR(K)*CONJG(AR(L) ) 

=AR(K)  AR  (L)  computes  cross  products  of  complex  amplitudes  for  all 
three  cases  of  angles  from  caustic  to  basement  limiting. 


IF  |AUX99|  =0  GO  TO  210 

AUX  l=i(<£>  -  <p  )  computes  complex  phase  differences. 

IF  REAL  (AUXl )*F2  4 -40  GO  TO  210  causes  terms  like  ^TT^SF^  to  be  ignored 


unless  greater  than  -40  in  which  case  the  sum  is  computed. 

N  „  ,  JFzWk-cP?)  if,  (4>k-4>1) 
V  =v  +ARt.K)AR*(i)  £ _ _ 

'  (  4k  -  4>*) 


This  makes  a  sum  of  contributions  from  the  Small  Angle  path.  Big  Angle  path 
and  basement  reflected  path. 


CONTINUE  to  209. 

AUX=2*REAL(V)  gives  the  average  value  of  a  complex  quantity  representing  the 
•sum  of  contributions  from  the  four  paths. 

IF  /AR(13)/  =0  go  to  209  Skip  next  computation  unless  transmitted  contribution 
of  basement  ray  intensity^  0. 

AUX  2=-4  7 T(^0S^  is  effective  path  length  for  basement  reflected  ray  in 
sediment. 

IF  -4 TT^qS^  Fg  -40  GO  TO  209  Ignore  contributions< -40 . 


AUX=2ReV+(4) 


f2_  e-4‘/reC^F/ 
-  4rr<c ©  % 


J 
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'  gives  contributions  for  basement  reflected  paths. 


AUX2=-4TE(^S^ ,  S^is  effective  path  length  in  sediment  for  Big  Angle  path. 


IF  -47T^  S  F,  /  -40  TO  TO  211 

o  c  2^ 

IF  AR(5)  =  0  GO  TO  211 


AUX 


2ReV  +-  («)£- 


Ignore  contributions  too  small  to  be 

4  TToC  _  Cj  P, 


'  physically  realistic. 
r4TTc*  oSd  Fz  ^4  F, 


(<4 


-  A-rrdL* 

T-s  c  P.  — Fl 


(3SI) 


—  Af-  rr  e>  ^5  £ 


adds  contributions  for  ray  along  Big  Angle  path. 
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AUX2=  -4frK<?-^ }  is  effective  path  length  in  sediment  for  Small  Angle  path. 
IF  -4.fr *0Sh  f^-^D GO  TO  212.  7  Ignore  contributions  that  are  not 

Cm) 


AUX=2ReV- 


IF  AR( 9)  =0  GO  TO  212. 

-477-cCc^^  -47TcC0-%  ^ 

r+( 4)  £ - -  e  01  ' 

'-4^0^ 

—  4-  TT&C  o  *^c-  •—  A/-TTc£  F, 

—  — - — - 


physically  realistic. 


V-C4-) 


—  4tt<£c>S  c- 

-4-ir<*.oSk>'Fz  -J-JT^o^bn 

k _ - _ ~  S 

-  4rr^ *  5 *> 


adds  contributions  for  ray  along  Small  Angle  path. 

AUX=AUX  +  4  /  AR(1)/2 
F/ 


G^vj 


The  first  term,  AUX  is  needed  as  part  of  the  frequency  averaging  and  the 

n 

second  term  4j[AR(l)|  adds  the  sum  of  intensities  of  the  4  layer  reflected  rays. 
IF  AUX>0  GO  TO  213  where  the  computation  is  completed.  Otherwise  set  the  value 
AL0SS( JJ , J ) =25 .  This  is  an  assumed  value. 

GO  TO  200  allows  continuation  of  computations. 


AL0SS( JJ , J ) =- 10  log  _  AUX 

10  4(cos0/R)‘z- 


(3SSj 


This  completes  the  computations  for  the  domain  where  Big  Angle,  Small  Angle, 
and  basement  reflected  paths  can  all  contribute.  The  4  in  the  denominator 
represents  the  sum  of  4  unit  intensities  and  (cos0/R)  is  the  inverse  squared 
slant  range  factor. 

CONTINUE 
N=J- 1 
NA=  J 

DO  701  1=1, N 
J=NA- I 
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IF  «KNSAP  GO  TO  999  Stops  computations  for  bottom  angles  too  small  for  ray  to 
travel  from  source  or  to  receiver. 

Statements  from  DR  to  last  DT  compute  same  quantities  as  in  case  for  angles 
above  critical  (caustic)  near  the  beginning  of  this  subroutine,  but  for  angles 
greater  than  NSAP. 

ALOSS ( JJ, J )  =- 1 0  log^DT  gives  loss  for  angles  between  SAP  and  TIR.  For  this 
case,  total  internal  reflection  occurs,  no  paths  penetrate  the  layer,  and  it 
is  assumed  no  losses  occur  in  the  sediment. 

CONTINUE 

GO  TO  999  Stops  the.  subroutine.  You  will  remember 

IF  J<NSAP  the  routine  goes  to  800. 

N= J- 1 
NA=J 

DO  801  1=1, N 
J=NA- I 

IF  J<NTIR  GO  TO  999  Computations  are  done  for  angles  between  TIR  and  basement 

limiting  ray  angle,  e.g.  ray  Rl.b.  shown  in  figure  2 

COMPUTE  BIG  TRANSMISSION  COEFFICIENT 

ALOSS ( JJ,J)=25  Loss  value  is  assumed. 

ALPH2=°G)  =  4rrDF( JJ)sinBA  Uses  Big  Angle. 

CZ 

V=V*BGI(J)  gives  the  two-way  transmission  coefficient  multiplied  by  the 
geometric  amplitude  as  done  before  using  Big  Angle. 

AUX3 =(vlc  gives  intensity  of  two-way  transmitted  ray. 

IF  AUX3=0  GO  TO  801  Skip  rest  of  computation. 

Otherwise  compute 
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DR  to  last  DT  statement  as  done  before  in  this  subroutine. 

AUX2=-4ftV(,(>S  log^S  where  S  represents  the  effective  Big  Angle  path  in 
the  sediment. 

IF  AUX2*F2  <-40  GO  TO  801.  causes  this  calculation  to  be  skipped  if 


-47FK0Slog10e<f-40 
Otherwise  compute 
AL0SS( JJ ,  J) =-10 


if  DT  1 1 

-VTTAoS  lo$loa  -h  lc,<5/oj~jco£pjzj  I  £ 


=  ~!0  /°3/( 


r 


lVl2(vr)e 


Here  only  sediment  refracted  paths  are  considered--no  layer  reflection. 


IF  AL0SS(JJ,J)>25  AL0SS( JJ, J)=25 


This  replaces  the  computed  value  when  greater  than  25  with  the  value  25  on 
the  assumption  that  the  greater  value  is  unrealistic. 

CONTINUE  Computations  are  made  in  steps  of  1  degree  for  the  allowable 

val ues . 

RETURN 

END  Stops  SUBROUTINE  LOSS(JJ) 


71 


III.  ANALYSIS  OF  DATA  AND  SELECTION  OF  VALUES  FOR  INPUT  TO  COMPUTER 


Program  BTLS  has  two  input  options.  The  choice  of  the  option  desired  provides 
for  the  computation  of  either  Geo-parameters  or  simulated  bottom  loss.  In 
either  case  the  environmental  geometry  is  important  and  the  following  inputs 
are  required. 

1.  T  =  thickness  of  the  sediment  (from  the  bottom  of  the  water  to  the  inter¬ 
face  with  the  basement)  is  characterized,  not  by  units  of  length,  but  by  the 
two-way  travel  time  for  sound  to  travel  from  top  of  the  sediment  to  the  base¬ 
ment  and  return.  This  time  should  be  entered  in  units  of  tenths  of  seconds. 
Two  decimal  places  should  be  used  (e.g.,  if  the  two-way  travel  time  is,  say, 
1.7  seconds,  it  would  be  entered  as  17.00.).  The  time  may  be  obtained  from 
sediment  thickness  charts  or  other  source. 

2.  NNFILE  =  Arbitrary  file  number--a  number  which  designates  the  geographi¬ 
cal  site  of  the  measurements  or  predictions. 

3.  IOPT  =  input  option,  0  for  Geo-Parameters 

1  for  Inversion  Parameters 

IOPT =0  causes  computation  (and  plotting)  of  simulated  bottom  loss. 

I0PT=1  causes  computation  and  printout  of  geo-parameters. 

4.  The  data  listed  next  is  required  for  all  computations,  and  it  should  not 
be  controversial  in  nature. 

ZB  =  water  depth  (of  bottom  below  surface).  Units  are  meters.  Entry 
should  contain  two  decimal  places,  e.g.  ZB=4800.00. 
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ZS  =  depth  of  sound  source  below  surface  of  water  (Standard  depth  is  244 
meters.)  Units  are  meters.  Entry  should  contain  two  decimal  places, 
e.g.  ZS  =  244.00. 

ZR  =  depth  of  receiver  of  sound  below  surface  of  water.  Units  are  meters. 

Entry  should  contain  two  decimal  places,  e.g.  ZR=305.00. 

CS  =  Speed  of  sound  in  water  at  the  depth  ZR  or  ZS,  whichever  is  smaller. 
Units  are  meters/second.  Entry  should  contain  two  decimal  places, 
e.g.  CS=1 507 . 00 . 

Cl  =  Speed  of  sound  in  water  at  bottom  of  water  column.  Units  are  meters/ 

second.  Entry  should  contain  two  decimal  places,  e.g.  C 1  =  1 54 3.70. 


GEOPHYSICAL  PARAMETER  COMPUTATION 


For  I  OPT  =  1  the  INPUT  INVERSION  PARAMETERS  required  are  listed  below. 
INPUT  is  format  free.  The  following  are  required  inputs. 


1.  DBL0SS 


8.  THC 


2.  FREQ 


9.  RATIO 


3.  THETA 


10.  G 


4.  ALOSS0 


11.  BETA 


5.  AL0SSM 


12.  GAM 


6.  FM 


13.  REF 


7.  TIR 
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The  statements  which  follow  are  intended  to  be  used  for  guidance  in  determining 
the  rest  of  the  necessary  input  data  for  computer  runs  of  Program  BTLS.  Data, 
tables,  and  articles  referred  to  are  meant  to  be  used  as  examples  only.  They 
may  not  represent  true  values  for  a  given  situation.  As  bottom  loss  work  pro¬ 
gresses,  it  is  expected  that  much  data  will  accumulate,  and  the  techniques 
and  "art"  used  in  interpretation  will  become  highly  sophisticated.  Therefore, 
it  is  essential  to  use  all  evidence  at  one's  disposal,  and  take  great  care  in 
choosing  assigned  or  estimated  values.  Otherwise,  large  errors  in  output 
results  should  be  expected.  With  this  caveat  in  mind  we  proceed  to  a  discussion 
of  ways  and  means  to  arrive  at  input  data. 

For  the  measurement  location,  the  bulk  sediment  type  is  to  be  obtained  from 
measured  core  data,  core  data  summaries,  or  other  source.  The  type  of  sedi¬ 
ment  must  be  known  in  order  to  choose  properly  some  assigned  values. 

The  quantities  DBLOSS,  FREQ,  and  THETA  are  used  to  compute  ALPH0,  the  logari¬ 
thmic  attenuation  coefficient  of  the  sediment  in  dB/m/kHz. 

DBLOSS  is  a  value  found  by  inspection  of  a  graph  (or  possibly  a  table  of 
data)  of  measured  bottom  loss  as  a  function  of  apparent  grazing  angle. 

Consider  Figure  18  as  an  example.  Look  on  the  curve  for  the  frequency  400  Hz. 
Find  a  reasonably  smooth  portion  of  the  curve  which  seems  to  be  as  nearly 
linear  as  possible.  Read  the  graph  and  determine  the  bottom  loss  in  DB  and 
the  apparent  bottom  grazing  angle  for  a  point  on  the  smooth  linear  portion 
of  the  curve.  For  example  DBLOSS  =  3.5  dB,  FREQ  =  400  Hz,  THETA  =  20  degrees. 
Entry  should  contain  two  decimal  places,  e.g.  DBLOSS  =  3.50. 
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For  other  examples  consider  Figure  19  where  DBLOSS,  =  4dB,  FREQ  =  400  Hz, 

THETA  =  20°  while  Figure  20  indicates  DBLOSS  =  8.5  dB,  FREQ  =  400  Hz, 

THETA  =  20°.  Note  that  usually  one  should  choose  mid-frequencies  and 
modest  grazing  angles  to  obtain  values  of  DBLOSS,  FREQ,  and  THETA^ 

ALOSS0,  ALOSSM,  and  FM  are  used  to  compute  the  density  of  the  layer  above 
the  sediment,  (2.=  RH2,  the  density  of  the  sediment,  =  RH3,  and  the  thick¬ 
ness  D  of  the  layer.  To  find  these  input  quantities  look  at  the  bottom  loss 
curves.  Look  at  the  portions  of  the  curves  for  large  angles.  Find  the  fre¬ 
quency  at  which  the  loss  is  smallest.  This  frequency  is  the  value  FM.  The 
units  to  be  used  for  this  quantity  are  Hertz  (Hz).  Now  read  the  minimum 
value  of  bottom  loss  for  the  curve  just  selected.  This  minimum  value  of 
bottom  loss  is  ALOSSM  in  units  of  dB  and  usually  occurs  at  relatively  high 
frequency.  This  quantity  is  characteristic  of  the  thin  layer  and  of  frequen¬ 
cies  greater  than  300  to  400  Hz,  therefore  data  for  frequencies  lower  than 
300-400  Hz  should  not  be  used  here.  The  value  of  ALOSS0  is  the  amount  of 
bottom  loss  where  the  curve  has  a  maximum  at  low  frequency  (say,  less  than 
300-400  Hz).  These  entries  should  all  contain  two  decimal  places,  e.g. 

ALOSS0  =  9.00,  ALOSSM  =  7.00,  FM  =  1600.00. 

The  angle  of  total  internal  reflection,  TIR,  is  the  largest  grazing  angle 
where  no  sound  energy  penetrates  into  the  layer  or  sediment.  All  the  sound 
is  reflected  from  the  bottom  for  bottom  grazing  angles  less  than  this  angle. 
This  phenomenon  can  occur  for  sediments  (or  the  top  layer)  in  which  the 
speed  of  sound  is  greater  than  the  speed  of  sound  at  the  bottom  of  the  water 
column.  For  the  angle  TIR  and  smaller  angles  there  is  no  bottom  loss. 
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HOT  T  OM  I.OTVS  iflOl 


Spofford,  Greene,  and  Hersey 
J.  Ac  oust.  Soc.  Am. 


Figure  20.  Bottom  loss  versus  grazing  angle  measured  in  six  one-octave 
bands  at  indicated  center  frequencies:  negative  high-angle 
frequency  dependence. 


After  Spofford  et  a1.(1982) 
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Consequently,  you  should  inspect  the  bottom  loss  curves  and  find  the  greatest 
angle  for  which  there  is  a  zero  value  of  bottom  loss.  This  angle  is  TIR  and 
the  units  are  degrees.  Entry  should  have  two  decimal  places,  e.g.  TIR  =10.00. 

The  caustic  angle,  THC,  is  the  apparent  bottom  grazing  angle,  at  which  there 
is  a  marked  increase  in  bottom  loss.  For  example,  in  Figure  18  the  caustic 
angle  THC  =  40°  as  indicated  by  spikes  in  the  curves.  Notice  Figure  19;  no 
caustic  effect  can  be  found  here.  The  sediments  are  not  thick  enough  to 
allow  the  caustic  to  form.  The  basement  grazing  ray  appears  to  correspond 
to  the  bottom  grazing  angle  of  about  25  degrees.  This  figure  is  discussed 
further  under  High  Angle  Losses,  Type  2.  In  Figure  20  the  caustic  angle  is 
at  34°.  Use  degrees  as  the  proper  units  for  THC.  This  angle  is  interpreted 
as  the  boundary  separating  strong  sediment  refracting  paths  from  weak 
reflected  paths.  Entry  should  have  two  decimal  places,  e.g.  THC  =  39.00. 

The  quantity  RATIO  is  the  value  of  the  speed  of  sound  in  the  top  layer  of 
the  sediment,  C2,  divided  by  the  speed  of  sound  in  the  water  at  the  bottom 
of  the  water  column.  It  is  computed  by  use  of  Snell's  law  from  the  angle 
TIR.  Alternately,  if  desired,  it  can  be  an  assigned  value.  It  is  a  number 
without  units.  Hamilton's  value  of  sound-speed  ratio  by  sediment  type  from 
Table  1,  reproduced  here  from  Spofford  et  al.(1982),  has  been  used  directly 
and  given  satisfactory  results  in  some  cases.  In  some  instances  it  should 
be  adjusted  to  account  for  unusually  low  or  high-speed  bottoms  which  manifest 
themselves  in  the  data  by  either  a  finite  loss  at  zero  or  minimal  grazing 
angle,  or  zero  loss  at  a  finite  grazing  angle,  respectively.  The  entry 
should  have  four  decimal  places,  e.g.  RATIO  3  1.0040. 
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Table  I.  Sound  Speed  profile  parameters  from  Hamilton 
by  sediment  type  and  corresponding  values  of  6  for  c(z) 
given  by  Equation  (15).  (units  -  km,  sec) 
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After  Spofford  et  al.(19H2) 


.  The  gradient  of  the  speed  of  sound  in  the  sediment  at  the  top  of  the  sediment 
is  designated  G.  This  quantity  is  computed  from  input  data  including  the 
apparent  caustic  grazing  angle,  THC,  especially.  If  desired,  an  assigned 
value  can  be  used  for  input  data.  Units  are  1/second.  Hamilton's 
tabulated  values  are  suggested  as  a  source  of  numerical  values  when  it  is 
desired  to  assign  a  value.  Table  1  may  be  used  as  a  source  of  possible 
values  with  the  caution  that  these  values  are  to  be  considered  as  examples 
only.  Entry  should  have  two  decimal  places,  e.g.  G  =  1.79. 

The  parameter  BETA=B  controls  the  rate  of  decrease  of  the  sound  speed  gradient 
in  the  sediment  with  increasing  depth  and  increasing  sound  speed.  It  controls 
the  shape  of  the  sound  speed  profile.  Its  values  lie  between  -1  and  +  in¬ 
finity.  This  parameter  is  a  pure  number  with  no  units.  The  value  of  B  must 
be  assigned.  Some  typical  values  of  B  are  given  in  Table  1.  As  illustrated 
in  the  table,  the  value  of  B  is  based  upon  the  type  of  sediment  at  the  site 
under  consideration.  Values  given  in  Table  1  are  typical  values  only  and 
may  not  be  characteristic  of  a  particular  site  (NNFILE  number).  Entry  should 
show  a  sign  and  two  decimal  places,  e.g.  BET=-.50. 

The  quantity  GAM=)fis  the  gradient  of  the  logarithmic  attenuation  coefficient 
of  the  sediment  and  has  units  dB/m^/kHz.  The  value  of  this  quantity  must 
be  assigned.  The  value  of  GAM  ranges  from  0  to  0.00040  and,  as  mentioned  in 
"Documentation  of  Bottom  Loss  Upgrade  Parameters,"  increments  may  be  typically 
0.00005.  Roughly  the  value  GAM  =  0.00005  corresponds  to  terrigenous  abyssal 
plain  sediments,  GAM  =  0.00010  to  terrigenous  continental  rise  sediments,  and 
GAM  =  0.00020  to  calcareous  sediments.  The  non-linear  loss  at  higher  angles 
and  lower  frequencies  may  help  in  estimating  GAM.  The  entry  should  show  5 
decimal  places,  e.g.  GAM  =  0.00010. 
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The  basement  reflection  coefficient  REF  is  a  pure  number  used  to  represent 
the  sound  pressure  reflection  coefficient  for  sound  penetrating  the  sediment, 
reflecting  from  the  basement,  and  returning  through  the  sediment  to  the  water. 
The  value  must  be  assigned.  This  is  a  pure  number  with  no  units,  -and  its 
value  is  in  the  range  from  0  to  1.  Basement  reflected  paths  may  make  important 
contributions  to  the  sound  level  after  a  bottom  interaction  occurs.  Sometimes 
this  effect  is  observable  in  bottom  loss  curves.  For  example,  a  look  at 
Figure  19  shows  a  rise  in  the  curves  at  a  bottom  grazing  angle  of  28°.  This 
has  been  interpreted  as  the  effect  of  reflection  of  the  sound  from  the  base¬ 
ment,  the  onset  of  basement  reflection  taking  place  at  28°.  Basement  reflec¬ 
tion  coefficients,  REF,  of  0.4  to  0.5  independent  of  frequency  and  grazing 
angle  appear  to  be  quite  common.  The  entry  should  contain  two  decimal  places, 
e.g.  REF  =0.50. 


•High  angle  losses  are  discussed  by  Spofford  et  al .  (1982)  as  follows: 

D.  High  Angle  Losses 

At  the  higher  angles  the  bottom-loss  data  show  several  characteristic 
patterns .  Very  little  significant  angle  dependence  is  apparent;  however,  four 
types  of  frequency  dependence  have  been  identified : 

Type  1.  Zero  (0)  Frequency  Dependence 

This  is  the  simplest  behavior  and  represents  the  frequency  and 
angle-independent  loss  characteristic  of  a  two- fluid  interface .  Figure  18  is 
a  good  example  of  this  type  where  the  loss  of  approximately  10  dB  corresponds 
to  an  impedance  of  3.2  or  a  sediment  density  of  2.0  g/cm^ ,  characteristic  of 
sandy-silt  found  on  continental  terraces. 

Type  2.  Positive  (+)  Frequency  Dependence 

When  the  sediments  are  not  so  thick,  at  the  high  angles  the  interface- 
reflected  path  is  augmented  by  a  basement-reflected  path  which  can  carry  signi¬ 
ficant  energy,  especially  at  the  lower  frequencies  where  volume  attenuation  in 
the  sediment  is  diminished .  Figure  19  illustrates  such  a  case.  The  sediments 
are  not  thick  enough  to  allow  the  caustic  to  form.  In  fact  the  basement 
grazing  ray  appears  to  correspond  to  the  bottom-grazing  angle  of  about  25 
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degrees  consistent  with  the  known  sediment  thickness  at  this  location  of 
approximately  90  meters  and  the  gradient  appropriate  to  this  type  of  calcareous 
sediment  of  1.7  sec”1.  The  higher  frequency  losses ,  dominated  by  the  interface- 
reflected  wave ,  tend  to  cluster  around  10  dB.  As  the  frequency  decreases ,  the 
amplitude  of  the  basement-reflected  path  increases  until  the  two  paths  have 
comparable  amplitude  at  200  Hz  (thereby  reducing  the  loss  by  3  dB)  .  At  50  and 
100  Hz  the  basement  return  dominates  and  attenuation  effects  are  nearly  negli¬ 
gible.  Basement  reflection  coefficients  of  6-8  dB  corresponding  to  REF  = 

0.4-0. 5  independent  of  frequency  and  grazing  angle  appear  to  be  quite  common. 

Type  3.  Negative  (-)  Frequency  Dependence 

A  third  type  of  frequency  dependence  at  high  angles  is  illustrated 
in  Figure  20.  Here  the  high  frequencies  show  less  loss  than  the  low  frequen¬ 
cies .  This  appears  to  be  modelable  as  a  thin  high- impedance  layer  at  the  top 
of  the  sediment.  The  shorter  wavelengths  respond  to  this  impedance  while 
longer  wavelengths  tend  to  "see"  through  the  layer ,  responding  more  strongly 
to  the  lower  impedance  characteristic  of  the  bulk  of  the  sediment  as  the  fre- 
queny  decreases.  In  this  case,  it  would  appear  that  the  thin  layer  has  an 
effective  impedance  ratio  of  1.50,  whereas  the  bulk  of  the  sediments  has  a 
near-surface  ratio  of  1.29.  The  higher  impedance  is  characteristic  of  sandy 
silt,  whereas  the  bulk  of  the  sediment  might  be  clayey  silt. 

Type  4.  Reversing  (±)  Frequency  Dependence 

A  fourth  type  of  behavior  is  occasionally  observed  when  the  thin-layer  (Type  3) 
effect  appears  to  peak  at  a  particular  frequency .  In  Figure  21  the  loss  de¬ 
creases  at  high  angles  until  200  Hz,  where  it  reaches  a  minimum  and  then  in¬ 
creases.  This  strongly  suggests  a  high-impedance  layer  (or  set  of  layers) 
which  acoustically  resonates  at  200  Hz.  If,  in  fact,  this  is  a  single  layer 
it  would  then  have  to  be  about  2  meters  thick  with  a  well  defined  lower  inter¬ 
face  with  the  sediment. 


RESULTS  OF  GEO-PARAMETER  COMPUTATION 

The  inversion  computation  yields  results  as  follows: 

RATIO  =  C2 
Cl 

RH2  =  density  of  top  layer  of  sediment  in  gram/cm3. 

RH3  =  density  of  sediment  below  the  layer  in  gram/cm3. 

D  =  thickness  of  the  top  layer  in  meters. 

G  =  sound  speed  gradient  at  the  top  of  the  sediment  in  units  sec"1. 
B  =  sound  speed  profile  parameter,  no  units. 
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ALPH0  =  logarithmic  sound  attenuation  coefficient  for  sound  in  sediment, 
units  DB/m/kHz. 

GAM  =  attenuation  coefficient  gradient  in  units  of  dB/m2/kHz. 

BZ  =  depth  of  the  sediment  in  meters. 

REF  =  basement  reflection  coefficient,  no  units. 

BOTTOM  LOSS  COMPUTATION 

If  bottom  loss  is  the  desired  result,  use  input  option  statement. 

IOPT  =  0. 

The  required  input  data  is  listed  below. 

This  example  illustrates  the  number  of  decimal  places  needed  for  each  entry: 


1. 

RATIO  =  1.0040 

6. 

D  =  .24 

2. 

G  =  1.79 

7. 

ALPH0  =  .029 

3. 

BETA  =  -.50 

8. 

GAM  =  .00010 

4. 

RH2  =  2.33 

9. 

REF  =  .50 

5. 

RH3  =  2.09 

10. 

BZ  =  1067.78 

These  quantities  mentioned  before  have  the  same  definitions  given  before. 

The  computation  yields  values  of  bottom  loss  in  dB  in  steps  of  1  degree  for 
all  applicable  apparent  bottom  grazing  angles  from  0  to  90°.  Computations 
are  done  for  the  following  frequencies  (in  Hz):  50,  100,  200,  400,  800,  and 
1600,  or  others  as  specified.  Frequencies  mentioned  are  the  frequencies  of 
the  centers  of  bands  one-octave  wide,  and  the  computed  bottom  loss  is  a 
frequency-averaged  value  for  the  one  octave  bandwidth.  The  results  are 
plotted  in  graphical  form. 
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For  runs  involving  a  basement  reflection,  if  90  degrees  is  reached  on  the 
loss  plot,  the  loss  value  increases  to  the  value  one  would  have  if  no  base¬ 
ment  return  were  present.  This  may  be  helpful  to  the  user  in  evaluating  the 
effects  of  basement  reflection  and  in  adjusting  the  value  of  the  basement 
reflection  coefficient  REF. 

Computed  bottom  loss  values  of  25dB  should  be  treated  suspiciously  because 
they  may  be  the  result  of  computations  which  result  in  loss  values  of  greater 
than  25dB.  These  are  thought  to  be  physically  unrealistic  so  are  set  arbi¬ 
trarily  at  25dB  by  the  computer.  See  statement  in  Subroutine  LOSS: 
AL0SS(JJ,I)=25- 

See  also  similar  statements  between  212  and  213  and  between  800  and  801  of 
Subroutine  LOSS. 
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IV.  CONCLUSIONS 


The  computer  program  directs  computations  of  geo-acoustic  ocean  bottom  para¬ 
meters  and  simulated  bottom  loss,  generally  in  accordance  with  the  criteria 
outlined  in  the  description  of  the  model  of  Spofford  et  al.(1982).  The  pro¬ 
gram  has  been  implemented  on  HP  1000  computers  inhouse  at  the  Naval  Oceano¬ 
graphic  Office  and  aboard  ship  at  sea.  At  present,  new  bottom  loss  data  is 
being  acquired  by  experiments  at  sea,  and  results  are  being  computed  by  use 
of  PROGRAM  BTLS.  Some  initial  problems  with  use  of  the  program  have  been 
overcome,  and  use  and  evaluation  are  continuing.  Some  problems  associated 
with  the  program  are  discussed  in  the  section  containing  recommendations. 
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V.  RECOMMENDATIONS 


It  is  recommended  that  changes  be  made  to  the  program  as  follows: 

In  SUBROUTINE  GEOPHYS  change  line  D=Cl/FM/4  to  read  D=C2/FM/4. 

Change  line  AUX=4*D*PI*FREQ/C1*$1N(THETA  to  read  AUX=4*D*PI*FREQ/C2*SI N(THETA) . 
These  change  the  speed  of  sound  in  the  layer  at  the  top  of  the  sediment  to 
the  correct  value  C2.  Change  AUX=AUX/T  to  read  AUX=AUX/T**2  in  order  to 
compute  the  square  of  the  two-way  layer  transmission  coefficient. 

Once  these  changes  are  carried  out,  the  program  should  be  run  again  with 
data  used  for  previous  computations  done  before  making  these  changes.  A 
careful  evaluation  of  the  old  and  new  results  should  be  made  to  determine 
the  effects  of  these  changes. 

The  program  computations  are  based  upon  an  assumed  bandwidth  of  1  octave. 

This  is  important  in  the  calculation  of  simulated  bottom  loss.  The  actual 
bandwidth  used  in  processing  measured  signals  to  provide  input  data  for  com¬ 
puting  geo-parameters  does  not  seem  to  be  important,  because  only  the  band 
center  frequency  is  designated,  and  the  bandwidth  does  not  enter  into  the 
calculations.  Nevertheless,  when  comparing  simulated  bottom  loss  plots  with 
actual  measured  bottom  loss  data,  the  bandwidth  used  for  the  simulation 
should  be  the  same  as  the  bandwidth  used  to  obtain  the  measured  data.  It 
seems  that  all  recent  measured  data  is  processed  with  a  1/3  octave  bandwidth. 
Therefore,  discrepancies  should  be  expected  in  a  comparison  of  simulated 
bottom  loss  with  1/3  octave  bandwidth  measured  loss.  It  is  recommended  that 
the  same  bandwidth  be  used  in  computing  simulated  bottom  loss  and  for  pro¬ 
cessing  measurements  of  bottom  loss. 
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APPENDIX  A 


-PROGRAM  BTLS 
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APPENDIX  B 


Flow  Chart  of  PROGRAM  BTLS 


/ 


BEGIN 


I 

I 


/ 


L 

PROGRAM  BTLS9 


I 

I - C  THIS  PROGRAM  SIMULATES  MEASURED] 

I  [BOTTOM  LOSS  DATA  ] 

I  C  INCLUDING  MEASUREMENT  ARTIFACTS] 

I 

I — - 1  BLANK  COMMON  CONTAINS  GEOACOUST] 
I  CC  PARAMETERS  ] 

I 
I 


COMMON  ZB,ZR,ZS,C1,C2,C3,RH1,RH2,RH3,G,D, 
alphO , cs , bz  ,  rf  e 

COMMON  /ALPH/GAM 
COMMON  /T/T 


I 

I---C  COMMON  PLOCK  AL  CONTAINS  The  ] 
I  [SIX  FREQUENCIES  DEFINED  ] 

I  [  IN  THE  FOLLOWING  DATA  STATEMENT] 

I  [IN  F,  AND  CORRESPONDING  1 

I  [  BOTTOM  LOSS  FOR  THE  SIX  ] 

I  [FREQUENCIES  AND  $0  ANGLES  IN  ] 

I  C  ALOSS.  THE  VALUSE  IN  ALOSS  ARE  } 

I  [CALCULATED  IN  SUBROUTINE  LOSS.  ] 

I 

I - [  ALL  VALUES  OF  ALOSS  A°F  ] 

I  [RETAINED,  WHICH  IS  ] 

I  [  NOT  NECESSARY  IN  THIS  VERSION  ] 

I  [OF  THE  PPO GRAM,  ] 

I 
I 


COMMON/A  1/ F (6) , ALOSS U ,90) 
DIMENSION  BRAYC546) 
EQUIVALENCE  (PRAY(1).F(1)),(BRAY(7), 
ALOSSdtl)) 


I 

I - [  LOOP  50P  READS  A  SET  OF  INPUT 

1  [DATA  AND  GENERATES  A  SET 
I  [OF  BOTTOM  LOSS  DATA  ON  EACH 
I  [PASS  UP  TO  A  MAX  1 M U M 
I 

I - [  SUBROUTINE  LCS2  CALCULATES  ALL 

I  [FREQUENCFY  INDEPENDENT 
I  r  INTERMEDIATE  VARIABLES 
I 
I 


1 

1 


1 

-) 


132 


•  • 

•  • 

•  • 

CALL  LOS2 

•  » 

•  » 

•  » 

I 

I 

I"4 

. \  TRUE 

< 

IF(T.EQ.999)  GO  TO  999 

> - 

\  .  .  . 

7 

I  FALSE 

7 

I 

7 

I - C  LOOP  41 

CALCULATES  BOTTOM  LOSS 

3 

7 

I  CAT  A  PARTICULAR 

3 

7 

I  C  FREQUENCY  IN  SUBROUTINE  LOSS 

3 

7 

I  CAND  PRINTS  IT 

3 

7 

I  C  OUT  ON 

EACH  PASS. 

3 

7 

I 

7 

I 

7 

7 

A  — 

—  > : 

DO  41  J=1 ,6 

7 

A 

• 

• 

7 

A 

I 

7 

A 

I 

7 

A 

7 

A 

;  : 

CALL  LOSS(J) 

»  • 

•  • 

Z 

A 

• 

•  • 

7 

A 

I 

7 

A 

I 

7 

A 

r4 

. / 

7 

A 

/ 

yPITE(6,2)(AL0SS(J,I),I  =  1 

,90)  / 

7 

A 

l . .  .  a 

. / 

7 

A 

I 

7 

A 

I - [  2  FORMATdx  ,8F8.2) 

] 

7 

A 

I 

Z 

A 

r 

41]  I 

7 

A 

7 

CONTINUE 

: 

7 

• 

. 

7 

I 

Z 

I 

7 

7 

2  \ 

CALL  E X E C( 1 4 , 2 , bRAY ,1092) 

:  : 

7 

•  • 

• 

7 

I 

7 

I 

7 

/ " 

— / 

7 

/ 

WRITE  (1 ,69) 

/ 

7 

/ . . . . 

7 

C  999] 


I  7 

I - [69  FORMAT ('  NORMAL  TERMINATION  -1  7 

I  [  BOTTOM  LOSS  MODEL')  3  7 

I  7 

0  < - o 

I 


STOP 
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\  SUBROUTINE  LOS2  / 

\ . / 

I 

I 


COMMON  ZB,ZR,ZS .C1.C2,C3,RH1,RH2,RH3,G,D , 
ALPHO.CS  ,BZ,REF 

COMMON  /DATA/  LZERO,  F  ZE R 0 , L 70 , F 7 0 , T H70 , 
L71 ,F71 .DBLOSS.FR E0» THETA,  A LOSS 0 , ALOSSM , 
FP ,TI R ,THC ,THRC  ,TH1B  ,B  ETA 


I 

I 


COUPON  /ALPH/GAP 
CO*MCN/T/T 
REAL  LZERO, L70,L71 
DATA  >1/3.141592655/ 

R  H  1  =  1 


I 

I - [  INPUTS  3 

I  t  1)  T  =  T W 0  WAY  TRAVEL  TIME  3 

I  [THICKNESS  OF  SEDIMENT  3 

I 

I - [  2)  NNFILE  =  ARBITRARY  FILE  3 

I  [NUMBER  3 

I 


I - [  3)  I  OPT  =  INPUT  OPTION,  0  EO*  3 

I  [GEO-PAR A*ETERS ,  3 

I 

I 


/ . / 

/  WRITE  C  1,701 )  / 

/ . / 


I 

I - [701  E  ORP  ®  T ( '  INPUT  OPTIONS  APE:  7 

I  [0  =  INPUT  GEO-PARAMETERS', /,  *  '3 

I  [1  =  INPUT  POTTOP  LOSS',/,  *  '  2  3 

I  [  END  PROGRAM',//,  *  '  ENTER  3 

I  [STATION  NUPPER  AND  INPUT  0RTI0N:3 

I  [)  3 

I 
I 


/ . / 

/  PEADd,*)  NNFILE, I  OPT  / 

/ . / 

I 

I 
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’  *  *  *  *  *  *  *  4 


********** 


************* 


<  I FCI0PT.EQ.2)  GO  TO  999  > 

\ . / 


/ 


I  FALSE 

I 

I 


***************************************** 


/ 


/  URITE(1,702)  / 

/ . / 


TRUE 


I 

I - [702  FORMAT (  '  ENTER  THE  TWO-WAY 

I  [TRAVEL  TIME  THICKNESS  IN  10THS 
I  [OF  SECON  *  D  S  :  ' ) 


I 

I 

/ . . . . ' . / 

/  RE  At)  Cl  »*)  T  / 

/ . / 

I 


I - [  BRANCH  TO  EXIT 

I 

I 


D 

D 

D 


D 


T  =  T  /  2  0  . 


I 

I - [  INPUT  experimental  geometery  T 


I  [PARAMETERS  D 

I 

I 

/ . / 

/  WRITE (1,705)  / 

/ . / 

I 


I [705  PORMATC'  ENTER  THE  WATER  1 

I  [DEPTH:')  ? 

I 

I 


/ . . / 

/  READ (1  ,*)  ZR  / 

/ . / 

I 

I 

!  **************  *************************** j 

/  WR ITE (  1  *  706 )  / 

/ . / 


I 

I [706  FORMAT (  '  ENTER  SCUPCE  DEPTH  3 

I  [AND  RECEIVER  DEPTH:')  D 

I 

I 


/ . . . "/ 

/  READC1 ,*)  ZS, ZR  / 

/ . / 

I 

I 


0 

7 

7 

Z 

Z 

z 

7 

Z 

z 

z 

z 

z 

7 

7 

7 

Z 

7 

7 

7 

7 

7 

7 

Z 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

•» 

4k 

2 

7 

7 

2 

7 

7 

2 

7 

7 

7 

7 

7 

7 
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/ . . . — '/ 

/  WPI TE ( 1  *707)  / 

/ . . 


I 

I - [707  formats  enter  sound  speed  3 

I  [AT  SOURCE  OR  RECEIVER  -  3 

I  [WHICHEVER  IS  SH  *ALLOWER:')  3 

I 
I 


/ . . / 

/  READd.O  CS  / 

/ . / 

I 

I 

/  . . . . . / 

/  WRITEC1.708)  / 

/ . / 


7 

7 

7 

Z 

Z 

7 

7 

Z 

Z 

z 

7 

7 

7 

7 

Z 

7 

z 

7 

Ar 


/ 

/ 

/.  . 


/ 

< 

\ 


/ 

/ 


/ 

/ 

/  .  . 


i 

I - C  70  8  FORMAT^  ENTER  SOUND  SPEED  3 

I  [AT  THE  BOTTOM  OF  THE  WATERS)  3 

I 

I 

- - - - 

READ  < 1 f*>  Cl  / 

. . . / 


I 

I - [  BRANCH  TO  INPUT  OPTION  3 

I 

I 

-  - . .  . . \  TPUF 

I  F  (  IOPT  .EG  .D  GO  TO  5  0  C  > - 

. / 

I  FALSE 
I 

I - [  INPUT  INVERSION  PARAMETERS  0 

I 

I 

. ' . . . . / 

WRITE  (  1  ,709)  / 

. / 

I 

I - [7C9  FORMAT ( *  ATTENUATION  3 

I  [CALCULATION',//,  ♦  '  FOR  A  POINT  1 

I  [BELOW  THE  CAUSTIC  VHERF  PL  CURVED 

I  [IS  NEARLY  LINEAR,'/  *,'  ENTER  A  1 

I  [LOSS,  FREQUENCY  AND  ANGLE:  (  ANGO 

I  [APPROX  20,  FREG  =  APPRO  *X  200) ' D 

I 
I 

. . . . / 

R  E  A  D ( 1  ,  *  )  DPLOSS  ,  FREC f THETA  / 

. / 

I 

I 


-0 

X 

X 

Y 
X 
X 
X 
X 
X 
X 

Y 

Y 

Y 
X 
y 
x 
x 

Y 

Y 

Y 

Y 
X 

Y 
X 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

■? 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 


/ . ' . . . / 

/  WRITE(1.710)  / 

/ . . . / 


I 

I 

I 

I 

I 

I 

I 

I 


/  PEA0C1.O  ALOSSO  ,  ALOSSM , FM  / 

/ . / 

I 

I 

/ . . . ' . . . ' . / 

/  WRITE ( 1  *711 )  / 

/ . / 


I 

I - C^ll  FORM  AT  (  '  ENTER  ANGLE  OF  3 

I  [TOTAL  INTERNAL  REFLECTION,',/,  *3 
I  t'  ANGLE  OF  CAUSTIC,',/,  *  '  3 

I  [RATIO  OF  SPEED  IN  SEDIMENT  TO  3 

I  [SPEED  AT  BOTTOM  OF  WATER:')  3 

I 
I 


/ . . . . . . "/ 

/  RE AD(1 ,*)TIR,THC, RATIO  / 

/ . / 

I 

I 

/ . . . . . . \ 

<  IF(THC.EQ.O.C)  then  > -  - 0 

\ . / 

92  9  9  I 
I 
I 

/ . . . . . / 

/  WRITE  Cl,  712)  / 

/ . / 


I 

I - [712  FORMAT ('  ENTER  GRADIENT  AT  3 

I  [TOP  OF  SEDIMENT:')  3 

I 

I 


/ . . . . . / 

/  RE«D(  1,0  G  / 

/ . / 

I 

I 


: 


—  [710  FORM  AT  (  '  ENT  FR  LOW  FREQ  LOS  S  3 
[LEVEL  AT  HIGH  ANGLES,',/,  *  '  3 

[the  minimum  loss  level  at  high  3 
[ANGLES',/,  *  '  AND  THE  FREQUENCY  3 
[OF  THE  MINIMUM  LOSS:')  3 


ENDI  F 


. / 

X 

7 

/  WRITE  < 1 *713) 

/ 

X 

7 

x 

7 

I 

X 

7 

I - ["H3  FORM  AT ( ' 

ENTER  CANONICAL 

3 

X 

7 

I  [CURVE  TYPE  AND  ATTENUATION 

3 

X 

7 

I  Cpro  FILE  GRADIE  *NT:'> 

3 

X 

7 

I 

X 

7 

I 

X 

7 

1*...,., . .. ...  """" 

X 

7 

/  RE  AO ( 1 «  * )  BETA, GAM 

/ 

X 

7 

x 

7 

I 

X 

7 

I 

X 

7 

/ . . . " . . . . . 

- - 7 

X 

7 

/  WRITE (1 ,714 ) 

/ 

X 

7 

/ . . . . 

.  .  .  / 

X 

7 

I 

X 

7 

I - [714  FORMATE' 

enter  basement 

3 

X 

7 

I  [REFLECTION  COEFFICIENT  IN 

3 

X 

7 

I  [PRESSURE  UNITS:'  O 

3 

X 

7 

I 

X 

7 

I 

X 

7 

/' . 

X 

7 

/  R  E  A  0  ( 1  ,  *  )  R  E  F 

/ 

X 

2 

/. . 

III/ 

X 

7 

I 

X 

7 

I - [  OPTIONS  AND 

ERROR  CHECKING 

3 

X 

7 

I 

X 

7 

1 

X 

7 

. . 

""\  FALSE 

X 

7 

<  IFCTIR.GT.C.  .OR.  RATIO  .LE.  0.) 

> - 

-I 

X 

7 

I 

X 

7 

I  TRUE 

I 

X 

7 

I 

I 

X 

7 

I 

X 

7 

:  RATlO=1./COS(TIfi*PI/130.) 

: 

I 

X 

7 

• 

• 

I 

X 

7 

I 

I 

X 

7 

A  ✓  — _ _ __ 

y 

7 

u  * 

I 

x_ 

7 

X 

7 

:  C?=C1*RAT10 

X 

7 

:  C3  =  C2 

X 

7 

X 

2 

I 

X 

7 

I - [  I F  GRADIENT 

G  IS  NOT  SPECIFIED 

3 

X 

7 

I  [CALCULATE  IT 

FRO**  THE 

3 

X 

7 

l 

I  [  CAUSTIC  ANGLE,  TH  C 

3 

X 

7 

I 

X 

7 

I 

X 

7 
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IF(G.*E.O.)  bO  TO  300 


\  TRUE 

> - 

/ 


I  FALSE 

I 

I 


I F ( THC • L  £. 0 ♦ ) 

•  •••••••••••• 

I  TRUE 


\  FALSE 
> - 

/ 


: 

THC  *30 

• 

1 

I 

I 

I 

I 

I 

•  : 

CALL  GCOMP 

•  # 

•  • 

I 

r\  s  — 

c  3003 

U  v 

I 

: 

CONTINUE 

: 

I 

I 

/  "" 

/ 

/ . 

we  I TE ( o, 1 40) 

. / 

/ 

. / 

I 

i---[  no 

i 

i 

FOR*AT(//,n  INFUT  DATAn, 

/)  3 

/ . . 

/ 

/ . 

Wf ITE (6 ,21 0)NNFILE 

. . .  f 

t  T  / 

. .  j 

I 

1 ---[  210 

I 

1 

FORMAT ( IX  1 14  v /•»  T  =o,F° 

.2)3 

/ . 

/ 

/ . 

WRITE (6t220)2btZS*ZRt 

- - - ! 

CS.Cl  / 

I 

I - C  220  FOR^AT(d  ZP=d,F8.2,p  7S  =  d,  2 

I  [Fg.2fn  Zp  =  DfF8.2  tD  CS  =  DtF8.2,D  ? 

I  CCl=DtF?.?  1)  ] 

I 

I 


WpITE(6,U1)DBL0SStFREGtTHETA 


•0 

U 

V 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

II 

u 

►0 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


7 

7 

z 

z 

7 

z 

7 

7 

7 

7 

7 

7 

7 

7 

7 

Z 

7 

Z 

z 

z 

7 

7 

7 

7 

7 

f 

Z 

7 

Z 

7 

7 

7 

7 

7 

7 
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r^t  f-j 


I 

I C  Hi  FOR*AT(b  DBlOSS=o,FP.2fD  3 

I  CFREC  =ntF8.2,n  THETA  *n,F8.2>  3 


1 

1 

/  . . . . . . . / 

/  WFlTE(6,U2)ALCSS0,ALOSSKfFH  / 

/ . . 

I 


I - [  142  FORMATCb  ALOSSO«n, F8.2,d  3 

I  [ALOSSM=ntF8.2,B  FM  *b.F8.2)  3 

I 

I 


/ . / 

/  WRITE(6«14?)TlRvTHCt RATIO  / 

/ . / 


I 

I - C  143  FORMAT  (b  tip  roje^.n  THC  3 

I  CntT8.2,n  RATIO  *BfF8.4>  3 

I 

i  —  C  change  degrees  to  radians  ] 

i 

i 


THETA  =  THETA*PI /I  80. 
THC  =THC*PI /I  80. 

T  I  R  =T  I R  *  P I / 1 £0. 


I 

I - C  CALL  ROUTINE  WHICH  CALCULATES  0 

i  [the  regaining  l 

I  [  GEOPHYSICAL  parameters  from  theo 

I  [INVERSION  INPUTS  ] 

I 
I 


;  ; 

CALL  GEC-PY 

•  • 

I 

I 

/  " 

: 

GO  TO  499  : - 

•-  0 

\ . . , 

u 

u 

I 

u 

I - [  OPTIONAL  PRANCH  TO  INPUT  GEO- 

3 

l» 

I  [0ARA",ETEPS  DIRECTLY 

0 

u 

I 

r  /  _  _ 

V 

V*- - - 

[  5  00  D 

I 

V 

V 

• 

• 

CONTINLE  : 

V 

• 

• 

u 

I 

u 

I 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

r 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

x 

x 

X 

X 

X 

X 

X 

X 

x 

X 

X 

X 

>0 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

1 

7 

7 

7 

7 

7 


J  20 


I 

U 

7 

'  . . .  ' . . . / 

U 

7 

/ 

URITEC1.715)  / 

U 

7 

/ . 

li 

7 

I 

V 

U 

L 

7 

I - C71 5  FORMAT ( *  ENTER 

RATIO  OF 

3 

U 

7 

I  [SPEED  IN  SEDIMENT 

TO  SPEED  AT 

3 

u 

7 

I  [BOTTOM  OF  VA  *TE  R : 

') 

3 

u 

7 

I 

u 

7 

I 

V 

7 

J +++++++ 

^  ^  /  /  /  /  /  /  /  /  y 

u 

l 

/ 

READd  ,  *  )  RATIO  / 

u 

*7 

* 

/,,, . 

1  f 

7 

I 

V 

L 

7 

I 

u  * 

7 

+  ++  +  +  +  +++++++  +  ++  +  +  +  +  ++  +  ++  +  +  +  +  ++  +  +  +  j 

u 

7 

/ 

WRITE(1*716)  / 

u 

7 

/ . 

. ../ 

(J 

7 

I 

u 

7 

I --- [71 6  FORM  AT ('  ENTER 

GRADIENT  AT 

3 

11 

7 

I  [TOP  OF  SEDIMENT  AND  CANONICAL 

3 

u 

7 

I  [CURVE  TYP  *E:') 

3 

u 

7 

I 

u 

7 

I 

u 

7 

/ . 

u 

7 

/ 

READd  ,*  )  6, BETA  / 

u 

7 

/ . 

. / 

D 

7 

I 

u 

7 

I 

t) 

7 

'  '  '  - - / 

l’ 

7 

/ 

ta^ITE (1 ♦ 7 1 7 )  / 

IJ 

7 

/ . 

. . . / 

LI 

7 

i 

u 

7 

I - [717  FORMATd  ENTER 

THIN  LAYER 

3 

u 

7 

I  [DENSITY,  SEDIMENT 

DENSITY  AND 

3 

u 

7 

I  [THIN  LAYFR  *  THICKNESS  IN  M  E  T  E  R  S  3 

i» 

7 

I  [') 

3 

u 

7 

I 

11 

7 

I 

u 

7 

/ . 

. .  '/ 

V 

7 

/ 

READd  » *  )  RH2,rh3,D  / 

u 

7 

/ . 

. / 

7 

I 

u 

7 

I 

u 

7 

f"'"" 

. . . ! 

u 

7 

/ 

WPITE (  1*716)  / 

U 

7 

L. 

/ . 

. / 

u 

7 

i 

u 

7 

I - [718  FORMAT ( *  ENTER 

SURFACE 

3 

J 

7 

I  [ATTENUATION  IN  DB/M/KHZ:') 

3 

u 

7 

I 

u 

7 

I 

u 

7 

/ — "" 

. . . . . . / 

l! 

7 

/ 

READd,*)  ALPHO  / 

u 

7 

/ . 

u 

7 
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I 

U 

7 

I 

U 

7 

. . . / 

u 

7 

/ 

WRITEC1 

,719)  / 

u 

7 

/ . . . . . . 

II 

7 

I 

u 

£ 

7 

I-- 

-[719  FORM  AT  (  '  ENTER  ATTENUATION 

3 

u 

7 

I 

[PROFILE  GRADIENT  AND  BASEMENT 

3 

u 

7 

I 

CREFLECTIO  *N  COEFFICIENT  IN 

3 

u 

7 

I 

[PRESSURE :*) 

3 

u 

7 

I 

u 

7 

I 

u 

7 

/ . 

" .  . .  l 

u 

7 

/ 

PE  A D ( 1 ,*) 

GAM, REF  / 

u 

7 

/ . 

u 

7 

I 

u 

L 

7 

I 

u 

7 

11 

*? 

C2=C1*RATI0 

u 

4. 

7 

C3=C2 

u 

7 

u 

7 

I 

u 

L 

7 

I-- 

-[  CALL  A  SUPRUUTINE  WHICH 

3 

u 

7 

I 

[CALCULATES  AND  PRINTS  VARIOUS 

3 

u 

7 

I 

[  CRITICAL  ANGLES  CHARACTERISTIC 

3 

t* 

7 

I 

[OF  THE  MEASUREMENT 

3 

u 

7 

I 

r\  s  «, 

u 

7 

U  * 

z 

[  4003 

1 

7 

7 

;  : 

CALLCANGLE  :: 

7 

• 

•  • 

7 

I 

7 

I - 

-[  PRINT  GEO-PAR AMETERS 

3 

7 

I 

I 

/ . / 

/  U'RITE(6,145)  / 

/ . / 


I 

I - C  145  FOR"AT(//,n  G E 0 -P A R AM F T E R S  3 

I  [  .//)  3 

I 

I 


/ . / 

/  VRITE(6,146)RATIO,Rh2,PH3,D,G,PETA,  / 

/  ALPHCtGAK,B7,C’EF  / 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 


. /  7 

I  7 

I - [  146  FOP*AT(n  R  A  T  =  c  ,  F  8 . 4  ,  /  t  o  Rh23  7 

I  [o,F8.2»n  RH3*b,F?.?,/,  1“  D  =n ,  3  7 

I  [F6.2,/,o  G  =°,F£.2,d  BET=b,F8.2,3  7 

I  C»  2°  ALO  =  b  ,  F  <?  .3  ,  AN  =  n  ,  Fl  1  .  5  ♦  /  ,  n3  7 

I  C  9  2  =ctF8.2,n  REF=n,FS.2)  3  7 

I  7 

I - [  CALL  ROUTINES  FOR  FREQUENCY-  3  7 

i  [independent  intermediate  3  7 

I  7 

I  7 
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CALL  REFL 


I 

I 


CALL  REFR 


I 

I 


CALL  BASE 


I 

0< - 

C  999}  I 

/ . . \ 

:  RE  T  URN  : 

\ . / 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

0 
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\  SUBROUTINE  GEOPY  / 

\ . / 

I 

I 


COMPON  Z8,ZR.ZS.C1.C2*C3,RH1,RH2,RH3,G*D, 
ALPHO.CS.BZ.REF 

COMPON/ DAT A /L ZERO, FZERO.L70 ,F70,TH70 ,171 , 
F71, DP  LOSS, FREQ  ,T  HET  A  ,  ALOSS  0  ,  ALOS  SM  ,  FP. , 
TIR  ,THC  ,THRC  ,THTB  .BETA 

COMMON  /ALPH/GAM 
REAL  L  ZE  RO  »  L70 ♦ L71 
DATA  PI /3. 141592655/ 


I 

I - COPPUTE  DENSITIES  RH2.RH3  3 

I 

I 


Z1  =C1  *  Rnl 
D=Cl/FM/4 

R  =  E  X  P (-ALOG (10. ) *AL0SS0/20) 
Z3=(T*P)/(1-R)*Z1 
R=EYP(-ALOG (10. )*ALOSSM/20) 
Z2=  ( 1  ♦R)  /<1-R>*  Z1*Z3 
Z2=SQRT ( 12) 
SPRAT=C3/C1 
RH2  =  7  2 / C2 
PH  3=  Z  3/ C3 


I 

I - COPPUTE  DISSIPATION  CONSTANT  3 

I  CALPHO  3 

I 

I 
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THETA=ACOS( C2/C1 *  C  0  S (THETA) ) 
eS=$IN(THETA) 
eC*COS( THETA) 
V32C(Z3-Z2)/(Z3+Z2) 
V21*(Z2-Z1)/(Z2+Z1) 

AUX=4*D*PI *FREQ/C1*SIN(THETA) 

V  =  C  A  B  S ( (V21 ♦V32*CPPLX ( COS (AUX) «  SIN(AUX) ) ) 
/(1+V?1*V32*CPPLX(COS(AUX)tSIN(AUX)))) 
T*CABS<1/<1+V21*V32*CWPLX(C0S(AUX),SIN( 
AUX)))>**2*16*ZT*Z3*<Z2/CZ1*Z2)/CZ2+Z3>>* 
*2 

AUX=10.**  (-DBLOSS/10. )-V**2 
AUX=AUX/T 

ALS=- ALOG1 C ( AUX ) *  1 0 
Fr=FREQ/1000 
S 1  G  =  G*  (1  +BE  TA  )  *  C  2 
A*U*< C2*C1+BETA) )**2 
ANU  =  BET A  *C  2 

ALPH0=ALS/FX-2*C2/G/(1+BETA)/PC*(GAP*C2** 
2/2/SIG/BC**3*(PS~BS**3/3)+(BETA*GAP*C2** 
2+2*GAP*ANU*C2)/4/SIG/bC**2*(THETA+BS*BC> 
♦  ( G  A  M  * (ANU**2~APU) +2*9ETA*GAP*ANU*C2 ) /2 / 
SI6/BC*BS+BETA*GAP*(ANU**2-AMU)/2/SIG* 
TMET A ) 

ALPH0  =  ALP  HO /(eS/BCA-P  ETA*  THETA)  *G*(1+-BETA) 
♦9C/2/C2 


I 

I 

/ . \ 

:  PET  URN' 

\ . / 
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\  SUBROUTINE  CANGLE  / 

\ . / 

I 

I 


COMMON  ZB,ZR,7S*C1.C2tC3,RHl,RH2,RH3,G,D, 
ALPHO»CS»EZ»REF 

COMMON  /ALPH/GAM 
COMMON  /T/T 


I 

I 


COMMON  /DATA/  L 7E RO , F Z ERO ,L 70 , F 70 ♦ T H7C , 
L71 , F71 »  DPLOS  S  t PREG*THETA*ALOSSO»ALOSSM» 
FM  ,T  IP  ,TH  C  tTHRC  ,THlB  ,B  ETA 


I 

I 


REAL  LZEp0,l70fL71 


I 

I 


DATA  PI  73.14159264/ 


I 

I 


D  I  =2  .*ZP-ZP-ZS 


I 

I - C  CALCULATE  sediment  thickness 

I 

I 


R  =P  F  T  A 
CZP=2000 

A  =  C2*(G*(1.-»B)*T+1.  +  P*ALOG<C2>) 


I 

0< 

C  603  I 


CZ  =  CZr+(A-C7P”6*C2*AL0G(C7P) )/(1  ,  +  P*C2/ 
C  ZP) 


I 

I 

/ . . . . . . \  TRUE 

<  I F ( ABS ( CZ-C 7P) . LT.  .01 )  GO  TO  ^0  > - 

\ . / 

I  FALSE 

I 

I 


3 


0 

X 

X 

X 

X 


0 

7 

7 

7 

Z 

7 

7 

7 

7 

7 

7 

7 

7 

7 
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CZP=CZ 


I 

I 

:  GO  TO  60  : 

\ . / 


0< 

C  7C3  I 


BZ=((CZ+B*C2)**2-C2**2*(1+e)**2)/(2*G*(1* 
B )  *  C2  ) 


I 

I - C  COMPUTE  REFINED  CRITICAL  ANGLE  3 

I 

I 


X 

X 

X 

-0 


; 

THRC  =  P  1/4 

: 

I 

I 

*  z 

CALL  TOTCI 

•  * 

•  .  •  ***  +  *  +  + 

•  • 

I 

I - C  COMPUTE  CPITICAL  RAY  DEPTH  CRD  3 

I 

1 


CAUX=C1 /C05(THRC) 

CRD=(CCAUX+PETA*C3)**2~(C3*(1+BEtA))**2) 
/2/G/(1+BETA) / C  ? 


I 

I [  CALCULATE  VELOCITY,  VELOCITY  3 

I  [GRADIENT,  AND  ATTENUATION  3 

I  [  CONSTANT  AT  DEPTH  OF  500  ME  T  E  R  S  3 
I 
I 


BSPD=$GRT(2*G*(1*9ETA)*C3*500*C3**2*(1* 

BETA)**2)-BETA*C3 

PGRAD=G*(1+BETA)/(0SDD/C?+BETA) 

ALBOT=ALPH0+GA»*500 


I 

I - [  COMPUTE  TWO  WAY  TRAVEL  TI**E  IN  3 

I  [SECONDS  3 

I 

1 
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*•••••••••• 

TWTT=2  *T 
SPRAT=C2/C1 


I 

I - C  COMPUTE  BASEMENT  CRITICAL  ] 

I  CANGLES  THl B » TH3B  ] 

I 

I 


C B= SORT (2*G*  <1  +  BE TA>*C3*B7+C3**2*(1  ♦BETA) 
**2>-BETA*C3 

THl P=ACOS (Cl /CB) 

TH3B=AC0S ( C3/CB) 

F=DI/TAN(THlP)42*C3*TAN(TH3B)/G*<BETA/( 
1+BFTA)) M 1*(1+TH3B/SIN(TH3B)/C0S (TH3B) ) / 
2/BETA) 


I 

I C  COMPUTE  SMALLEST  APPARENT  ANGLE] 

I  CEOR  BASEMENT  REFLECTION  ? 

I 

I 


THP=ATAN( 01/ P) 


I 

I - C  CHANGE  RADIANS  TO  DEGREES  ] 

I 

I 


THP=THR*1 SO / P I 
A=T I R  *  1 8C/P I 
A 1=THC  *1 tC/PI 
A2=THRC*1 80/PI 
A3=THlP*1 80/ PI 


I 

1 - [  PRINT  OUT  CRITICAL  ANGLES  ] 

I 

I 

/ . . . / 

/  WRITE (6 ,1 )  / 

/ . / 


I 

1 - 1  1  FOP*AT(1X»//»D  critical 

1  [ANGLESo,/)  T 

I 

I 


/ . . . . . '/ 

/  WRITEC6, IDA, SPRAT, A1,A2,THB, A3, TVTT,/ 

/  BZ  ,  CRD  ,PGRAD , PSPD ,ALB0T ,GAM  / 

/ . / 
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I 

I - C  11  FORMAT  (n  TIR  *n,F8.2,n  RAT  =  3 

I  C,F8.4,n  THC  *n,F8.2?n  THRC  »n,F83 

I  C  2 »  1/,n  THE  *n,F8.2,*»  THRB  =n?F83 

I  C  2 ,  2/,n  TWTT  *o,F8.2,n  BZ  en,F8.3 

I  C2,n  CRD  =b,F8.2,/,  3«»  G500  =  «*,f83 

I  C2,n  C500  *n,F8.2,n  *500  =n,F9.3,3 

I  C  6AM  *n,F11.5>  3 

I 

C  9993  I 

/ . . . \ 

:  PETURN  : 

\.... . / 
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J 


\  SUBROUTINE  TOTCI  / 

V . / 

I 

I 


COUPON  ZB,ZR,ZS»C1,C2,C3,RHl,RH2tftM3,G,D, 
ALPHO,CS,BZ.REf 

COUPON  n> AT* /LZEPO,FZEftG,L70,F70,TM70,t71  , 
F 71, DB LOSS, FREQ  ,TMETA , ALOSSO , ALOSSP , FP, 
Tift,  THC.THRC.THlB, BETA 

REAL  L  ZERO  ,  L70  ,L71 


1 

I - C  CHECK  THAT  C3  LESS  THAN  C? 

I 

I 

/  *************  ""  ""\  TRUE 

<  IF (C3  .LE.C2 ) GO  TO  1C  > - 

\ . / 

I  FALSE 

I 

I 

/ . . . . . . .  t 

/  WRITE  (1  ,1 1 )  / 

/ . / 


I 

I - [  11  FORMAT (n  C3  MUST  BE  LESS 

i  [than  or  equal  to  c2«j> 

I 

I 


/ . \ 

:  stop 

\ . / 


0<  — 

[1C3  I 


PI=3.U1592654 


3 


I 


I C  COMPUTE  TIR  =  ANGLE  OF  TOTAL 

I  [INTERNAL  REFLECTION 

I 

I 


3 

3 


TIP  =  0 


I 

I 


. . . . "  \  TRUE 

<  IF(C2.LE.C1)GG  TO  100  > - 

\ . / 

I  FALSE 

I 

I 


• 

» 

TIR=AC0S(C1/C2) 

• 

• 

• 

I 

r / •••••••••••. 

C  100  3 

1  v 

1 

: 

CONTINUE 

m 

• 

• 

I 

I - C  COMPUTE 

I  CCAUSTIC) 

THC=CRITICAL  ANGLE  ( 

3 

3 

I 

I 


0 

z 

7 

Z 

z 

z 

z 

z 

z 

0 


PI =3 • 1 A 1 592654 
P I  2  -  P I  /2 
D  EL  T  =  P I /ISO 
M =2*7B-ZR-ZS 
D2=2*C3/G 

D2  =  D  2*BETA /(I ♦BETA) 

FAC  =  C3  /Cl 

TH3=AC0S(FAC*CC$(THRC)) 

R=Dl /TAN (THRC) ♦e2*TAN(TH3)*(1M1*TH3/SIN( 
TH3)/COS(TH3))/?/PETA) 


I 

I - C  FIMFI  Zc  EXPRESSION  FOR  RANGE  1 

I  CPY  CONJUGATE  DIRECTIONS  3 

I 

I 


: 

DO  200  1=1 ,1000 

: 

I 

I 

2 

AUX=THRC+  DELT 

• . 1 

A  I 

A  I 


A  / . V  TRUE 

A  <  IE(AUX.GE,PI2)G0  TO  5CC  > - 

A  \ . / 

A  I  FALSE 

A  I 

A  I 

A  / . \  TRUE 

A  <  I  E  (  AUX . LE . TIf ) GO  TO  5C0  > - 

A  \ . / 


A 

A 

A 


C 

Z 

Z 

Z 

7 

7 

V 

7 

Z 

7 

Z 


I  FALSE 

I 

I 


131 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


:  T  HR  C=  A  UX 

:  TH3=AC0S( FAC* COS (THRC) ) 

:  RAUX=Dl/TAN(THRe)+D2*TAN(TH3)*(1*(1+TH3/ 
:  SIN(TH3)/COS(TH3))/2/BETA) 


<  IF(RAUX.GT.R)bO  TO  500  > 

\ . . / 

I  FALSE 


I 


TRUE 


R=RAUX 


I 

1 

/ . . 

:  GO  TO  200  : 

\ . / 


0< 

C  500  3  I 


R=P AUX 

D  ELT=-DELT/3 


I 

I 

. . . . . . TRUE 

<  I  F  (ABS CDELT  3.LE  .1  .E-6)  GO  TO  999  > - 

\ . . . / 

I  FALSE 
I 

0< - 


C  2003  I 


CONTINUE 


0 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

V 


X 

X 

X 

0 


I 

I 

/ . . .  . / 

/  WR ITE (1 f 1 )  / 

/ . / 

I 


I C  1  FORMAT ( o  CRITICAL  ANGLE  3 

I  [FAILED  TO  CONVERGE^)  3 

I 

I 

/ . . . . \ 

:  STOP  : 

\ . / 


7 

7 

7 

7 

7 

7 

7 

7 

7 

V 

7 

7 

Z 

7 

7 

y 

L. 

7 

7 

7 

7 

7 

7 

7 

o 


0 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 
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z 

I - [  BOTTOM  ANGLE  OF  REFRACTED  RAY  D  Z 

I  CAT  CRITICAL  ANGLE  J  Z 

I  Z 

0  < - 0 

C  999D  1 


THC=ATAN( D1 /R) 

CB=SQRT(2*G*(1*BETA)*C3*BZ+(C3*(1*BETA))* 

*2)-BETA*C3 

TH1 P=  A C 0 S C  C 1 / C B ) 


I 

I 


RETURN 
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•  ••••••••••••••••  •  ••••••••••••••••  ••••••••••• 

\  SUBROUTINE  REEL  / 

\ . / 

I 

I 


COUPON  ZB,ZR,ZS,C1,C2,C3,RH1,RH2,RH3,G,D, 
alpho,cs,pz,ref 

COMPON/ DATA /LZERO,FZERO,L70, F70,TH70,L71 , 
F71 ,DBLOSS, FREQ  *THETA, ALOSSO, ALOSSM, FK, 
TIR ,THC ,THRC ,THlB ,BETA 

REAL  LZER0,L70,L71 

COMPON  /R/NSAP,RV21,RV32,RSIN2,RDTS,RDTR, 
R  T  ,  R  G I 

DIMENSION  RV21 (90) ,PV32(90) ,RSIN2(9P)  , 
RDTS (90) ,RDTR (90) ,RT ( 9  u) ,RGI (90) 


I 

I  [********************************3 

i  [***  : 

I  C  RETURNS  IN  COMPON  T 

I  C  NSAP=NSAP=SPALLEST  ANGLE(DEG)  T 

I  CTO  REACH  SURFACE  DUE  TO  SOUND  3 

I  [PROFILE  T 

I  C  RV21 ,RV32*REFLECTI0N  COEFFICIENT 

I  CS  FOR  RESPECTIVE  INTERFACE  T 

I  C  RSIN2=STN  OF  GRAZING  ANGLE  IN  T 

I  CHIGH  SPEED  LAYER  T 

I  C  RDTS ,RDTR=RFFLECTION  PATH  DELTA} 

I  CT  FOR  SOURCE, RECEIVER  RESP.  T 

I  C  RT=TRAVFL  TIME  FOP  REFLECTION  T 

I  [PATH  T 

I  C  RGI=GEO^ETRIC  INTENSITY  OF  T 

I  [REFLECTION  PATH  T 

I  £********************************} 

I  [***  } 

I  [  COMPUTE  NSAP  T 

I 
I 


•  •  •••••••••••• 

PI =3. 1415 92654 
N  S  A  P  =  1 


I 

I 

. . . . TRUE 

<  IF (Cl .GE .CS)GO  TO  100  > - 

\ . / 

I  FALSE 

I 

I 


0 

7 

Z 

7 

7 
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SAP=ACOS (Cl / C  S)*180/PI 
NS AP=SAP+1 


I 

I 

t  . . . / 

/  URITE(I.I)  SAP  / 

/ . / 

I 


I - C  1  FORMAT ( **  SMALLEST  ANGLE  TO 

I  [REACH  SURFACE***,  F8.2) 

I 

0< - 

C  1  CO  3  I 


] 

] 


7 

7 

7 

7 

7 

7 

1 

7 

7 

7 

7 

7 

7 

0 


CONTINUE 


I 

I  - 
I 
I 
I 


C  COMPUTE  REFLECTION  COEFFICIENTS] 
[FOR  ANGLES  GREATER  THAN  TIR.SAP  ] 


NTIR=TIR*1P0/PI+1 
N  =  "A  IC ( NT  I R  ,  N  S  A  P ) 
71  =RH 1 *C  1 
Z2=RH2*C2 
Z3  =  RH3*C3 
DELT  =  PI /  180 
T H  =  (NSAP-1 ) * D  E L T 
FAC  2  =  C  2 /C 1 
FAC3  =  C3  /Cl 
FACS=CS/C1 
AR  =  2*Z«?/C  S 
AS=2*ZS/CS 
PI =2*  ZP-ZR-ZS 


I 

I 


00  200  I=NSAP,90 


I 

I 


7H=TH+DELT 
C  TH  =  COS ( 7 H ) 


I 

I - [  SURFACE  SCATTERING  TRAVEL  TIME 

I  [DIFFERENCES 

I 

I 


] 

5 
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A 

A 

A 

A 


S  =  S  0  R  T ( 1“ < FACS*CTH) **?) 
R  DTS (  I  ) *A$*S 
RDTR( I)=AR*S 


A 

A 

A 

i 

i 

A 

/'■ 

. \  TRUE 

A 

< 

I F ( I • LT • N ) GO  TO  200 

> - 

A 

\  .  . 

. / 

2 

A 

I  FALSE 

Z 

A 

I 

Z 

A 

I - [  SIN? 

Z 

A 

I 

Z 

A 

I 

Z 

A 

2 

A 

• 

PSIN2(I)=SORT(1-( FAC2*  CTH) *  *  2 ) 

• 

» 

Z 

A 

• 

• 

2 

A 

I 

■ 

Z 

A 

I - C  REFLECTION 

COEFFICIENTS 

Z 

A 

I 

Z 

A 

I 

Z 

A 

2 

A 

A  1 =Z1 /SIN (TH) 

Z 

A 

A2  =Z?/RS  IN?  ( I  ) 

Z 

A 

A3=Z3/SQRT(1-(FAC3*CTH)**2) 

Z 

A 

RV3?(I)=(A3-A2)/(#3+A?) 

Z 

A 

R V2 1 (I)*(A2-A1 )/ ( A2  +  A1 ) 

Z 

A 

2 

A 

I 

Z 

A 

I - C  TRAVEL  T I P  E 

l 

Z 

A 

I 

Z 

A 

I 

Z 

A 

2 

A 

F=Dl/TAN(TH) 

Z 

A 

R  T ( I)=R/(CTH*C1) 

Z 

A 

2 

A 

I 

Z 

A 

I - C  GEOMETRIC 

INTEN'SITY(AMPLITUOE) 

Z 

A 

I 

z 

A 

I 

z 

A 

7 

A 

• 

R  r,i  (  I  )  =  CTH/P 

• 

Z 

A 

• 

m 

7 

A 

I 

7 

A  S  _ 

A 

A 

c 

2or  i  i 

A 

CONTINUE 

m 

• 

• 

I 

I 

/ . \ 

:  RETURN  : 

\ . / 
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\  SUBROUTINE  RE  FR  / 

\ . / 

I 

I 


COMMON  ZB  ,ZR  ,ZS  ,C 1 ,C2,C3»RHl ,  RH2 » »H3 , G , D ♦ 
ALPHO,CS,BZ,REF 

COMMON/DATA/LZE  *0, F  ZER0,L70 , F 70 » T H 70  ,  L 7 1  , 
F  71 t  DBLOS  S  »  FREQ *T  H  E  T  A ,ALOSSO»ALOSSMt  FM  * 
TIR,THC,THRC,TH1B,BETA 

COMMON  /ALPH/GAM 
REAL  LZERO, L70.L71 

COUPON  /B/BV21,BV32tBSlN2fBTRF,BGIfBLS, 
BT,PDTS  *B  P  T  R 

COMMON  /S/SV21»SV32«SSlN2»$TRF,SGI»SLSt 
ST»SDTStSDTR 

DIMENSION  BV21(90) »BV3^(9C) »6 S  I N  2 (9D  « 
BTRF(OC)tPGI(90),BLS(9u),BT(90),BPTS(90), 
B  DTE (90) 

DIMENSION  SV21(90),SV3Z(90),SSIN2(90), 
STRF(9C),SGI(9O),SLS(9C>,ST(90)tSDTS(9O>, 
S  PT  P ( 90  ) 


I 

I— f«*******«***********************n 
1  [****************#***************3 


I  [***  3 
I  [  FOP  EACH  SPECULAR  BOTTO"  ANGLE  3 
I  [BETWEEN  THE  ANGLE  OF  TOTAL  3 
I  [INTERNAL  3 
I  [  REFLECTION  AND  THE  CRITICAL  3 
I  [ANGLE  THIS  SUBROUTINE  RETURNS  3 
I  [  (SMALL  BOTTOM  ANGLE  PREFIX  S,  3 
I  [BIG  BOTTOM  ANGLE  PREFIX  P)  3 


I  [  V21  ,V32  =  RFFLECTI0N  COEFFICIENTS} 


I  [AT  INTERFACE  3 

I  [  SIN2=SIN  OF  GRAZING  ANGLE  IN  3 

I  [HIGH  SPEED  LAYER  3 

I  [  TRF=TR ANSNISSION  COEFFICIENT  3 

I  [FACTOR  3 

I  [  GI=GEOMETRIC  INTENSITY  3 

I  [  L  S  =  A  R  C  LENGTH  OF  RAY  IN  THE  3 

I  [BOTTOM  3 

I  [  T=TRAVEL  TI*E  3 

I  [  DTS ,PTR=DELTA  T  FOR  SURFACE  3 

I  [REFLECTION  FOR  S OUR C E , P E C E I  V E R  3 

I  [RESF.  3 

I  [***********************+********3 

I  [***  3 

I 

I 
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•  •••••••••••••••••••••••••••••••••••••••A 

PI=3. 141 592654 
NT  IR  =T  IR  *1  80/PI  ♦! 
NTHC=THC*180/PI+1 
N*NTHC~NTIR 
DE  L  =PI / 1 80 
TH= (NTHC-1 ) *  D E L 
7l=RHl*Cl 
Z2=RH2*C2 
Z3  =  RH3*C3 
FAC  2  =  C  2 /C  1 
FAC  3  =  C  3  /C  1 
FAC  S  =  C  S  /Cl 
$EST=THRC 
BE$T=THRC 
D 1 =2  *ZP”Z  h~Z  S 
D2=2*C3**2/G/C1 
D3  =  2*C3/6*BETA/(1<-8ETA> 
D4=BETA/(1*BETA)/G 
07  =  1  *B  E  TA 
D5=D1 / Cl 
06  =  2/ C  S 
ANU=BETA*C3 
SIG*G*(1+BETA)*C3 
AP»U  =  C3**2*(1  +  bETA>**2 
D8=2*C3/G/(1*3ETA)/AlpHC 
D9  =  -GAM*C3**2/SIG/6 

Dl0=(PETA*GAM*C3**2*GAM*ANU*C3*2)/4/SIG 

Pi  1  =>3*09 

D12=(GA**(ANU**?-A*U>+2*ALPH0*SIG*2*BFTA* 

ANU*C3*GA^) 12  /S  IG 

D14=BFTA*GA*,*(ANU**?-AKU)/2/SIG*PFTA* 

ALPHC 


I 

I 

: 

DO  10  1=1,90 

: 

I 

I 

: 

PTRF (I ) =0 

: 

I 

c 

1  CD 

I 

: 

STR  F ( I ) =0 

i 

•  T 

I 

I 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


>  :  DO  100  1=1  f N 


I 

1 

J=NTHC-1 

R=D1 /TAN (TH) 

I 

1 

:  :  CALL 

REFPN(TH*SEST»9EST,SA,SA3,BA, 

•  • 

:  :  B  A  3  ) 

•  • 

•  • 

■ .  ■ . 

•  • 

I 

I 

/ . . . . . . . . . \  TRUE 

<  I F (P A  .GT.THlB ) GO  TO  998  > - 

\ . / 

I  FALSE 

I 

I 


BC1  -COS(cA) 

BS1*$QPT(1-3C1 **2) 
BC2=FAC?*PCl 
eS?=SGPT(1-oC2**?) 

BC3=FAC3*PC1 
FS3=SQRT (1-BC3**?) 


I 

I - C  BIG  ANGLE  REFLECTION  COEFFICIENT 

I 

I 


A 1 =Z1  /BS1 
A?=Z?/PS? 

A3=Z3/BS3 

BV3?(J)=( A3-A?)/(A3*A?) 

B V  ? 1  (J)=(A2-A1  )/( A2  +  A1) 
BTRF(J)=16*A1*A3*(A?/(A34A2)/(A?4A1))**? 
PSIN?f J)=BS? 


I 

I - [  GEOMETRIC  AMPLITUDE  1 

I 

I 


DRDT=-Dl/PS1**?*D?*CB$1/eS3/BC3**?>*(1.* 
E A3*  TANCBA3) / D7 ) 

BGI  (J)=PCT/ (R*PS1*ADS(DRDT)) 
PCI(J)=SQRT(LGI(J)) 


•  o 
7 
7 
7 
7 
7 
7 
7 
7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


I 

1 

I - C  ARC  LENGTH 

3 

7 

I 

7 

I 

7 

7 

:  BLS(J)=D8/BC3*(D9/BC3**3*BS3**3+D10/BC3* 

L 

7 

:  BS3+(D11/PC3**3-*D12/BC3)*BS3+(D10/BC3**2* 

7 

:  DU)*BA3) 

7 

j 

I 

7 

I - C  TRAVEL  TIME 

3 

7 

I 

7 

I 

7 

l 

:  BT(J)=D5/PSl*D4*(ALOG((1*BS3)/(1-es3))+2* 

7 

:  BA3/PC3/BETA) 

7 

2 

I 

7 

I - C  DELTA  T  SOURCE 

RECEIVER 

3 

7 

I 

7 

I 

7 

2 

:  BSS=SQRT  0-(FACS*ECl)**2) 

7 

:  BDTR  (J )  =  D6* ZR*PSS 

7 

:  PDTS  (  J  )  =  Dfc  *  ZS  *  PS  S 

7 

2 

I 

n  r _ __  __  __ 

7 

C  9983  I 

SC1=C0S(SA) 

• 

I 

I 

\  TRUE 

/ . . . . . . . 

<  I  F ( SA  .GT. TH1 B ) GO  TO  9©9 

-  - 

N,  _  _  _ 

X 

V 

/ 

l 

I  FALSE 

7 

I 

7 

I 

7 

7 

:  AUX=FACS*SC1 

• 

7 

. 

• 

7 

I 

7 

I - C  IF  PAY  DOES  NOT 

REACH  SURFACE 

3 

7 

I  [NOT  COMPUTED 

3 

7 

I 

7 

I 

7 

/ . 

<  I F ( AUX .GT .1 . ) GO  TO  999 

\  TRUE 

7 

^  —  —  —  —  —  —  —  —  —  —  —  —  — 

/ 

7 

I  FALSE 

7 

I 

7 

I - [  DELTA  T  SOURCE 

RECEIVER 

3 

7 

I  7 

I  7 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 


I  7 

I C  $*ALl  ANGLE  REFLECTION  3  7 

I  [COEFFICIENT  3  7 

I  7 

I  7 


7 

7 

7 

7 

7 

7 

7 

7 

7 


I  7 

I - [  GEOMETRIC  AMPLITUDE  D  7 

I  7 

I  7 


7 

7 

7 

7 

7 

7 


1  7 

I - [  ARC  LENGTH  1  7 

I  7 

I  7 


7 

7 

7 

7 

7 


I  7 

I - C  TRAVEL  T I  *  E  D  7 

I  7 

I  7 


7 
7 
7 
7 

I  7 

0  - C 

C  990]  I 


ST(J)=C5/SS1*D4*(AL0G((1+SS3)/(1-SS?>>+2* 

SA3/SC3/BFTA) 


SLS(J)=D8/SC3*(r9/SC3**3*SS3**3*D10/SC3* 

SS3+(M1/SC3**3*012/SC3)*SS3*(D10/SC3**2* 

D14)*SA3) 


DRDT  =  -Cl/SSl**2+D2*(SS1/SS3/SC3**2)*n  + 
S A3*TAN ( S A3) / P7  ) 

SGI(J)=SC1/(R*SS1*AoS(dRDT)) 

SGI ( J )=SGPT ( :gi ( J ) ) 


A1  =71  /  S  S 1 
A2=Z2/SS2 
A3  =  Z3/SS3 

SV3?(J)=(A3-A2)/(A3+A2> 

SV21 (J)*(A2-Al )/<A2*A1) 
STRF(J)*16*A1*a3*<a2/<a3+A2>/<A2*a1))**2 
SSIN2(J)=SS2 


SSS=SQRT(1-AUX**2> 
SDTR (J)=D6*ZR*SSS 
SDTS (J  >  =  D6* ZS*SSS 
SSl*SQRT(1-SCl**2) 
SC2=FAC2*SC1 
Ss2=$QRT (1-SC2**2) 
SC3=FAC3*SC1 
SS3=SQRT ( 1 -SC3**2) 
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( 


A 
A 
A 
A 

(  A 


( 


( 


m 

• 

TH=TH-DEL 

• 

• 

• 

I 

c 

1003 

I 

CONTINUE 

• 

• 

•  ,  , 

• 

I 

I 

I"" . . \ 

:  RETURN  : 

\ . / 
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\  SUBROUTINE  BASE  / 

\ . / 


I 

1  c*** ******* ************** *««*«««« 3 
I  [***  3 
I  C  FOR  EACM  SPECULAR  BOTTOM  ANGLE  3 
I  [SUB  RETURNS  3 
I  C  BSV21 *BSV32*RE FLECTION  3 
I  [COEFFICIENTS  AT  INTERFACE  3 
I  [  BSSIN2=SIN  OF  GRAZING  ANGLE  IN  1 
I  [HIGH  SPEED  LAYER  3 
I  C  BSTRF=TRANSPISS ION  COEFFICIENT  3 
I  [FACTOR  3 
I  [  BSGI-6E0METRI C  AMPLITUDE  3 
I  [  BSLS=ARC  LENGTH  OF  RAY  IN  3 
I  [BOTTOM  3 
I  [  BST=TRAVEL  TIME  3 


I  [  BSDTS ,BSDTR=DELTA  T  FOR  SURFACE3 

I  [REFLECTION  FOR  SOURCE  RECEIVER  3 

I  [***  3 

I 

I 


COMMON  Ze,ZR,ZS,Cl,C2tC3,PHl,PH2,PH3,G,D, 
ALPHO.CS.PZ  *  R  E  F 

COMMON  /DATA/LZERO,FZERO,L70,F70,TH70, 

L71 , F  7 1 , DPLOSS • FREQ  t  THtTA  « A  LOS  S  0 1  ALOSSM  « 
FM,TIR,THC,THPC,THlP,Bt.TA 

COMMON  /BS/PSV21,PSV32,BSSIN2,BSTPF,BSGI, 
BSLS, PST, BSDTS, PSDTP 

DIMENSION  BSV21  (90>,BSV32(90),PSSIN2(90>, 
BSTPF(90),FSGI(90),PSLSC9n),BST(90), 
BSDTS(90) *  E  S  D  T  P (90) 

COMMON  /ALPH/GAM 
REAL  L7ERO  ,L7C,L71 
PI =3.1 41 592654 
NTIR=TIR*1 8  C / P I ♦ 1 
DEL=PI /  180 
T  H  =  P 1  It 
Z1 =RM1 *Cl 
?2=RH2*C2 
?3=RH3*C3 
FAC2=C2/C1 
FAC3  =  C3  /C 1 
FACS=CS/C1 

CP=$QRT(2*G*(1+BETA)*C3*BZ+(C,*(1*PETA))* 

*2)-PETA*C3 

FACP=CB  /Cl 
M  =2*  ZR-Z k-2 S 
D2=2*C3**2/C1 /G/C1+BETA) 

D3=2*C?/G  /(  1+BETA  ) 
D4=?/G/(1+BETA) 

D5*D1  /  Cl 
D6  =  2/C  S 

D7=2*C3*PETA/G/ (1*BETA)/C1 
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I 

I - C  THE  FOLLOWING  2  LINES  ADDED  5/  ] 

I  C  29/ 79  ] 

I 

I 


•  •••••••••••••••••••••••••••••••••«•••••• 

DA=2.*C3/6/ ( 1 •♦BETA) 

DB  =  D  A*B  ETA 
ANU'BET  A  *C3 
S I G=  6* (1 ♦BETA)*C3 
AHU«C3**?*(14BFTA)**? 
D?=?*C3/G/(1+8ETA>/ALPHO 
D9=-6AM*C3**2/SIG/6 

Dl0=(BETA*GAM*C3**2+GAM*ANU*C3*2)/4/SIG 

D1 1 =-3*D9 

D12=(GAM*(ANU**2-AMU)+2*ALPH0*SIG+2*BETA* 
ANU*C3*GAP) /2/S IG 

Dl4rBETA*GAP*(ANU**?-ANU)/2/SIG+BETA* 

alpho 

B  SA=TH-1 . E-6 

TH3  =  ACOS( FAC  3* COS (TH1B)  ) 

R  =  Dl/TAN(THlB)*D3/COS(TH3)W*(TH3/?+(. 
?5+PETA/2)*SIN(2*TH3) ) 

AU  7  =  A  T  A  N  (  D  1  /  P  ) 

N=  AUX*1 80 /PI +1 
N=90-MAX0(N,NTIR) 


I 

I 


A - > : 

DO  10  1=1,90  : 

A 

. 

A 

I 

A 

I 

A  .  .  . 

A  : 

PSSIN2(I)=0.0 

A  : 

PSV32(I) =0.0 

A  : 

BSV21  ( I )  =0. 0 

A  : 

BSG I ( 1 ) =C.O 

A  :  .  . 

A 

I 

A  l 

100 

I 

A  .  .  .  < 

PSTRF(I)=C  : 

I 

I 

A - >  : 

DO  10°  I =1 ,N  : 

A  , 

• 

A 

I 

A 

I 

A  .  .  .  , 

A  : 

J  =  9  0  ~  1  : 

A  : 

TH=TH-DtL  : 

A  :  .  .  , 

A 

I 

A  I 

A 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


CALL  PA  SANG  (TH  tBSA  «8SA3) 


I 

1 

j  +  +  +  +  +  +  +  +  ++  +  +  +  ++  +++*  +  +  +  +*  +  ++  +  **+  +  +  +  +  +  +  +  +  +  +  y 

<  I F (BSA3. EQ.O.O)  GO  TO  100  > 

\ . / 

I  FALSE 

I 

I 


TRUE 


B  Cl =COS (bSA) 

PSl  =  SQRT(1-dCl**2) 
BC2-FAC?  *BCl 
PS2=SQRT ( 1 -6C2**2) 
BC3=FAC3*BC1 
BS3=SQRT(1-BC3**2) 
BCB=FACB*BC1 
PSP=SQRTM-bCB**2) 


I 

I - C  DELTA  T  SOURCE  RECEIVER  3 

I 

I 


AUX=FACS*BC1 


I 

I - [  IF  RAY  DOFS  NOT  REACH  SURFACE  3 

I  CMOT  COWPUTFD  1 

I 

I 

/ . . . . . .  \  TRUE 

<  IF  (AUX  .GT.1  .)  60  TO  999  > - P 

\ . /  X 

I  FALSE  X 

I  X 

I  X 

. . y 

:  SS=SQRT(1-AUX**2)  :  x 

PSDTP(J)=C6*7R*SS  :  X 

:  PSDT$(J)=D6*ZS*SS  :  x 

. . x 

I  x 

I - C  LAYER  REFLECTION  COEFFICIENTS  1  X 

I  X 

I  X 

. X 

:  *1  =  Z1/BS1  :  x 

:  A2=Z?/B32  :  x 

:  A3  =  Z?/Bb3  :  x 

:  BSV32  (J  )= ( A3-A2) / (A3  +  *2)  :  x 

:  ESVtl (J )  =  (A2~Al)/ C  A  2 ♦ A 1 )  :  x 

:  BSTRF  (J)  =  16*A1* A3*(A2/ ( A 3 ♦ A 2 ) / ( A 2 ♦ A  1 ) )**  :  X 

:  2  :  x 

:  PSSIN'2(J)=BS2  :  x 

. . x 


r. 

7 

7 

7 

Z 

7 

7 

Z 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

Z 

Z 

7 

7 

7 


Z 

7 

z 

*T 

7 

7 

7 

7 

7 

7 

7 

7 

7 

Z 
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A  I  X  7 

a  i — C  geometric  amplitude  3x7 

a  1  x  z 

a  1  x  7 

A  . X  7 

A  :  PS AB'ACOS  (BCB)  :  X  7 

A  . . - . .  X  7 

A  I  X  7 

A  I - C  FOLLOWING  4  CARDS  ADDED  5/59/793  X  7 

A  I  X  7 

A  I  X  7 

A  .  X  7 

A  :  DRDT3=DA* (B$3/BC3**3* ( BSA3-BS AB*BS3*BC3-  :  X  7 

A  :  BSB*BCB>+(1 . - ( C P/ C3 > *  * 2*B S3 / BS B ) )*DB/PC3*  :  X  7 

A  :  (1./BC3>BS3*(BSP/BC3«(CB/C3)**2*BC3/9SB))  :  X  7 

A  :  DRDTs“Dl/PSl**2+(C3/Cl)*BSl/B$3*DRDT3  :  x  7 

A  :  R=Dl/TANCTH)  :  X  7 

A  :  PSGI (J)*PCT/ (R*BS1*ABS (DRDT) )  :  X  7 

A  :  BSGI ( J )=SQRT ( BSGI ( J) )  :  X  7 

A  . . x  7 

A  I  X  7 

A  I - C  ARC  LENGTH  3X7 

A  I  X  7 

A  I  X  7 

A  .  X  7 

A  :  PSLS(J)=Df/PC3*fD9/PC3**3*(PS2**3-PS9**3)  :  X  7 

A  :  ♦DlC/PC3**2*(PS3*FC3-BoB*BCB)*(Dl1/PC3**  :  X  7 

A  :  3*D12/PC3>*(BS3-BSB>*(D10/BC3**2*M4)*(  :  x  7 

A  :  BSA3-PSAB))  :  X  7 

A  . . x  7 

A  I  X  7 

A  I - C  TRAVEL  TINE  3X7 

A  I  X  Z 

A  I  X  7 

A  .  X  7 

A  :  BST ( J ) =D5/PS1 +D4* ( (PS A3-BSAF) /BC3+BET A*  :  X  7 

A  :  AL0G((1*BS3)/(1-BS3>*(1-BSB>/<1*BSB>)/2>  :  X  7 

A  . . x  7 

A  I  X  7 

A  0< - 0 

A  C  10PD  I  X 

A  .  X 

- :  CONTINUE  :  X 

. . x 

I  x 

0< - 0 

C  9993  I 


RETURN 
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\  SUBROUTINE  REFRN(THfSEST»BEST»S^tSA3,  / 
\  BA  »B A  3  >  / 

\ . . . / 

I 

I 


COMMON  ZB,ZR,ZS,C1,C2,C3tRHl,RH2,RH3,G,D, 
ALPHO,CS,PZ  » R E  F 

COMB ON  /DATA/LZERO,FZER0,L70,F70,TH70, 

L71  , F71 tDBLOS St  FREQ, THETA ,ALOSS Ot ALOSSM, 
FM,TIPtTHC  ,THPC  ,TH19  .BETA 

REAL  LZEPO,l70 ,L71 


I 

I  — ~[********************************3 

I  [***  3 

I  C  6IVEN  3 

I  [  T I R= ANGLE  OF  TOTAL  INTERNAL  3 

I  [REFLECTION  3 

I  [  T  H  =  G  R  A  Z I NG  ANGLE  3 

I  [  SEST=ESTIHATE  SMALL  ANGLE  3 

I  [  6EST=ESTI*ATE  BIG  ANGLE  3 

I  [  RETURNS  3 

I  [  S  A  =  SM ALL  ANCLE  REFRACTED  RAY  3 

I  C  B  A  =P I G  ANGLE  REFRACTED  RAY  3 

1  [***  3 

1 

I 

/ . """" . "  " . . ' . . \  TRUE 

<  IF (TH.GT.TIR)GO  TO  100  > - 

\ . / 

1  false 
1 
1 

/ . / 

/  WRITE (1 ,1 )  / 

/ . / 

1 


1 [  1  FORMATCo  SPECULAR  ANGLE  OF  3 

1  [REFPACTED  RAY  LESS  THAN  TIRn)  3 

I 

I 


/ . \ 

:  STOP 

\ . / 


0< 

r  ico3  i 


o 

7 

7 

7 

7 

7 

7 

7 

7 

7 

Z 

7 

7 

7 

7 

7 

7 

7 

0 
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Dl=2*ZB-ZR-ZS 
PI2=3. 141592654 
Pl2=Pl2/2 

D2=2*C3/G* (BETA/(1 ♦BETA) ) 

D3  =1 / ( 2*B  ET  A  ) 

FAC  3  =  C  3 /C 1 
R=D1  /TAN  (TH) 
T3=ACOSfFAC3*COS(SEST)) 

VAL  *D1/TAN($E$T)+D2*TAN(T3)*(1M14T3/ 
SIN(T3)/COS(T3) >*D3)-R 

DEL= (SEST-TI R )/1 0.0001 


I 

I 

I 

C  1023  I 


A - >  :  DO  200  1  =  1,500 

A  . . . 

A  I 

A  I 

A  /  ’"\  TRUE 

A  <  I  F  C  APS  (  DEL  ).LT  .  1  .  E-6  )  GO  TO  201  > - 0 

A  \ . /  7 

A  I  FALSE  Z 

A  I  Z 

A  I  Z 

A  .  7 

A  :  SA  =  SEST-[)EL  :  7 

A  . . .  Z 

A  I  Z 

A  I  7 

A  /" . . . TRUE  ? 

A  <  I  F  (  S  A  .  LE  .TI  R  )  oO  TO  211  > - 0  7 

a  \ . /  x  7 

A  I  FALSE  X  7 

A  I  X  7 

A  I  y  7 

A  / . '\  TRUE  y  7 

A  <  IF(SA.GE.THcC)GO  t0  211  > - V  7 

a  \ . /  y  7 

A  1  FALSE  X  7 

A  1  X  7 

a  i  y  7 

a  .  y  z 

A  :  T3=ACCSC FAC3* COS CSA) )  :  y  7 

A  :  V ALP=D1 /T AN (S A  ) ♦ D2 * T A N ( T 3 ) * ( 1 + ( 1 ♦T 3 / S I N C  :  X  7 

A  :  T 3 ) /COE( T 3 ) ) *  C 3  )-R  :  x  7 

A  :  AUX  =  V*L*  wALp  :  y  7 

A  :  S  E  S  T  =  S  a  :  y  7 

A  :  VAL=VALP  :  y  7 

A  . . .  x  7 

a  i  y  7 

A  I 
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7 


A 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


/ . . . . \  TRUE 

<  IF(AUX  •6T#0»)GO  TO  200  > - 

\ . . . .  ... . / 

I  FALSE 
I 

C  21  CD  I 


•»»...««»» 

DEt*-DEL/3 


I 

I 

/ .  . \ 

:  GO  TO  200  : 

V . / 


0< 

[  2113  I 


PLL  =  PEl  /3 


I 

0< 

[  2  OC  3  I 


x 

X 

0  X 

U  X 

U  X 

U  X 

U  X 

u  x 

U  X 

U  X 

U  X 

u  x 

U  X 

V  x 

l)  X 

U  X 

l'  X 

- 0 

u 

u 

It 

u 

u 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 


CONT I NUE 


I 

I 

/ . . . / 

/  WRITE  (1,2)  / 

/ . / 


I 

1 - C  2  FOR^AT(o  SMALL  BOTTOM 

I  [REFRACTION  ANGLE  DOES  NOT 
I  [CONVERGE^) 

I 

I 

/ . \ 

:  STOP  : 

\ . / 


0< 

[2013  I 


S  A  3=  T  3 


3 

3 

3 


7 

7 

i 

7 

7 

7 

7 

7 

7 

7 

7 

7 


7 

7 

7 

7 

7 


I 

I - 1  BA  INITIALIZATION  D 

I 

I 
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T3*ACOSfFAC3*CO$(BEST)> 

VAL  =Pl/TAN(eEST)*D2*TAN(T3>*(1+( 1+T3/ 
SIN(T3) /CCS (T3) >*D3)-R 

DEL=<PI2-BEST)/10.0001 


I 

I 

I 

C  3  02  D  I 


A - >  :  DO  400  1*1* 50C  : 

A  . . . 

A  I 

A  I 

A  . .  . \  TRUE 

A  <  I F ( ABS ( DEL ).LT. 1 .  E-6  )  60  TO  401  > - 0 

A  \ . . . /  7 

A  I  FALSE  7 

A  I  Z 

A  I  Z 

A  .  Z 

A  :  B A=BEST+DEL  :  Z 

A  . . .  Z 

A  I  Z 

A  I  Z 

A  / . .  . . . . \  TRUE  Z 

A  <  IF(BA.LF.THRC)GO  TO  411  > - 0  7 

A  \ . /  X  Z 

A  I  FALSE  X  Z 

A  I  X  Z 

A  I  X'  '  7 

A  / . \  TRUE  X  7 

A  <  I  F  (bA  .GE.PI2)  GO  TO  411  > - V  Z 

A  \ . /  X  Z 

A  I  FALSE  X  ? 

A  I  X  Z 

A  I  X  Z 

A  .  X  ? 

A  :  T3=ACOS(FAC3*COS(BA))  :  x  Z 

A  :  VALP  =  M/  TAN(BA  >  +  D2*TAf*(T3)*(1  +  <1+T3/SIK(  :  X  Z 

A  :  T3)/crs(T3))*t3)-P  :  X  Z 

A  :  AUX=VAL*VALP  :  X  Z 

A  :  B  E  S  T  =  B  A  :  X  Z 

A  :  VAL=VALP  :  x  z 

A  . . .  X  Z 

A  I  X  Z 

A  I  |  X  Z 

A  / . \  TRUE  X  Z 

A  <  JF(AUX.GT.O.)oO  TO  400  > - C  X  Z 

A  \ . /  U  X  Z 

A  I  FALSE  U  X  Z 

A  I  V  X  Z 

A  C  4 1C]  I  U  X  Z 
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A 

A 

• 

• 

DEL=-DEL/3  : 

u 

u 

X 

X 

z 

7 

A 

• 

. . • 

u 

x 

1 

A 

I 

u 

X 

Z 

A 

I 

u 

X 

z 

A 

/' 

u 

X 

z 

A 

• 

• 

GO  TO  400  - - 

X 

z 

A 

\  . 

u 

X 

z 

A 

u 

X 

z 

A 

A 

A  ^  —— —  ________ —  _ 

U 

X 

z 

7 

A 

L 

A 

t 

4in 

I 

u 

1 

A 

2 

A 

: 

DEL  =  DEL/3  : 

s 

7 

A 

. 

• 

u 

2 

A 

I 

u 

7 

A 

0  < - 1 

Z 

A 

c 

4003 

I 

2 

A 

2 

-  ; 

CONTINUE  : 

7 

• 

• 

2 

I 

7 

I 

7 

/  ' 

. . / 

7 

/ 

WRITE  Cl  ,3)  / 

7 

/ 

. . / 

7 

I 

Z 

I - [  3  FORMAT ( n  LAPGE  BOTTOM 

3 

7 

I  [REFRACTION  ANGLE  DOES  NOT 

3 

7 

I  [CONVERGED) 

3 

7 

I 

7 

I 

7 

r 

. . . \ 

7 

: 

STOP 

7 

\. 

7 

7 

7 

f>< - 

[  401  3 

I 

P  A  3=  T  I 

z 

•TTTTTTTtTTTTT 

>rrrrrrrrr* 

I 

I 

/ . . \ 

:  RETURN  : 

\ . / 
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\  SUBROUTINE  B A S » NG  (TH  , B S A  , B S A3 )  / 

\ . / 

I 

I 


COMMON  ZB,ZR ,ZS *C1 ,C2,C3,RHl ,RH2, RH3,G,D, 
ALPHO.CS.BZ.REF 

COMMON  /DATA/LZER0,FZER0,L70 ,F70,TH70, 

L71 , F71 ,DBLOSS, FREQ, THETA, ALOSSO.ALOSSM, 
FM,TIP,THC  ,THRC  ,TH1B,B  ETA 

REAL  LZER0,L70 ,L71 


/ 


*  +  +  + 


< 


\ 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

1 

I 

I 

I 


[**************** 


[***  3 

C  GIVEN  3 

C  BZ=SEDI*ENT  DEPTH  3 

C  TH1B=BASEMENT  CRITICAL  ANGLE  3 

C  TIR=ANGLE  OF  TOTAL  INTERNAL  3 

[Reflection  3 

[  TH=GRAZING  ANGLE  3 

[  RETURNS  3 

C  BSA=BASCMFNT  REFLECTED  RAY  3 

[ANGLE  IN  WATER  3 

[  BSAT=BASEMENT  REFLECTED  RAY  3 

[ANGLE  IN  SEDIMENT  3 

[  +  *+  +  +  **+*+******  +  +  +  +  +**+++*****  +  3 

[******♦****♦*+**+***+*+*****+***3 
[***  3 


"""  \  TRUE 


IF (TH.6T.TIR)G0  TO  100  > 


/ 


0 


I  FALSE 


I 

I 

/ . . . . . . . . . . ""/ 

/  WR I T  E  ( 1 , 1 )  / 

/ . / 

I 


I - [  1  FORMAT ( n  SPECULAR  ANGLE  OF 

I  [BASEMENT  RAY  LESS  THAN  TIRn) 

I 

I 

:  STOP  : 

\ . / 


0< 

[  i or  3  i 


3 

3 


7 

7 

Z 

z 

7 

7 

Z 

Z 

Z 

Z 

7 

7 

7 

* 

z 

7 

7 

0 
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D1=2*ZB-ZR-ZS 
PI *3.141 592654 
P I 2=P I /2 

D2=2*C3/G/(1*BETA) 

FAC3=C3/C1 

CB=SGPT(2*G*( 1*BETA)*C3*BZ*(C3* (1+BETA) > * 

FACB=CB/C1 
P  *t>  1  /  TAN  (TH) 

BEST=BSA 
AUX'COS (BEST) 

THP= ACOS ( FACB*AUX ) 
TH3=ACOS(FAC3*AUX ) 

VAL=Dl/TAN(8ESTm2/(FAC3*AUX)**2*((TH3- 

DEL  =  F>i  /  180 
A UX  =  T  I  R 


I 

I 

/ . ' . . .  "  ' .  \  FALSE 

<  IF (THlB.6I.TIR)  > - 

\ . / 

I  TRUE 
I 


AUX=TH1B 


C  1013 


>  +  +  +  +  +  +  ++  +  +  +  +  +  +  + 


I 

0< - 

I 


TRUE 


<  IF ( (BEST-DEL) .6T.AUX)60  TO  IP?  > 

\ . .../ 

I  FALSE 

I 

I 


DEL=DEL*.9 


I 

I 

/ . . . \ 

:  60  TO  101  : 

\ . / 


[  102  3 


I - C  COMPUTE  BS  A 

I 

0< - 

I 


I 

I 

I 

I 

I 

I 

I 

I 


0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 


X 

X 


3  x 


X 

0 


0 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

0 


A --> :  DO  200  I s  1  f 5 °C 

A  . . 

A  I 

A  I 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


B  SA  =  BE  ST-DEL 


I 

I 

r . . . . . . TRUE 

<  IMBSA.LE.AUX>G0  TO  211  > - 

\ . .../ 

I  FALSE 

I 

I 

/ . . .  ""  TRUE 

<  IF(BSA.6E.PI2)G0  TO  211  > - 

\ . . . / 

I  FALSE 

I 

I 


AUXl *cos ( es A) 

THP=ACO$(FACB*AUXl) 

TH3  =  A  COS ( FAC2*AUX1  ) 

VALP=D1/TAN(6SA>*D2/(FAC3*AUX1)**?*((Th3- 
B AUX  =V  At*  VAlp 

best=bsa 
V AL=VALP 


I 

I 

. . . .  true 

<  I F (B AUX  .GT .0. ) GO  TO  20C  > - 

\ . . . . . / 

I  false 

I 

C  2 103  I 


0EL=-DEl/3 


\  TPUr 
> - C 

/  u 

V 

u 

u 

u 


u 

V 

V 

c< - 

[  2113  I  U 

.  V 

PEL- DEL/3  :  V 

. V 

I  l' 

0< - 


I 

I 

. . . . . 

<  I F C APS ( DEL  ).LT .1 *E-6>  GO  TO  201 

\ . 

I  FALSE 

I 

I 

r" . . . \ 

:  GO  TO  200  : - 

\ . / 


0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

V 

X 

X 

X 


X 

X 

X 

X 

X 

0 


r\ 

2 

2 

2 

2 

2 

V 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

7 

7 

2 

7 

2 

7 

7 

7 

7 

2 

7 

2 

7 

7 

2 

7 

7 

7 

y 

1 

7 

? 

7 

7 

2 
7 


154 


A  t  ?OPD  1  u 

A  •••••••••••••••••••••••••••••••••••••••••••••  U 

— :  CONTINUE  :  y 

. . .  V 

1  U 

I  u 

.  u 

:  PSA=PI /2.0-1 .OE-6  :  U 

:  TH3=0 • 0  :  v 

: . .  U 

I  u 

I  r 

,.*1  tl 

/  WRITE (1*2)  /  U 

/ . /  V 

I  u 

I - C  2  FORPATCd  BASEMENT  AN6LE  DOES  f  V 

I  [NOT  CONVERGED)  3  U 

I  U 

0  < - 0 

c  2on  i 


6SA3=Th3 


I 

I 


RETURN 
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\  SUBROUTINE  LOSS(JJ)  / 

V . / 

I 

I 


COMMON  ZB,ZR.ZS,C1,C2.C3.RH1,RH2,RH3,G,D, 
ALPH9.CS.BZ , REF 

COMMON  /DATA/LZERO.F  ZERO.L70.F70.TH70, 

L71 ,  F71 , DBLOSS , FREQ , THETA . ALOSS 0, ALOSSK  , 
FM,TIR,THC,THRC,TH1P.BETA 

COMMON /A L/  F  (f  )  ,  ALOSS  (6,90 
COMMON  /R/NSAP*RV21 ♦ R V32 * RS I N2 * R ots ♦ R DTR , 
RT ,  RG I 

DIMENSION  RV21  (90)  ,PV3  2(90)  .RSIN2  (90  » 
RDTS(90),PDTR(90),RT(90),PGI(90) 

COMMON  /B/BV21.BV32.BSIN2.BTRF.BGI.ELS. 
BT.PDTS.BDTR 

COMMON  /S/SV21 .SV32.SS 1N2.STRF  ,  SGI , SLS . 
ST.SDTS.SDTR 

DIMENSION  BV21(«0),PV32(90),BSIN2(90)» 
BTRF(9C),BGI(90),BLS(90),BT(90,BDTS<90), 
BDTR  (90 

DIMENSION  SV21 (90) »SV32(9H) ,SSIN2(90) , 
STRF(9C),SGI(90),SLS(9u),ST(90),SDTS(90), 
S  DTP  (90 

COMMON  /BS/PSV21,BSV32,BSSIN2,BSTPF,BSGI, 
BSLS.PST.FSf'TS.PSDTP 

DIMENSION  ESV21 (90) »BSV32(9C) »BSSIN2(90) , 
BSTPF(90),BSGI(90).PSLS(90),PST(90), 
BSDTS(90)  ,BSDTR  (90 

COMPLEX  AUX99,V,AIM,AUX1,AR(16)*PH(16) 
DIMENSION  toP ( 4 ) 

REAL  LZEPO.L70.L71 


I 

.*..***.*****. ****************1 
I  [*.**.***************************] 

I  [***  7 

I  [  THIS  SUBROUTINE  COMPUTES  3 

I  [APPARENT  LOSS  DUE  TO  PROPAGATION} 

I  [A*******************************]] 

I  [***  3 

I 

I 


J  56 


PI  =3. 1415  52654 

ALPH0=A10G(T0.) /4/PI/1000r*AlPH9 
EX=A10G10(EXP(1.)) 

PI2=2*FI 

A  I M  =C  MPL  X ( 0 • , 1  •  ) 
SQRT2=SQRT<2.) 

F  1  *=F  (  J  J  )  /SORT  2 
F2*F  <JJ)*SQRT2 
FPI2D=4.*PI*F (JJ)*D/C2 
NAVI  *  £ 

FPIINC  =  FPI2D* (SQRT2-1.0/SQRT2)/NAV1 
FPIIN=FPI2P-(NAVl-1)*FPIINC/?. 
NTHC=THC*1 tiO/Pl 
N=90-NTHC 


I 

I C  COMPUTE  LOSS  FOR  ANGLES  ABOVE 

I  [CRITICAL  (CAUSTIC) 

I 

I 


A - >  :  CC  10  1  =  1 ,90 

A  . . 

A  I 

A  [103  I 

A  . 

- -  ALOSS (J J . I )  =  25 . 


I 

I 


A - >  :  DO  100  I  =1  , N 

A  . . 

A  I 

A  I 

A  . 

A  :  J  =  9 1  -  i 

A  :  FPIAUX=FPIIN 

A  :  A  U  x  =  0  • 

A  . . 

A  I 

A  I 

A  . 

A  P - >  :  DC  20  NAV  =  1,NAVl 

A  P  . . 

A  P  I 

A  P  I 

A  P  . 

A  P  :  ALPH2= FPIAUX* nSIN2(J) 

A  P  :  C  =C  0  S (ALPH2) 

A  P  :  S=SIN(ALPH2) 

A  P  :  W=  (  RV21  (  J  )  ♦R  V3  ?(  J  )  *  C  )  **2+  (  °V32  (  J  )  *  S  )  **  2 

A  B  :  VV=VV/((1+RV21(J)*RV3?(J)*C)**2+(RV21(J)* 

AP  :RV3?(J)*S)**2) 

A  p  :  AUX=AUX+VV 

A  P  :  FPIAUX=FPIAUa+FdI INC 

A  p  . . 
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A  B  I 

A  B  C  203  I 

a  e  . . . . 

A  - :  CONTINUE 

A  . . . . 

A  I 

A  I 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


A  p - >  :  DC  110  K-1  ,3 

A  P  . . 

A  P  I 

A  B  I 


VV= AUX/NAVl 
DR*PI2*RDTR< J> 

DS-PI2*RDTS(J) 

AU  X  =  Dfi ♦ DS 

DT=SIN(F?*AUX)-SIN(F1*AUX) 

DT  rDT / AUX 
A  U  X  =D S-DR 

DT=DT+SIN(F2*AUX)/AUX-SIN(F1*AUX)/AUX 
DT=DT-2*($IN<F2*DS)-SIN(F1*DS>)/D$ 
DT=DT-2*($IN(F2*DR)-SIN(F1*DR))/0R 
DT=DT /  2  /  F 1  ♦  1 
DT=VV*  DT 

AtPH2=FPl2D*BSSIN2(J)/2 

C=C0S(ALPH2) 

S=SIN (ALHH2) 

V=BSTRF(J)/CWPLX( (1*BSV32(J)*PSV21 (  J  )  )  *  C  » 
(BSV32(J)*BSV21  f  J  )  -1  )  *  S)  *  *2 

V  V  r  C  A  B  S  (V  )  *  P  S  G I  <J)*REF 
A  U  X  r  0 


I 

I 

/ . T*UE 

IF (VV. EQ.O. )G0  TO  120  > - 

\ . . . . . / 

I  FALSE 

I 

I 


BR  <  1  )  =  VV 
PR ( 2 ) =- VV 
BR  ( 3 ) =- VV 
BR  <4 )  =  VV 

V=CMPLX(BST(J)fALPHQ*BSLS<J)) 

P  H ( 1 ) =p I  2* V 

PH(2)=PI2*(V*3SDTS< J)  ) 
PH(3)“*l2*(V*bSDTR (J) ) 
PH(4)=PI2*(V*BSDTR  (J)^BSDTS(J)) 

v=o 


I 

I 


n 

7 

7 

Z 

Z 

Z 

Z 

7 

Z 

Z 

z 

z 

7 

7 

7 

Z 

z 

z 

7 

z 

7 

7 

7 

7 
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A  B  Z 

A  B  :  NA=K*1  :  z 

A  B  Z 

A  B  I  Z 

A  B  I  Z 

A  B  Z 

A  B  C-->:  DO  110  L=NA,4  :  Z 

A  B  C  . . z 

ABC  I  Z 

ABC  I  Z 

A  B  C  .  Z 

ABC  :  AUX=BR(K)*BR(L)  :  Z 

ABC  :  AUX1=CPH(K)-C0NJG(PH(L)))*AI^  :  Z 

A  B  C  . . Z 

ABC  I  Z 

A  P  C  I  Z 

ABC  I" . . \  TRUE  Z 

ABC  <  IF(REAL(AUXl)*F2.LT.-40.)GO  TO  110  > - 0  Z 

ABC  \ . . . /  X  Z 

ABC  I  FALSE  X  Z 

ABC  I  X  Z 

ABC  I  X  Z 

A  B  C  .  X  Z 

ABC  :  V  =  V* AUX* ( CEXP( F2* AUXl  ) -CEXP (  FI  * AUX1 ) ) /  :  X  Z 

ABC  :  AUXl  :  X  Z 

ABC  : . .  x  z 

ABC  I  X  ? 

A  P  C  0  - - 0  Z 

ABC  C  1 1 C  D  I  Z 

A  B  C  .  Z 

A - :  CONTINUE  :  Z 

A  Z 

A  I  Z 

A  I  Z 

A  Z 

A  :  AUX=2*»E AL<V)  :  7 

A  :  AUX2=-2*PI2*ALPH0*PSLS(J)  :  Z 

A  Z 

A  I  Z 

A  I  Z 

A  / . V  T  R  U  c  Z 

A  <  I  F  (  AUX2*  F  2.  LT  .-40.  )  GO  TO  120  > - V 

A  \ . /  Z 

A  I  FALSE  7 

A  I  7 

A  I  Z 

A  Z 

A  :  AUX=AUX+4*(EXP(F2*AUX2)-EXP(M*AUX2)>/  :  Z 

A  :  AUX2*PF (1 )**2  :  Z 

A  Z 

A  I  Z 

A  0< - 0 

A  [  12T]  I 


159 


A 

A 

A 

A 

A 

A 

A 


DT=DT+AUX/(2*RGI(J))**?/M 
ALOSS ( J J . J>=-10. *A LOG  10  COT) 


C  1001 


A  — 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


• 

9 

CONTINUE 

• 

9 

9 

I 

I 

9 

9 

• 

NM  A  X  =91  -N 

9 

9 

^  9 

I 

I - C  COMPUTE 

LOSS 

FOP 

I  CT1R  AND 

I 

I 

THS 

NTIR=TIR*180/PI+1 

N  =NM  A X -NT  I R 

I 

I 

• 

9 

DO  200  I  =1  »N 

• 

• 

• 

I 

I 

• 

J  A  X-  I 

• 

IF  (J  .LT.NSAP)uO  TO  800 
FALSE 


\  TRUE 
> - 

/ 


I 
I 

I - C  C  0  M*>U  T  E  REFLECTION  COEFFICIENT  ] 

I 

I COMPUTE  FPEG  AVERAGE  ABSOLUTE  ] 

I  CREFLECTION  ] 

I 

I 


A  : 

F  PI AUX=FPIIN 

A  : 

AUX  =  0 

A  : . 

A 

I 

A 

I 

A  . 

A  R —  >  : 

DO  130  N A V  = 1 , N A V 1 

; 

A  e  . . 

•c 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 


B 

B 

B 

B 

P 

B 

P 

B 

B 


ALPH2=FPIAUX*RSIN2( J) 
C=C0S(ALPH2) 

S-S  IN(ALPH2) 

VV=(RV21(J)+RV32(J)*C)**2+(RV32(J)*S)**2 
VV*VV/ ( (1 +  RV21 ( J)  *RV32  <J)*C)**2+ (RV21 (J) * 
RV32 (J)*S>**2) 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


A  B  : 

AU  X  *  A  UX  ♦VV 

A  B 

FPIAUX=FPIAUX+FPHNC 

A  B  . . . 

A  B 

I 

A  B  C  1 30D 

I 

A  B  . . 

A - : 

CONTINUE 

• 

• 

A  . . . 

+ 

A 

I 

A 

I 

A  . 

VV=  AUX /NAVI 
ALPH2=FPI?D*RSIN2 (J) 

C=CPS ( ALPH2 ) 

S  =S IN ( ALPH2 ) 

V=CPPLX(RV21(J)*RV32(J)*C,RV32(J>*S> 

V=V/CKPLX(1*RV2T(J)*RV32(J)*C»RV21(J)* 

PV32(J)*S> 

V=V/CAB$(V)*RC-I(J)*SQRT(VV) 

AR (1 )=V 
AR (2 ) =-V 
AR(?)=-V 
AR (4 )  =  V 

PH(1)=PI2*RT(J) 

PH(2)=PI2*(RT(J)ARDTS(J)) 

PH(3)=ri?*(RT(J)+RDTR(J)) 

PH(4)=PI2*(RT(J)4RDTP(J)*R0TS(J)) 


I 

I-' 

I 

I 

I 


COMPUTE  BIG  TRANSMISSION 
CCOEFFIC I  ENT 


A  R  (  5  )  =  u 
AR  (6)=C 
A  R  (  ?  )  =  0 
AR  C  *  >  =  C 


I 

I 

/ . .  \  TRUE 

<  IF(BTRF( J) .E0.0.)60  TO  205  > - 

\ . / 


A 

A 

A 


I  FALSE 

I 

I 


D 

D 


0 

y 

x 

y 

y 


7 

7 

7 

■* 

i 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

4. 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 


161 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


/ . . \  TRUE 

<  IF(BT(J).LT.RT(J))GO  TO  205  > - 

V . / 

I  FALSE 

I 

I 


ALPH2»FPI2d*8SIN2(J  > /2 
C=COS (ALPH2) 

S=SIN(ALPH2) 

V=BTRF( J ) /CWPLX f ( 1 +  PV3 2 ( J ) * BV2 1 (J))*C«( 
B V32 ( J )*BV21 ( J)-1 >*S>**2 

VSV*PGI  (J ) 

A  R  ( 5  )  =  V 
AP  (  6 )  =  -V 
AR (7) =-V 
AR (8 ) = V 

V=CPPLX(6T(J) ,ALPHO*BLS<J) ) 
PH(5)=PI2*V 
PH(6)=PI2*(V+6DTS(J)> 
PH(7)=PI2*(V+&DTR(J)) 
PH(8)=PI2*  CV  +  BDTS (J)*BDTP(J) ) 


I 

o< - 

[  2C5  3 

i 

z 

CONTINUE 

: 

m99999999m 

I 

I - C  COMPUTE  SPALL  TRANSMISSION 

I  [COEFFICIENT 

I 

I 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

0 


3 

3 


ALPH?=FPl2D*SSlN2(J)/2 
C=COS (ALPH2) 

S=SIN(ALPH2) 

V=STRF(J)/CMPLXf(1+SV32(J)*SV21(J))*C«( 
SV32(J)*SV21  (J>-1  )  *  S  )  *  *2 

V  =  -V*SGI (J ) *  A  I M 
AR(9)=V 
AR (10) =-V 
AR  <1 1  )=-V 
AR (12) =V 

V=CVPLX (ST( J) r  ALPHO*SLS(J) ) 
PH(9)=PI2*V 

PH(1 0) =P 12* ( V  +  SDTS  (  J)  ) 

PH (1 1 ) =PI2* ( V+SDTR  (  J)  ) 
PH(12)=PI2*(V*SDTS(J)*SDTR(J)) 


I 

1 C  COMPUTE  BASEMENT  RAY  TRANSMISSI3 

1  CN  COEFFICIENT  3 

I 

I 


7 

1 

1 

1 

1 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 
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A  Z 

A  :  ALPH2=FRI2d*BSSIN2(J>/2  :  Z 

A  :  C=C0S(ALPH2)  :  Z 

A  :  S  *S I N ( A  LPH2 )  :  Z 

A  :  VsB5TRF(J)/CMPLX((1+BSV32(J)*BSV21(J))*C*  :  Z 

A  :  (eSV32 ( J >*BSV21 <J )-1 ) *S> **2  :  Z 

A  :  V  =  V*BS6I ( J  )*REF  :  Z 

A  :  AR(13)=V  :  Z 

A  :  AR ( 1 4 ) *-V  :  Z 

A  :  AR  ( 1 5 )  *  -V  :  Z 

A  :  A  R  ( 1 6 )  =V  :  7 

A  :  V=CMPLX(PST(J) ,ALPHO*BSLS(J))  :  Z 

A  :  PH(13)*PI2*V  :  Z 

A  :  PH(14)=Rl2*(V*bSDTS(J))  :  Z 

A  :  PN(15)=PI2*(V*BSDTR(J>)  :  7 

A  :  PH(16)=PI2*fV*BSDTS(J)+BSPTR(J))  :  Z 

A  :  V=0  :  Z 

A  Z 

A  I  Z 

A  I  Z 

A  Z 

A  P - >  :  DO  210  K  B 1  *  1  5  :  Z 

A  P  Z 

A  P  I  Z 

A  B  I  Z 

A  S  Z 

A  P  :  N  A  =  K  + 1  :  Z 

A  p  7 

A  P  I  Z 

A  P  I  Z 

A  P  Z 

A  P  C - >:  DO  21G  L x  N  A , 1 6  :  Z 

A  B  C  . . Z 

A  P  C  I  Z 

ABC  I  Z 

a  e  c  .  z 

APC  :  AUX99=AR(K)*C0NJG(AR(L>)  :  Z 

A  P  C  Z 

APC  I  Z 

APC  I  Z 

APC  / . . . . . \  TRUF  Z 

APC  <  I  F  (  C  ABS  (  AUX90)  .EQ  .0.  )G0  TO  210  > - 0  Z 

APC  \ . /  X  Z 

APC  I  FALSE  X  Z 

APC  I  X  Z 

ARC  I  X  Z 

A  P  C  .  X  Z 

ABC  :  AUX1  =  (PH(K)-C0NJG ( P  H ( L ) )  )  *  A  I  M  :  X  Z 

A  p  C  . . .  x  z 

APC  I  X  Z 

APC  I  X  Z 

APC  / . . . . .  "  '  \  TRUE  X  Z 

APC  <  I  F  (REAL(  AUX1  )*  F?.LT  ,-40)GO  TO  210  > - V  7 

A  B  C  \ . /  X  Z 

APC  I  FALSE  X  z 

APC  1  X  Z 

APC  I  X  Z 
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A  P  C  . .  X  7 

ABC  |:  V=V*AUX99*(CEXPfF2*AUX1)-CEXP(Fl*AUXl))/  :  X  7 

ABC  j:  AUXl  :  X  7 

A  B  C  |: . .  x  7 

ABC  1  X  7 

ABC  0< - 0  7 

ABC  C  2 1  CD  I  7 

A  P  C  . 7 

A - :  CONTINUE  :  7 

A  7 

A  1  7 

A  1.7 

A  7 

A  :  AUX=2*RE Al<V)  :  7 

A  7 

A  1  7 

A  1  7 

A  / . . . . \  TRUE  7 

A  <  IF(CAPS(AP(T3)).EO.0.)6O  TO  209  > - 0  7 

A  \ . /  X  7 

A  I  FALSE  X  7 

A  1  X  7 

A  1  X  7 

A  X  7 

A  :  AUX2=-2* PI2*ALPhO*PSLS ( J )  :  X  7 

A  x  7 

A  I  X  7 

A  I  X  7 

A  / . . . ' . . . \  TPUE  X  7 

A  <  IF(  AUX2*F2.LT.-4u)G0  TO  2  09  > - V  7 

A  \ . /  X  7 

A  I  FALSE  X  7 

A  1  X  7 

A  I  X  7 

A  X  7 

A  :  AUX= AUX+4* ( EXF ( F2*AUX2 )-EXP( FI *AUX2) ) /  :  X  7 

A  :  AUX2*CABS(AR(13))**2  :  x  7 

A  x  7 

A  I  X  7 

A  0< - 0  7 

A  [  20QD  I  Z 

A  7 

A  :  AUX2=-2*PI2*ALFH0*6LS(J)  :  Z 

A  :  AUX2=-2*PI2*ALPhO*6LS(J>  :  Z 

A  Z 

A  I  Z 

A  I  Z 

A  / . . .  \  TRUE  Z 

A  <  I  F  (  AUX2*F  2.LT  .-40)GP  TO  211  > - 0  Z 

A  \ . /  X  Z 

A  I  FALSE  X  Z 

A  I  X  7 

A  I  X  7 

A  / . '\  TRUE  X  7 

A  <  IF(CABS(ARC5)).EQ.0.)G0  TO  211  > - V  7 

A  \ . /  X  7 

A  I  FALSE  X  7 

A  I  X  7 

A  I  X  7 
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A  X  7 

A  I:  AUX'AUX-UMEXF^  F2*AUX2)-EXP(F1*AUX2)  )/  :  X  7 

A  p  AUX2*CABS (AR (5) )**2  :  X  7 

A  . . x  7 

A  I  X  7 

A  0< - 0  7 

A  c  2113  I  7 

A  .  7 

A  :  AUX2=-2*PI2*ALPH0*SLS(J)  :  7 

A  . . .  7 

A  I  Z 

A  I  Z 

A  t"  . . . .  . *~~\  TRUE  Z 

A  <  IF(AUX2*F2.LT. -40)60  TO  212  > - - - 0  7 

A  > . /  X  Z 

A  I  FALSE  X  7 

A  I  X  7 

A  I  X  7 

A  f  . \  TRUE  X  7 

A  <  I  F  (  C  ABS  (  AR  (9)  )  .EQ  .0.  )60  TO  212  > - - - V  7 

A  \ . /  X  7 

A  I  FALSE  X  7 

A  I  X  7 

A  1  X  7 

A  .  X  7 

A  :  AUX  =  AUX  +  4* ( F XP ( F2* AUX2 )-E XP ( FI  * AUX2 ) ) /  :  X  7 

A  :  AUX?*CABS (AP (9) )**2  :  X  ‘  7 

A  X  7 

A  I  X  7 

A  0<  - . — . . 0  7 

A  C  2  1  2  D  I  7 

A  Z 

A  :  AUX=AUX/F1*4*CARS(AR(1  )  )  *  *  2  :  7 

A  Z 

A  I  Z 

A  I  Z 

A  t" . .  . . . . . . . .  "  \  TRUE  7 

A  <  I F  ( AUX  -GT  .0.  )  GO  TO  213  > - 0  7 

A  \ . t  X  7 

A  1  FALSE  X  7 

A  I  X  7 

A  I  X  7 

A  .  X  7 

A  :  ALOSS ( J J , J  )  =  2 5 .  :  X  7 

A  . . .  x  7 

A  I  X  7 

A  I  X  7 

A  / . \  X  7 

A  :  60  TO  20P  : - 0  X  7 

A  \ . /  U  X  7 

A  U  X  7 

A  0X7 

A  0< - 0  7 

A  C  2 1 3  D  I  U  7 
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A  .  U  7 

A  :  ALOSS(JJ,J)=-10.*ALOG10(AUX/(2.*RGI (J))**  :  u  z 

A  :  2)  :  U  7 

A  . . .  u  7 

A  I  U  7 

A  0< - 0  7 

A  C  20C  3  I  7 

A  7 

- -  CONTINUE  :  7 

. . 7 

I  7 

C  7003  I  7 

•••••••••••••••••••••••••••••••••••••••••••••  7 

:  N=J-1  :  7 

:  NA=j  :  7 

. . 7 

I  7 

I  7 

.  7 

A - > :  DO  701  1=1, N  :  7 

A  7 

A  I  7 

A  I  7 

A  7 

A  :  J  =  N  A  - 1  :  7 

a  7 

A  I  7 

A  I  7 

A  / . . . \  TPUF  7 

A  <  IF(J,LT,NSAP)GO  TO  909  > - 0  7 

A  \ . /  V  7 

A  I  FALSE  X  7 

A  I  X  7 

A  I  X 

A  X  7 

A  :  DR=PI2*RDTP(J)  :  X  7 

A  :  DS =P I2*P DTS  (  J  )  :  X  7 

A  :  AUX=DR*DS  :  X 

A  :  DT  =  S I N ( F 2* AUX ) -SIN ( FI  * AUX )  :  X  7 

A.  :  DT  =  DT/aux  :  x  7 

A  :  AUX=DS-DR  :  x  7 

A  :  DT=DT+SIN(F2*AUX)/AUX-SIN(F1*AUX)/AUX  :  X  7 

A  :DT»CT-2*(5INCF2*DS)-SIN(F1*DS))/DS  :  X  7 

A  :  DT=CT-2* ( SI K ( F2 *DP)-S IN( F 1* OR) ) / D°  :  X  7 

A  :  DT=  DT 121  H ♦ 1  :  X  7 

A  :  ALCSS(JJ,J)=-10.*ALOG10(DT)  :  X  7 

A  X  x 

A  I  X  7 

A  [  7013  I  X  7 

A  X  7 

- -  CONTINUE  :  X  7 

. . x  7 

I  x  ? 

I  x  7 

/ . . . \  X  7 

:  GO  TO  99°  : - V  7 

\ . /  x  7 
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A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


0< - 

c 

eOD3  I 

N  =  J-1 

NA  =  J 

I 

I 

• 

• 

DO  801  1=1 ,N 

• 

• 

• 

I 

I 

• 

J  =  NA-  1 

• 

• 

^  • 

I 

I 

- \  TRUE 

/  ' ‘ 

< 

I F  (J.LT.NTIR)GO  TO  999 

> - 

\ . , 

I  FALSE 

I 

I - C  COMPUTE  BIG 

I  [COEFFICIENT 

TRANSMISSION 

I 

I 


AlOSSCJJ , J )  =  25 • 
ALPH2=FPI2D*b$IN2 ( J) 
C=COS(AtPH2> 

$  =S  I  N  (  A  l  PH  2  ) 

V=BTPF(J)/C*PLX((1+FV32(J)*BV21(J))*C»( 
BV32(J)*BV21 <J)-1>*$>**2 

V  =  V  *BG I  CJ ) 

AUX3=CAPS ( V ) ♦ ♦ 2 


I 

I 

/ . . . . . . \  TRUE 

<  I  F  (  AUX ? • EQ • 0 ) GO  TO  801  > - 

\ . / 

I  FALSE 

I 

I 


DR=PI2*PD7R(J) 
D$=PI2*BDTS ( J) 

AUX=DP+LS 

DT=SIN(F2*AUX)-SIN(F1*AUX) 
DT=DT /AUX 
AUX=D$-DP 

DT=DT+SIN(F?*AUX)/AUX-SIK(F1*AUX)/AUX 

DT=rT-?*(SIN(F2*DS>-SlN(Fl*DS))/DS 

DT=rT-2*(SIN(F2#DP)“SlN(Fl*DR))/D° 

DTcDT/2/FUI 

AUx2=-2*PJ2*ALPHj*BLS(J)*EX 


167 


A  I  X  7 

A  I  X  7 

A  / . .  . . . . . . TRUE  X  7 

A  <  IF  <  AUX?*F  ?.LT  .-40.)G0  TO  801  > - V 

A  V . /  X  7 

A  I  FALSE  X  7 

A  I  X  7 

A  I  X  7 

A  .  X  7 

A  :  ALOSS( J J . J > =-10*< AUX2+ ALOGIO ( AUX3*DT/RGI (  :  X  7 

A  :  J)**2>)  :  X  7 

A  . . x  7 

A  I  X  7 

A  I  X  7 

A  / . . . . . . . . . \  FALSE  X  7 

A  <  IF(AL0SS(JJ,J).GT.?5.)  > - 1  X  7 

A  V . /  I  X  7 

A  I  TRUE  I  X  Z 

A  I  1X7 

A  .  I  X  7 

A  :  ALOSS ( J J  »  J ) =25  :  I  X  Z 

A  . .  I  x  z 

A  I  1X7 

A  0< . 0 

A  C  8C 1 3  I  X 

A  .  X 

:  CONTINUE  :  X 

. . x 

I  X 

0< - o 

C  C99D  I 


RETURN 


•  n  L  I  U  'I  • 

\ . / 
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SUBROUTINE  GCOMP 


* 

A 

A 

A 


\ 

\ 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


/ 

/ 

— C  THIS  SUBROUTINE  COMPUTES  THE  1 
[GRADIENT  (G)  ? 

[  AND  ITERATES  TO  FIND  THE  ACTUAL? 


[CRITICAL  ? 
[  ANGLE  (THC).  BOTH  ARE  THEN  ? 
[RETURNED  TO  THE  ? 
[  MAIN  PROGRAM.  ? 


COMMON  ZB,ZR,ZS»C1,C2,C3,RH1,RH2,RH3,G,D, 
ALPHO,CS,BZ,REF 

COMMON  /DATA/LZERO,FZERO,L70tF70,TH70» 
L71.F71.D PL OSS, FREQ, THETA,  ALOSSO.ALOSS*, 
FM,TI«>fTHC,THRC,TH1B,BfcTA 

REAL  LZEPO,L70,L71 
PI=3. 1  4  1  5  9265 
THRC  =3  0. 

THC0=THC 
B=  P  E  T  A 

THRC=THRC*(PI/1PO.) 

THINC  =  1  . 

M  =2*ZB~ZR~ZS 

TH3=ACOS((C?/C1)*COS(THRC)) 

F  A  C  T 1  =  C2.*C3*6)/(1 .♦?> 

SIN3=SIN ( TH  3 ) 

COS3=COS (TH3) 

SIN7SG=SIN3*S IN3 
C0S‘rSQ  =  C0S3*CPS3 

FACT?=1.+(1./(2.*B))*(1.+(TH3/(SIN3*COS3) 

)  ) 

GZEPO=(FACTl*(1./COS3Sj)*FACT?+<r*/(1.+P) 
)  *(1./COS3SQ)*(1.-((Th3*(COS?SQ-SIN3SG)) 
/(SIN3*COS3>>>>  *((C3/C1 )*(SIN(THCC)/ 
SIN3))/(D1/SIN(THRC)**2) 

R  =  (M/TAN(THRC))+(FACT1*TAN(TH3)/GZERC)* 

F  A  C  T  2 


THC1=ATAN(D1 /P) 

TH  C 1 =THC1 * (1 60. /PI ) 

MA  X I T  =  1  60 

I 

[  ico? 

I 

* 

DC  600  I  T  =  1  ,MA  X  I  T 

: 

I 

I 
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THRC=THRC*(180./PI) 

THR  C=THRC  +  THINC 
THRC=THRC*(PI/180.) 

TH3  =  A  COS  (  <C  3/Cl  )*COS(THPO) 
FACT1=(2.*C3*B)/(1.*B) 

SIN3=SIN( TH3) 

CO  S3=C  OS  ( TH3) 
SIN3SQ*S1N3*SIN3 
C0S3SQ=C0S3*C0S3 

FACT2=1.+(1./(2.*B))*(1.+(TH3/(SIN3*COS3) 
)  ) 

GZEPO=(FACTl*(1./COS3SG)*FACT?+(C3/(1.+e> 
)  *(1./COS3SQ)*f1.-((TH3*(COS3SQ>SIN3SQ)) 
/(SIN3*COS3)>>>  *((C3/C1)*(SIN(THPC>/ 
SIN3))/(Dl/SIN(THRC>**2) 

R=(Dl/TAN(THRC) )  ♦  (  F  A  C  T  1  *  T A  N  (  TH3  )  /  G  Z E  RO) * 

F  ACT2 

THC  ?  =  AT AN U1  /R) 

THC2  =  THC2* (1 80.  /PI  ) 


1 

I 

/ . ******************'******* . -«***<r**y  TRUE 

<  IF (APS (THC0-THC1 ) »LE •  1.E-A)  GO  TO  7C0  > - 

\ . / 

I  FALSE 

I 

I 

/ . '\  TRUE 

<  I F ( ( (THC0-THC1 ) *(THC0-THC?) )  .GE.  0.)  GO  > - 

<  TO  EDO  > 

\ . / 

I  FALSE 

I 

I 


: 

THINC=  -THINC/2. 

: 

I 

0< - 

c 

5003 

I 

* 

THC1 =THC2 

• 

» 

I 

c 

6003 

I 

: 

CONT 1 NUE 

• 

• 

• 

• 

I 

I 

/"' . . . . . / 

/  W  F  I T  E  ( 1  *6  50)  l*AXIT,THC0,THC2  / 

. / 


/ 


z 

I  z 

I - C  650  F0R»AT(1X,I5»2XfnlTERATI0NS3  Z 

I  [AND  STILL  NO  SOLUTIONo,/,  1  1Xt«>3  Z 

I  CTHC0*ntFl0.6,5X,DTHC2*D,F10.6)  3  Z 

I  Z 

I  Z 

/ . . . \  Z 

:  GO  TO  800  : - 0  Z 

\ . /  X  z 

x  z 

X  z 

0< - .0 

C  7003  I  X 

.  X 

THC=THC2  :  X 

G  =G  Z  E  R  0  :  X 

. X 

I  X 

o< - 6 

[  8003  I 

:  RETURN  : 


\ . / 
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/ . \ 

:  BEGIN  : 

\ . / 

I 

I 


BLOCK  DATA  NULL 

COMMON  /DATA/ALZERO, F Z ERO , A L 70 . F 70 , TH70 , 
AL71,F71,DBL0SS,FREQ.  TH E T A , ALOS SO , 
ALOSS*,FM,TIR,T«C,THRC,THlB,BETA 
CCPPON  /ALPH/GAM 
COMMON  /T/T 

COMMON  /R/NSAP,RV21(90),RV32(90),RSIN2( 
90)  ,PDTS (90) ,RDTR (90) ,  RT (90) ,RGI (90) 
COUPON  /B/BV21 (90) «PV32(90) »ESIN2(90) « 
BTRF (90) ,PGI (90 ) ,BLS (9 j) ,  B T ( 90 ) , B D T S ( 90 ) 
,BDTR (90) 

COP PON  /S/SV21 (90) «SV32(90) »SSIN2(90) « 
STRF(9C),SGI(90),SLS(9C),  S T  ( 90 )  ,  S DT S ( 90 ) 
,  SDTP (90) 

COUPON  /B$/PSV2T(90),BSV32(90),PSSIN2(90) 
, BSTRF(90) .BSGI (90) .  BSLS(90),BST(90), 
BSDTS (90)  ,BSDTR  (90) 

COPMON/AL/F (6)  ,  AL OSS (6,90) 

DATA  F/50. ,  IOC.,  2 00., 4 00., 800., 1600./ 


I 

I 

/ . \ 

END 

\ . / 
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APPENDIX  C 


RUN  STREAM 

f 

An  example  of  the  computer  input  run  stream  for  each  of  the  two  input  options 
follows . 

:RUN,BTLS9 

INPUT  OPTIONS  ARE:  0  =  INPUT  GEO-PARAMETERS 

1  =  INPUT  BOTTOM  LOSS 

2  =  END  PROGRAM 

ENTER  STATION  NUMBER  AND  INPUT  OPTION: 

8  0 

ENTER  THE  TWO-WAY  TRAVEL  TIME  THICKNESS  IN  10THS  OF  SECONDS: 

4.58 

ENTER  THE  WATER  DEPTH: 

1800.0 

ENTER  SOURCE  DEPTH  AND  RECEIVER  DEPTH: 

244.0  305.0 

ENTER  SOUND  SPEED  AT  SOURCE  OR  RECEIVER  -  WHICHEVER  IS  SHALLOWER: 

1511.72 

ENTER  SOUND  SPEED  AT  THE  BOTTOM  OF  THE  WATER: 

1535.11 

ENTER  RATIO  OF  SPEED  IN  SEDIMENT  TO  SPEED  AT  BOTTOM  OF  WATER: 

.996 

ENTER  GRADIENT  AT  TOP  OF  SEDIMENT  AND  CANONICAL  CURVE  TYPE: 

1.7  .86 

ENTER  THIN  LAYER  DENSITY,  SEDIMENT  DENSITY  AND  THIN  LAYER  THICKNESS  IN  METERS: 
2.33  2.33  .04 
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ENTER  SURFACE  ATTENUATION  IN  DB/M/KHZ: 

.021 

ENTER  ATTENUATION  PROFILE  GRADIENT  AND  BASEMENT  REFLECTION  ^COEFFICIENT  IN 
PRESSURE: 

.00015  .708 

RUN  ,BTLS9 

INPUT  OPTIONS  ARE:  0  =  INPUT  GEO-PARAMETERS 

1  =  INPUT  BOTTOM  LOSS 

2  =  END  PROGRAM 

ENTER  STATION  NUMBER  AND  INPUT  OPTION: 

8  1 

ENTER  THE  TWO-WAY  TRAVEL  TIME  THICKNESS  IN  10THS  OF  SECONDS: 

4.58 

ENTER  THE  WATER  DEPTH: 

1800.0 

'  ENTER  SOURCE  DEPTH  AND  RECEIVER  DEPTH: 

244.0  305.0 

ENTER  SOUND  SPEED  AT  SOURCE  OR  RECEIVER  -  WHICHEVER  IS  SHALLOWER: 

1511.72 

ENTER  SOUND  SPEED  AT  THE  BOTTOM  OF  THE  WATER: 

1535.11 

ATTENUATION  CALCULATION 

FOR  A  POINT  BELOW  THE  CAUSTIC  WHERE  BL  CURVE  IS  NEARLY  LINEAR, 

ENTER  A  LOSS,  FREQUENCY  AND  ANGLE:  (ANG=APPROX  20,  FREQ  =  APPROX  200) 

5  20  200 

ENTER  LOW  FREQ  LOSS  LEVEL  AT  HIGH  ANGLES, 
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*  THE  MINIMUM  LOSS  LEVEL  AT  HIGH  ANGLES 
AND  THE  FREQUENCY  OF  THE  MINIMUM  LOSS: 

10  6  80 

j 

ENTER  ANGLE  OF  TOTAL  INTERNAL  REFLECTION, 

ANGLE  OF  CAUSTIC, 

RATIO  OF  SPEED  IN  SEDIMENT  TO  SPEED  AT  BOTTOM  OF  WATER:  ' 
55  34  .996 

ENTER  CANONICAL  CURVE  TYPE  AND  ATTENUATION  PROFILE  GRADIENT 
.86  1.7 

ENTER  BASEMENT  REFLECTION  COEFFICIENT  IN  PRESSURE  UNITS: 
.708 

/BTL14  *RUNTIME  ERROR*  22UN  @  40126 
/BTL14  *RUNTIME  ERROR*  22UN  @  4022 


.175 


176 


rpEOuEmjv  •  m  hz 


m  o 


o  r 


*J  O  t/>  l/i 


177 


CMZPIG  AWCIE 


FPEOtlEMCV  •  8W  HZ 


•o  O  #-  *-  O  C 


©  m 


GRAZING  ANGLE 


rPEOUtHtV  *  HZ 


r 
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GBAZruG  AWGiE 


FREOU€NCV  •  800  HZ 


180 


CRAZING  ANGLE 


FffEOUCWV  ■  16M  MZ 


r 


* 

w 


m  ©  *-  o 
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G*A2JHC  AMCU 


r.ie 

INPUT  :pt:oh:*  APE  :  0  •  INPIT  GE0-PARAF£T£9$ 

;  •  INPUT  BCTTC r  LOSS 
2  •  END  PRCGRAI" 

•«"IP  -TATKN  hUfliER  AND  INPUT  OPTION : 

:•  o 

ErrEP  ^H£  TUm-Um>  TPauEU  T ITE  ThIC<NESS  In  ie*HS  0r  SECONDS: 

J.Sc 

EN*EP  "“HE  UA^ER  DEP^H: 

1  nV 

ErTER  SOURCE  DEPTH  ANC  9ECEIUE9  D£PTH •* 

BA4  30S 

ENTEP  SC-LNr  SPEED  at  SOURCE  09  RECEIUER  -  UHlCW£U£P  1$  SHAUGLER: 

1C1  ‘ 

EN"EP  SOUND  SPEED  AT  THE  BOTTON  Of  TH£  LATER: 

1S3S.S 

ENTER  RATIO  OP  SDEED  IN  SEDIflEN^  TO  SPEED  AT  BOTTon  OT  LATER: 

.9% 

ENTER  GRADIENT  AT  TDp  Of  SEDIP1ENT  AND  CANONICAL  CUPUE  TYPE: 

1  . 97  .86 

ENTER  THIN  IAVE9  DENSITY,  S£ D I N£ NT  DENSITY  AND  THIN  LAYER  THICKNESS  IN  PETE9S : 

2.33  2.33  .04 

EN^ER  SURFACE  ATTENUATION  IN  DB/P/KNZ : 

.021 

EN~ER  ATTENUATION  PPOfILE  GRADIENT  AND  BASENENT  R£f LECTION  COEf f ICIENT  IN  PRESSURE 

.eeeis  .?es 


182 


OUTPUT  OF  COMPUTER  RUN 


CRITICAL  ftHGirES 


-996-0-1HC-- - 34,-74-  WR€- 


TIR  — - 6 , 00— R8T-  - 

THB  =  30.6?  THRB  =  45.0? 

TWT-T— = - ,46  6Z--= 422.  <3  CRD  = 47-7.^5-?- 


—46  .-97 - 


G500  = 


1 .34  C500  =  2279.83  A500  = 


. 096  GAM  = 


. 00015 


PAT  = 
-RW2— 
D  = 


.9960 
-2,33  -RH3-= 


2^33 


04 


U  — S - 

AL  0= 

+  Ut.  1  - 

. 02 1 GAM= 

.  00015 

.  p?  =  422 

.13-  REF= 

.71 

-1.25 

-  .32 

.  05 

-.11 

-  .86 

-1  .65 

-1  .68 

-1  .  79 

_ -g . Q£ 

_  -1 .97 

-1 .81 _ 

-1-45 

- .  fit 

.  15 

2.22 

_ 2 .2 1 

1  .42 

t  .29 

1  .47 

1.19 

.3? 

1  .  09 

1  .71 

1  .42 

_ _ _ _ l  .24 

I  .3i  .. 

.77  .. 

__1  ,*U _ 

_  2,03  „ 

_  _8 ,07_ 

_ 2j  35 

7.07 

6.95 

6 . 75 

6 . 74 

6 . 67 

6.54 

6.35 

6.31 

A  .  47 

A  313  . 

6.10 

6,12 

5.99 

5.82 

5.97 

6 . 09 

5.98 

5.89 

5 .82 

5.82 

5.83 

5 . 86 

5.79 

5.76 

_ 5.75  _ 

_ 5.69 

5  .  S3 _ 

_ 5 . 56 

5 . 57 

5 . 48 

_ __5,39  . 

5.40 

5.43 

5.45 

5.4? 

5.45 

5.46 

5.49 

5 . 5  0 

_ _ _ 5 .50 

_ 5^-SJ _ 

_  _5_.52_  - 

.  5.52 _ 

_ 5J£„ 

_ 5_52 _ 

.5.52 

_ 5,52 

5.51 

5.49 

5 . 48 

'■D 

in 

5.45 

5 . 45 

5.45 

5.45 

4.78 

7,87- 

'-1  .47 

-.14 

-.10 

-.97 

-  .83 

-  .51 

-  .85 

-  .91 

-  .  39 

■  54 

_ 1  .81 

_ 2,65 

_ 2*0.0,. _ 

„1*41 

.  J  ,25 

_ 1,  66 

1  .41 

1  .30 

1  .96 

2.71 

1  .52 

2.57 

2.26 

2.36 

2 . 6S 

_ 2*80 

3.39 

3.68 

4 . 07 

4.35 

7.93 

6,13 

7.99 

7.94 

7.89 

7.69 

7.55 

7,27 

7.24 

7.23 

7.16 

7.16 

7 . 06 

7. 00 

6.85 

6.97 

6 . 89 

6 . 87 

6 . 85 

6.91 

6.75 

6.72 

6.57 

6.57 

6.58 

6 . 6i 

6.62  .. 

6 .58 

.  6.58  . 

6.5S 

£  *50- 

_ £  *53 

_ 6,62 

6.57 

6.57 

6 . 55 

6.56 

6 . 56 

6.56 

6.56 

6.53 

_ 6-43 

5-42 

_ 6.43 

_  .6,41  .. 

6,35 _ 

.6  ..32  _ 

_ 6,33 

6.34 

6 . 36 

6.36 

6.35 

6.35 

6 . 35 

6.37 

6.37 

6.36 

5.2ft 

7 . 95 

1  .64 

1  .27 

-  .  05 

-.20 

-  .  46 

-  .  32 

.  07 

1  .  74 

_  3-44 _ 

3.30 

_ 2. 16 

_ 1_*88 

...  2.20 _ _ 

2^69  . 

2.23 

2.46 

3.00 

2.99 

3.13 

4 . 45 

4.38 

3.74 

4.26 

4.71 

_ 5  .  <11 

_ 5*58 

6 . 02 

_ £.41_ 

_  7.02 

7.55 

8.11 

_ 8 , 1  4 

8.17 

7.95 

8 . 00 

8.06 

7.95 

7.98 

7.99 

7.90 

7.61 

7.82 

7.91 

7.76 

7.80 

7,82 

7.74 

7.62 

7.65 

7.62 

7.65 

7.56 

7.62 

7.65 

7 . 66 

7 .58 

_ 7*5 1 _ 

_ 7.52 

7.50 

_ 7,57 _ 

7.50 _ 

7 . 49 

7.50 

_ 7 , 5  0 

7.51 

7.45 

7 . 45 

7 . 39 

7 . 39 

7.41 

7.41 

7.40 

_  7*48  . 

7,42 

7,41 

7^36 

7,42 

7,42 

7,  4J 

_ Z,_4 1 

7.42 

7 . 43 

7 . 44 

7.43 

7 . 43 

7 . 44 

7 , 44 

7 . 45 

f*  .  OS 

8 . 08 

*  1.43 

1  .39 

2.27 

1  .  75 

2.38 

O  .  b  i' 

4 . 64 

2 . 78 

..  2.64 

2 . 79 

3.62 

.7.66 

3*  S£ 

4,  Ob 

4 . 73 

4  .79 

4 . 94 

5.45 

5 . 76 

6 . 77 

6 . 06 

6.60 

6 . 84 

7  .  1  6 

7.46 

7.64 

7 . 88 

8 . 03 

8  ■  03 

8.17 

8,13 

8 . 07  .  

8  .  OS 

8.10 

8 . 08 

8 . 06 

8 . 04 

8.12 

8. 03 

8 . 05 
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- 0-rft? — 

8 . 05 

-3  qr. 

— 6-r-&fi - 

7.97 
y  e,f 

8 , 02 — 

8. 03 

?  Ch 1  . 

— 8,02 

8. 02 
?  ClK, 

- 6-r47 - 

8. 03 

->  a  a 

-7,99 — 

7.95 

7  h£ 

— 6-1—86 - 

7.96 

7  on 

8.  03 

7.96 

7  A  f 

i  t  J  J 

7.93 

-I  A  4 

f  1  70 

7.97 

Q  O 

r  *  7*r 

7.91 
y  aq 

T  t 

7.95 

7  97 

f  i  7*t 

7,96 

_ 2^a.-3- 

r  . 

7.93 

*»  Q, 

f  »  77 

7.90 

7  Q*s 

r  t  9  1  - 

7 . 92 

_ _ 7  Q  4 

/  <  7  1 

7.90 

6  97 

r  i 

7.88 

8  00 

7,89 

r  i  7  a 

7.90 

r  i  7w 

7.90 

r  f  7T1 

7.90 

7,90 

■  — r  i  7  i - 

7.90 

3.44 

-&r+4. 

3,01 

5-r33- 

3. 18 

4. 06 

3.94 

3 , 73 

4 . 08 

-  4,83  — 

5,13 
—  -7-— 07 

7,37  7.55  7,60  7,52  7,92  8,04  8,11  8,16 

8-r44 - 8, 4-6 — — &T46 - 8^+2 - fti-13 - 8^4-0 - 8.44 - ft^t-O  — 

8.08  8.11  8.08  8,05  8,08  8.06  8,09  8.06 


- ¥ i - 

- <7  T  -VO - 

— g  >  yj - 

- ¥  t  - 

- - 

- - ¥  *  - - 

- g ■ - 

- ¥  ,  47^ - 

8 , 02 

8 , 05 

8 , 03 

8, 04 

8.02 

8. 03 

8. 01 

8. 02 

-£.-01  - 


8, 02 

8.00 

-Z-7£ 

4 , 49 

8,31 

8,18 


8.01 

8 , 02 

8. 03 

8  ,  01 

8. 01 

8. 02 

8 . 01 

8  02 

8 , 02 
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